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AN APPARATUS FOR THE STUDY OF . 
SULTRASONIC PROPAGATION IN LIQUIDS UNDER PRESSURE 
by R. I. Tarr 

King’s College, Newcastle upon Tyne, England 


A description is given of an apparatus designed for the measurement of velocity and 


absorption of glee scontcs in liquids. The apparatus has a capacity of two litres and can 
coyer the range 0 to 703 kg/cm? (0 to 10000 lb/sq. in.). A pulse method is used for both 


velocity and attenuation and either fixed or variable path-length can be used. 
Measurements in water at temperatures 15° C to 80°C are presented and compared 
if with recent measurements by others and with theory. 


Sommaire 


On décrit un appareil imaginé en vue de la mesure de la vitesse et de l’absorption des 


" ultrasons dans les liquides. Cet appareil a une capacité de deux litres et permet de faire 
\ des mesures entre 0 et 703 kg/cm? (0 et 10000 livres/pouce?). On utilise une méthode 


étre fixe ou variable. 


dimpulsions pour déterminer la vitesse et l’affaiblissement; la longueur de trajet peut 


On donne les résultats de mesures faites dans l’eau, aux températures 15° 4 80° C, et 
p 9 


on les compare ayec des résultats récents d’autres auteurs ainsi qu’avec la théorie. 


Zusammenfassung’ 


Es wird eine zur Messung der Geschwindigkeit und Absorption von Ultraschall in Fliis- 


sigkeiten entwickelte Apparatur beschrieben, die 2 Liter Fliissigkeit aufnehmen kann und 
4 in einem Druckbereich von 0 bis 775 kp/cm? (0 bis 10 000 Ib/sq. in.) verwendbar ist. Bei 


festgehaltener oder veranderlicher Lange des Mefiweges kénnen sowohl die Schallgeschwin- 
digkeit als auch die Dampfung mit Hilfe einer Ta puiceisiinds gemessen werden. 


MeB8ergebnisse fiir Wasser im Temperaturbereich von 15° bis 80° C werden angegeben 
und mit der Theorie und mit den in letzter Zeit durchgefiihrten Messungen anderer 


Autoren verglichen. 


1. Introduction 

Although a great deal of research on the pro- 
pagation of sound in fluids has been carried out in 
the past two decades, comparatively little attention 

s been paid to the study of sound velocity and 
nore particularly of sound absorption at high pres- 
sures. Apart from the early work of Swanson [1] 
ind Biauarp [2], the most noteworthy contributions 
in this field are those of Horron [3] and of Smiru 
and Lawson [4], who investigated the pressure de- 
pendence of velocity of sound in water up to pres- 
sures as high as 9000 kg/cm”. More recently Lrrovirz 
nc CarnevALe [5] reported both velocity and ab- 
orption measurements in water to a pressure of 
2000 kg/em?. 
_ The purpose of the present work is to study the 
ffect of pressure on velocity and absorption in other 
liquids besides water. Although the pressure range 
ised is appreciably smaller than that used by some 
of the investigators mentioned above, it is sufficient 
to establish velocity and attenuation isothermals 
from which valuable information may be deduced. 


- 


Measurements have been made in water, glycerine, 
benzene, carbon tetrachloride, pentane, ether, carbon 
dioxide (liquid phase) and mono-fluor trichlor- 
methane (freon 11). The results for water alone are 
presented in this paper; measurements on the other 
liquids are described elsewhere [6]. 

In some respects the velocity data of Smirn and 
Lawson differ from those of Hotron, the principal 
discrepancy involving the behaviour of the maximum 
in the velocity of sound as a function of temperature 
as the pressure is increased. In the present work 
particular stress was laid on accuracy of measure- 
ment and as the results obtained in general favour 
Smita and Lawson’s data it is considered of value 
to report them here. No other work on the effect of 
pressure on absorption in water appears to have 
been reported in the literature, apart from the paper 
by Lrrovirz and CarnevaLe whose measurements 
were made at temperatures of 0° and 30° C. The 
data given here is for 16° and 23° C. 

Both variable and fixed path acoustic interfero- 
meters were built and used with a pulse system for 
the measurement of velocity and absorption. Velo- 
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city measurements were made using a pulse inter- 
ference technique, a method whereby the received 
echo is mixed with continuous R.F. in such a way 
that the amplitude of the displayed echo depends on 
the relative phase of the two signals. In this manner 
the velocity can be measured to the same degree of 
accuracy as can be achieved using continuous waves 
in a conventional Pierce interferometer. As is shown 
below, the method can be used with a fixed path 
acoustic system and is in general superior to the 
usual method used by workers at high pressure, i. e. 
that of measuring directly the pulse transit time 
through a known path length. 


2. Apparatus 


2.1. Pressure system 


A schematic diagram of the pressure system is 
shown in Fig. 1. The apparatus was designed for a 
normal maximum working pressure of 10 000 Jb/sq. 
in., obtained with a hand operated screw intensifier. 


Pressure gauge 


Reservoir 
Water jacket 


Fig. 1. Schematic lay-out of apparatus. 


Early measurements were made with a moving path 
interferometer which is an integral part of the main 
pressure chamber shown in Fig. 2. The screw car- 
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Fig. 2. Main pressure chamber and variable-path inter- 
ferometer. 


rying the reflector passes right through the pressure 
cylinder, being sealed at either end with S.E.A. ring- 
fittings. A graduated head and revolution counter is 
attached to the upper end of the screw shaft and a 
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power drive is applied at the other end to control 
the reflector movement. The crystal holder is moun- 
ted on three spring-loaded screws which enable the 
necessary parallelism adjustments of crystal and 
reflector to be made. The cylinder end caps are 
secured at either end with large nuts and pressure 
obturation is achieved using dural wave rings. 
Provision is made in the upper end cap for the 
admission of an electric lead seal to enable R.F. to 
be fed to the crystal. For frequencies up to 6 Me/s 
Barium Titanate disc transducers were used. Owing 
to the low impedance of these crystals it is possible 
to obtain satisfactory results without any isolation 
of the live crystal lead from the liquid under investi- 
gation. Higher frequencies (above 6 Mc/s) demanded 
the use of high impedance transducers such as 
quartz, and for these crystals the insulation required 
is achieved by using an oil-filled crystal holder which 
was fitted with a diaphragm to allow pressure 
equalization. 


Although the measurements on water described 
here were made using the moving-path acoustic 
system, for most purposes a more convenient fixed- 
path interferometer was employed. This device con- 
sists simply of a mounted crystal and reflector sepa- 
rated an accurately known distance by fixed spacers, 
the whole being mounted in a tube sealed at one end 
with a sylphon bellows. The tube fits inside the 
pressure cylinder in place of the moving interfero- 
meter screw, and the cap holes are closed by screwed 
plugs. 


Comparative tests of the two types of interfero- 
meter showed that comparable accuracy in the 
measurement of velocity can be achieved and in 
many respects the fixed-path instrument is to be pre- 
ferred for the work in hand. The unit is easy to 
handle and may be quickly removed from the pres- 
sure cylinder for cleaning and refilling. The time 
taken in measurement is considerably reduced; more- 
over many more points on a velocity isotherm may be 
obtained than would in practice be possible with the 
variable-path instrument. It also has the great 
advantage of appreciably reducing the volume of 
liquid required for measurement, while at the same 
time isolating the liquid from that used only for 
production of pressure. 


The use of a fixed-path interferometer for ab- 
sorption measurements presents some difficulty. A 
technique sometimes used consists simply in obser- 
ving the change in the amplitude of the echoes as the 
pressure is varied. This method is however prone 
to large experimental errors. The arrangement finally 
adopted was to mount in the beam path an auxiliary 
reflector somewhat smaller than the main reflector. 
A proportion of the beam is thus intercepted to pro- 
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Air vent 
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Main reflector 
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Fig. 3. Fixed-path interferometer. 


2.2. Electrical apparatus 


The circuit used for velocity measurements is 
shown in Fig. 4. The conduction of a gating valve 
Vz, is controlled by an asymetrical multivibrator V, 
which allows short bursts of R.F. from the Hartley 
Oscillator V, to pass to the crystal. The received 
echo is amplified by V, and before detection is mixed 
with a pre-determined amount of R.F. from V,. The 
echo amplitude then displayed on the oscilloscope 
thus depends on the relative phase of the echo com- 
pared with the continuous R.F. When the circuit is 
used in conjunction with a moving path interfero- 
meter, the echo pulse will pass through a maximum 
and minimum for each half wavelength traversed by 
the reflector. The wavelength and hence the velocity 
can thus be determined. By counting a large number 
of echo maxima and monitoring the frequency by a 
crystal check wavemeter it is possible to measure the 
velocity to better than 0.1%. The process is repeated 
at different pressures and temperatures to give the 
required results. 

When used with the fixed-path interferometer the 
pressure is slowly varied and the number of pulse 
maxima plotted against pressure. Suppose n such 
maxima are counted in increasing the pressure from 
atmospheric, py, to a value p. Then the velocity 
v, may be shown to be related to the atmospheric 
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Fig. 4. Circuit for velocity measurement. 


vide a reference echo and measurements are made 
by comparing echo levels and noting the manner in 
which the difference in echo amplitudes changes with 
pressure. Fig. 3 shows the fixed path interferometer 
with the subsidiary reflector mounted on the crystal 
holder, an arrangement which permits both reflectors 
to be accurately aligned parallel to the transducer. 


pressure velocity vg by the relation 


- oi 1 

Vy 2fd (1) 
where d is the distance between crystal and reflector 
and f the frequency. The sign of n is dependent upon 
whether the velocity is increasing or decreasing: for 
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all liquids n is negative, i.e. the velocity increases 
with pressure. 

It is possible to measure vy with a fixed-path 
acoustic system by tuning the crystal through its 
resonant frequency from, say, /; to fg and observing 
the number of pulse maxima, m. vg is then given 
by the expression 

2d 


m 


lh—fel- (2) 


However only a small range of frequency can be 
covered within the resonant peak of the crystal and 
in general the experimental accuracy is inferior to 
that which may be obtained by more conventional 
methods. In this work a moving path interferometer 
was usually employed to establish vg when necessary. 


V9 


Electrical apparatus used for absorption measure- 
ments comprised a pulse generator and amplifier, 
calibrated attenuator, high gain receiver and an 
oscilloscope display unit, all of which have now 
become standard items of equipment for this type 
of investigation. In order to obtain greater accuracy 
in the monitoring of the received pulse amplitude a 
method of pulse suppression was devised whereby 
the top section only of the echo pulse envelope was 
displayed on the oscilloscope. This expedient enabled 
small changes in echo amplitude (and hence attenua- 
tion) of the order of 0.1 dB to be readily discerned, 
and it was of particular value when used with the 
fixed path interferometer where the amplitude chan- 
ges of the main and reference echoes due to an 
applied pressure were relatively small. 


R.I. TAIT: ULTRASONIC PROPAGATION IN LIQUIDS 


‘ACUSTICA 
Vol. 7 (1957) 


For the control of temperature in the apparatus a 
“closed circuit” system is used whereby water is 
circulated from a small auxiliary heating tank to the — 
cylinder water jacket (see Figs. 2 and 5). The heat- — 
ing tank contains a double immersion heater and 
thermostat by which the temperature may be set to 
any value from 16° to 95° C. Two thermometers — 
were employed to measure the water temperature; 
one in the heater tank and the other in a sealed 
pocket built into the water jacket. In addition — 
a thermocouple was mounted inside the upper part 
of the pressure cylinder in close proximity to the 
crystal holder. With the pump in operation the ther- 
mostat gave a control of temperature inside the 


é oats 1 : 
heating tank to within +,°C over long periods, 
and owing to the great thermal capacity of the pres- 
sure cylinder this variation was reduced to less than 


ac z °C in the cylinder itself. 


Fig. 5 shows the complete equipment as used for 
velocity and absorption measurements. The pressure 
apparatus was solidly mounted on wooden plinths 
resting on a concrete floor with the electronic equip- 
ment arranged conveniently behind. Provision was 
made for the pressure cylinder to be used in either 
a horizontal or vertical position: it is shown in the 
vertical with the variable-path interferometer in 
place. The drive motor and gearbox may be seen 
under the cylinder. Coarse and fine speed controls 
together with a reversing switch are mounted on the 
control panel seen at the left of the figure. The 
heating tank is to the left of the oscilloscope. 
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Fig. 5. Photograph of complete apparatus. 
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3. Results 
3.1. Velocity 


The measurements of ultrasonic velocity in air- 
free distilled water are plotted as isothermals on 
Fig. 6. A frequency of 3 Mc/s was used. Isobaric 
curves are given in Fig. 7: it will be seen that the 
shape of the isobars in the neighbourhood of the 
-yelocity maximum depends critically on the slope of 
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Fig. 6. Velocity isothermals; water, frequency 3 Mc/s. 
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Fig. 7. Velocity isobars; water, frequency 3 Me/s. 


the isotherms, particularly those at the three highest 
temperatures. The experimental error involved in the 
temperature range was estimated at 1m/s and 
1.5 m/s at zero and maximum applied pressure re- 
spectively. The curves of Fig. 6 are linear to within 
these limits and their slopes show a slight increase 
with temperature; an effect more pronounced at the 
higher temperature resulting in the two curves at 
68.5° and 78° C overlapping. The extent of overlap 
lies barely outside the experimental error but re- 
peated measurements did justify the results given. 
The direct inference from this is immediately shown 
in Fig. 7 where the temperature of maximum velocity 
is displaced with pressure to a higher temperature 
and in fact at 78°C a maximum is not reached 


above about 3000 Ib/sq”. 


3.2. Absorption 


Measurements of absorption coefficient (a) were 
carried out in distilled air-free water using the 
variable path interferometer at a frequency of 
12 Mc/s. The results obtained at two temperatures 
16° and 23° C are given in Fig. 8 in the usual form 
of a/f? vs. pressure. The value of a//? obtained for 
atmospheric pressure are in reasonable agreement 
with published results for absorption in water. The 
experimental error was estimated at 3% for the 
lower temperature and 5% for the higher and within 
these limits a/f? exhibits a linear decrease with pres- 
sure, the 23°C curve having a somewhat lower 
slope. 


Applied pressure ——= 


Fig. 8. Absorption in water at two temperatures, fre- 
quency 12 Me/s. 


4. Discussion of results 


4.1. Velocity 

The velocity isobar for atmospheric pressure is 
in good agreement with published results, as indeed 
it should be in view of the accuracy expected. A 
convenient form of presentation of the pressure re- 
sults is to plot the isothermal curves in the form of 
v/v» against pressure. This is done in Fig. 9 where 
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all the isotherms of Fig. 6 lie within the shaded area. 
Curves representing. Brquarp’s, Hotron’s and 
Smiru and Lawson’s results are also given. It is ap- 
parent that the relatively large change in slope with 
temperature observed by Biquarp is not confirmed 
by the more recent results. From the present work 
it appears that around 30°C the slope of the v/v, 
curve shows a minimum and this conclusion is also 
supported by Smita and Lawson’s observations. This 
phenomenon is probably related to the minimum 
observed in the isothermal compressibility in the 
same temperature region. 
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Fig. 9. v/vp in water vs. pressure. 
Ba Biquarp, 
S. and L.: Smrrx and Lawson, 
He Hotton, 
eas Tair, 
Ba Bripeman, theor. (500 kg/cm?), 
Bs: Bripeman, theor. (1000 kg/cm?). 


The discrepancy between the results given in 
Fig. 9 of the various investigators amounts to a 
maximum of 4 or 5m/s at 10000 lb/sq’, a value 
in excess of that accounted for by the stated ex- 
perimental errors. The method of measurement 
used for the data presented in this work yielded ab- 
solute velocity values: this was not the case with the 


R.I. TAIT: ULTRASONIC PROPAGATION IN LIQUIDS 


_ACUSTICA 
Vol. 7 (1957) 


other measurements where the measured values are 


relative to a reference velocity (1510 m/s at 30° C_ 


and 1 atm). However the main source of the discre- 
pancy probably lies in the measurement of pressure. 
The lack of a uniform pressure standard is one of 
the most unsatisfactory aspects of precision high 
pressure work and many conflicting results in the 
high pressure field may be directly attributed to this 
cause. The gauge used for the present investigation 
was ‘a 12” sub-standard Budenberg Bourdon type 
(full scale reading 13 000 Ib/sq’’), which was cali- 
brated at the National Physical Laboratory. Both 
Hotron and Smita and Lawson used manganin re- 
sistance gauges which were designed to cover a much 
greater pressure range. It is therefore felt that the 
data obtained in this work are, within the range of 
pressure covered, more reliable. 


Perhaps the greatest value of sound velocity mea- 
surements lies in the fact that they provide a means 
of directly determining the adiabatic compressibility 
fy of a fluid. Then, provided the isothermal com- 
pressibility is known, it is possible to calculate the 
ratio of specific heats y, a constant which otherwise 
can only be derived indirectly with far less certainty 
from thermodynamic formulae. 


Table I gives the y values calculated from the 
various data available for pressures up to 1000 atm. 
Bripeman’s figures derived from purely thermo- 
dynamical considerations are given for comparison. 
It is also of interest to calculate the expected change 
with pressure in v/vg from Bripeman’s data alone; 
values so obtained are shown in Fig. 9. By represents 
the calculated value for a pressure of 500 kg/cm?. 
Better agreement is obtained at 1000 kg/cm? and 
B, represents the value extrapolated from this pres- 
sure to 500 kg/cm?, (assuming a linear law govern- 
ing v/vyvs. pressure over the range in question). 
In view of the errors involved in calculating v/v,» 
from thermodynamic data, the agreement between 
the experimental and theoretical values is quite satis- 
factory. It will be seen from Table I that although 
the figures for y obtained from the velocity measure- 
ments are undoubtedly to be preferred to the thermo- 
dynamie values, the discrepancy between the results 
of the different investigators indicates that a greater 
degree of experimental accuracy, particularly with 


Table I: Values of y; Water. 


Present results | HOLTON | SmirH and Lawson BRIDGMAN (theory) 

cape Temperature Temperature Temperature “patapeapaees 
30°C 50°C 30° C 50°C 30°C | 50°C 30° C 50°C 
0 1.020 1.040 1.019 1.039 1.020 1.039 1.01 1.04 
500 0.997 1.020 0.993 1.015 0.999 1.016 1.02 1.05 
1000 1.008 1.034 1.002 1.026 1.013 1.027 1.03 1.05 
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respect to pressure measurement, is required before 
the existing thermodynamic data can be amended 
with impunity. 

To return to the isobaric curves of Fig. 7, it ap- 
pears that the temperature maximum is displaced 
by pressure toward higher temperatures. This result 
confirms the findings of Smira and Lawson and as 
far as the pressure is taken in this work the results 
are in excellent agreement with these authors’. Hot- 
ton, in following the isotherms at 30° and 50° C, 
found the curves to converge and cross at about 
5.300 kg/cm? from which he concluded the direction 
‘of displacement of the velocity maximum was to- 
wards a lower temperature. This result has not been 
confirmed by either the present work or that of 
‘Smirx and Lawson and must therefore be taken to 
be in error. The reason for the discrepancy is dis- 
cussed at length by Smiru and. Lawson, who con- 
clude that it is due to inaccuracies in the measure- 
ment of pressure. 


4.2, Absorption 

The results of Lrrovirz and Carnevae for tem- 
peratures of 0° and 30°C are given in Table II, 
together with those obtained here at 16° and 23° C. 
The two sets of results fall nicely together giving 
a more general picture of the behaviour of «/f/? with 
pressure. A decrease of slope with temperature is 
apparent in the curves a/f? vs. pressure, the effect 


being more marked at the lower temperatures, but 
the relative change in a/f? appears to be roughly 
constant. 

In discussing the results, the logical approach is 
firstly to calculate the expected change of attenuation 
on “classical” grounds. Considering shear viscosity 
alone, a//? is given by the Stokes equation 


a 817, 8x? 9.V. 


P ai ov 3 3 ; 
where 7), is the coefficient of shear viscosity, V, the 
specific volume and v the sound velocity. The pres- 
sure coefficient of v has already been measured and 
the effect of pressure on 7, and V, is known from 
the work of Bripaman and others. 


(3) 


If however it is assumed that the observed excess 
absorption above the classical value is due to a bulk 
or compressional viscosity 7, then the appropriate 


a 2 x? 4 
(a) = oa (m ae 3 | > (4) 


expression is 


from which 
Eire eee wee 
My = (1 ey ) = 3 Ns (5) 
or 
4 Ao Ss. 
y= 3 Ns (Sa = 1] . (6) 


Values for CHER and el ea > Vobis | elaas and 


Table II 
| | ; | 

T p Om Bae | (i) - [o-17 | %obs. — Melass. | cobs. | ae | iim 
(°C) (10 Ai a) hice: a | ?? Jobs. . | PEA 1017 aclass. | Laing s | pa 

) 0 1731 | 57.4 40.3 3.36 1.792 5.64 

0 2 15.7 54.5 38.8 3.47 1.752 5.76 

0 4 14.5 51.6 ihe 3.56 1.716 5.85 

0 6 13.5 48.7 35.2 3.62 1.691 5.91 

Ox 8 12.7 | 45.9 | 33.2 3.62 1.674 5.85 

0 10 11.9 43.3 | 31.4 3.64 1.664 5.86 

OM eA 12 11.3 | 40.8 29.5 3.62 1.658 5.78 
16 0 9.3 30.8 | 21.5 3.31 ‘han 3.42 
16 2 8.8 29.6 20.8 3.36 1.104 3.47 
16 4 8.3 28.4 20.1 3.42 1.098 3.54 
16 6 7.9 27.1 | 19.2 3.44 1.091 3.55 
16 8 7.4 25.9 18.5 3.50 1.084. 3.61 
16 10 6.94 24.6 17.66 3.54 1.077 3.64 
164) 12 6.55 23.4 16.85 3.57 1.070 3.67 
23. ~C | 0 8.32 26.0 17.68 3.13 0.938 2.66 
23 2 7.93 24.9 16.97 3.14 0.938 2.68 
23 4 7.52 23.8 16.28 3.16 0.938 2.70 
23 6 7AS 22.8 15.62 3.17 0.938 2.71 
23 8 6.77 21.6 14.83 3.18 0.938 2.72 
23 10 6.44 20.6 14.16 3.19 0.938 2.73 
rN 12 6.12 19.5 13.38 3.19 0.938 2.73 
30 0 6.15 18.5 12.35 3.01 0.800 2.15 
30—s 2 5.9 17.6 7 2.98 0.806 2.14 
30 4 5.6 16.7 11d | 2.98 0.812 2.13 
30 6 5.4 15.8 10.4 2.93 0.817 2.10 
30 8 5.15 15.0 . 9.85 2.91 0.823 2.10 
30 10 4.9 | 14.2 9.3 2.90 0.829 2.10 
30 12 4.7 ) 13.4 B7 2.85 0.835 2.06 
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yy are given in Table II. It will be observed at 
all temperatures that both the compression and shear 
contributions to the losses decrease with pressure. ' 


The application of the Stokes equation to different 
liquids shows that (a/f?) class, usually exhibits a 
minimum at some particular pressure. This pheno- 
menon can in general be expected for all liquids as 
at a sufficiently high pressure the change in viscosity 
predominates over the changes in sound velocity 
(compressibility) and density. One might therefore 
expect the observed absorption to show a minimum 
value. The minimum will naturally occur at dif- 
ferent pressures for different liquids: in some, 
(water for example) the freezing point is reached 
before the minimum is attained. 


Within the pressure range of interest here the 
shear viscosity of water behaves abnormally in that 
at the two lower temperatures 77, exhibits a minimum 
as the pressure is raised. BripGMAN gives a satis- 
factory interpretation of this phenomenon in terms 
of an association effect. As can be seen from the 
table the compressional viscosity derived from the 
sound absorption data also exhibits anomalies, but 
whereas at 0° C a minimum is observed in 7,, 
appears to show a maximum. At 16° C no maximum 
was observed within the pressure range available 
but there is some evidence of one being reached at 
higher pressure. The shear viscosity at 30° C rises 
with pressure: Lrrovirz and CaRNevALE report a 
slight minimum in 7), at this temperature. 


Although it is undoubtedly connected with the 
association effect in water, the rather striking anti- 
thetic behaviour of the shear and bulk viscosities 
cannot be satisfactorily explained since the concept 
of compressional viscosity is a hypothetical one, and 
the mechanism giving rise to it is not fully under- 
stood. The most successful theory is that due to 
Hatt [7] who considers the excess absorption in 
water as being due to a structural relaxation mech- 
anism. Hatt assumes that the water molecule can 
exist in two states: in state 1 the volume occupied by 
a molecule is larger and the energy is lower than in 
state 2. Compression will therefore cause an increase 
in the relative number of state 2 molecules. Based on 
the above argument Hatt derives an expression for 
the excess absorption and by making some assump- 
tions, theory and experiment are shown to be in 
agreement. The effect of pressure on Hatt’s equa- 
tions is discussed in detail by Lirevirz and Carne- 
vaLE who conclude that the theory requires modi- 
fication before the observed pressure results can be 
explained. It becomes necessary to consider the open 
packing as being the higher free energy state. Rea- 
sonable agreement between experiment and theory is 
then obtained for both sets of pressure results. 
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However it must be allowed that some assumptions 
are made in Lirovitz’ treatment of the pressure ef- 
fect, as well as in the basic theory of Hatt, which 
are not altogether justifiable. The modification to — 
Hatw’s theory introduced by Lrrovirz is incompatible 
with most modern theories of the structure of water 
which require the closer packed molecule to be in the 
higher energy state. Recent investigations using X- 
ray diffraction indicate that the two-state model is 
over-simplified and that the molecular arrangements — 
in water are more complicated and that the mole- 
cular arrangements in water are more complicated. 
It is not surprising therefore that theories based on 
the two-state model do not give a complete inter- 
pretation of the physical properties of water. 


5. Conclusion 


An appreciable amount of data on the effect of 
pressure on the velocity of sound in pure water is 
now available through the work of several investiga- 
tors. This information is potentially of value in the 
field of thermodynamics but from an examination 
of the results obtained it is evident that with the 
present techniques, the accuracy which may be ob- 
tained is limited by the pressure measurement. This 
aspect of the work is considered worthy of further 
development. 

Absorption measurements have not hitherto been 
extensively undertaken for liquids under pressure. 
Great accuracy in pressure measurement is not es- 
sential for this work as the limit is set by the ultra- 
sonic technique employed. The results which may 
be obtained from such an investigation yield infor- 
mation on the effect of pressure on compressional 
losses in liquids and are valuable in that they provide 
a new line of approach for the study of anomalous 
absorption. 
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SCHALLDRUCKVERLAUF UND SPEKTRUM 
BEI DER SCHWINGUNGSKAVITATION 


von L. Boun 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Die Spitze eines magnetostriktiven Sondenmikrophons grofer Bandbreite (10 bis 1000 
kHz) wird in den Brennpunkt einer Schallstrahleranordnung von der Frequenz 15 kHz 
eingeftihrt. Es werden der Druckverlauf und das Spektrum der an dieser Spitze einsetzen- 
den Kavitationsvorginge gemessen. 

Bei schwacher Kavitation bilden sich Luftblasen, die durch die Kavitationssté%e in ihren 
Kigenresonanzen angeregt werden. Sie stabilisieren sich bei bestimmten Radien, die vom 
Luftgehalt des Wassers abhangen und bei untersattigtem Wasser unterhalb der Resonanz- 
groKe fiir die Anregungsfrequenz liegen. Bei starkerer Anregung zerplatzen sie. 

Bei starker Kavyitation bilden sich scharfe Druckimpulse aus, die das Spektrum weit- 
gehend bestimmen. 

Weiter werden Angaben iiber die Phasenlage der Impulse in bezug auf die Anregung, 
iiber Impulsform und -dauer und iiber zeitlich statistische Schwankungen um die periodi- 
sche Lage gemacht. 


Summary 


The pressures and sound spectra for cavitation set up in the focus of a transducer array 
of 15 ke/s are measured with a magneto-strictive probe. 

In weak cavitation, bubbles are formed and set in oscillation at their natural frequen- 
cies. They stabilise themselves with a definite radius (depending on air content), which 
can lie below that for resonance with the source, and burst under strong excitation. In hard 
cavitation strong pressure pulses are set up, which mainly influence the noise spectrum. 

Other matters considered are relationship between pulse and excitation; pulse form 
and duration; statistical fluctuations from the mean position. 


Sommaire 


On a introduit, au foyer d’un dispositif émetteur acoustique de fréquence 15 kHz, la 
pointe d’un microphone sonde a magnétostriction, et on a relevé la variation de la pres- 
sion et le spectre des phénoménes de cavitation ayant lieu a cette pointe. 

Dans le cas d’une faible cavitation, il se forme des bulles d’air excitées a leurs résonan- 
ces propres par les chocs de cavitation. Leurs rayons se stabilisent a des valeurs détermi- 
nées qui dépendent de la concentration de l’eau en air et, dans le cas d’une eau non 
saturée, sont inférieures aux valeurs de résonance relatives 4 la fréquence d’excitation. 
Quand l’excitation est plus forte, ces bulles crévent. 

Dans le cas d’une forté cavitation, il apparait des impulsions de pression aigués qui 
déterminent l’allure du spectre dans une large mesure. 

On donne en outre des informations sur la phase des impulsions par rapport a l’ex- 
citation, sur la forme et la durée de ces impulsions, et sur les variations statistiques dans 


le temps autour de la phase périodique. 


Einleitung 


In dieser Arbeit wird der Versuch unternommen, 
mit Hilfe einer akustischen Mefmethode Aufschlu8 
iiber den Ablauf des einzelnen Kavitationsprozesses 
bei periodischer Anregung durch ein starkes Schall- 
feld zu erhalten. Der Grundgedanke der Unter- 
suchungsmethode besteht darin, daB in einem durch 
eine fokussierende Schallstrahleranordnung erzeug- 
ten Schall-,,Brennpunkt* eine Grenzflache angebracht 
wird, an der bevorzugt oder sogar ausschlieflich 
Kavitation einsetzt. Diese Grenzflache wird durch 
die empfindliche Spitze eines Sondenmikrophons 
groBer Bandbreite gebildet, das den Druckverlauf 
beim Kavitationsproze8 unmittelbar registriert. Dab 


: 


Bio 5. 


man auf diese Weise nur Kavitation an Grenzflachen 
beobachten kann, diirfte wohl keine entscheidende 
Einschrankung darstellen, da nach verschiedenen frii- 
heren Untersuchungen angenommen werden kann, 
daf8 der charakteristische Ablauf des individuellen 
Kavitationsvorgangs an einer Grenzflache sich nicht 
grundsatzlich von dem im freien Medium unter- 


scheidet. 


1. MeBanordnung und Apparatur 


1.1. MeBanordnung 


In Abb. 1 ist die verwendete MeBanordnung skiz- 
ziert. In einem mit Glaswanden versehenen Wasser- 


a ee 
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Abb. 1. Schema der MefSanordnung. 1 Fokussierende 
Anordnung von Magnetostriktionssendern, 2 
Mikrophonspitze, 3 Nickeldraht, 4 Permanent- 
magnet, 5 Abnahmespule, 6 Absorptionsmate- 
rial, 7 Anordnung zur Gaserzeugung durch 
Elektrolyse. 


gefas von 60cm x 60cm x 80cm befindet sich als 
Schallsender 
21 parallelgeschalteten Magnetostriktionsschwingern 
(Nickelblechpakete von 45 mm x 45 mm x 105 mm), 
die auf 14,6 kHz abgestimmt sind. Die Leistungs- 
aufnahme betragt auf der elektrischen Seite maximal 


eine fokussierende Anordnung von 


1 kW. Ein magnetostriktives Sondenmikrophon (s. 
Abschnitt 1.2.) wird senkrecht zum Schallstrahl in 
den Fokus eingefiihrt. Infolge der induktiven Ab- 
nahme der Mikrophonspannung wird primar die 
zeitliche Druckanderung angezeigt, so daf} man erst 
noch eine Entzerrung durch einen Integrator vor- 
zunehmen hat, um den Druckverlauf direkt zu mes- 
sen. Die Mikrophonspannung kann entweder einem 
Analysator (selektiver Pegelmesser, Bandbreite 80 Hz, 
Frequenzbereich 3 bis 620 kHz) mit nachfolgendem 
Pegelschreiber oder den vertikalen Ablenkplatten 
eines Oszillographen zugefiihrt werden. Die Anzeige 
auf dem Oszillographen wird mit periodischer Ab- 
lenkung oder bei einmaliger Strahlablenkung beob- 
achtet. Da im allgemeinen der Kavitationsverlauf 
innerhalb einer Periode der Anregungsfrequenz in- 
teressiert, ist im Falle der einmaligen Ablenkung 
wegen der sehr kurzen Grundperiode von 68 us 
eine photographische Registrierung notwendig. Dazu 
wurde eine Kamera mit Vorsatztubus dem 
Oszillographen angebracht. Unter Verwendung einer 
Nachbeschleunigung am Oszillographen konnte der 
Druckverlauf tiber ein bis zwei Perioden der Grund- 
schwingung auf hochempfindlichen Réntgenfilm 
(,,Fluorapid“) abgebildet werden. Ein elektroni- 
sches Schaltgeraét (eine aus der Sinusspannung des 


vor 
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Generators hergestellte Impulsfolge wird durch einen 
Druckknopf-Kontakt auf einen bistabilen Multivibra- 
tor geschaltet, der einen einmaligen Rechteckimpuls 
regelbarer Dauer erzeugt) dient zur Auslésung der 
einmaligen Strahlablenkung und zur Hellsteuerung 
des Strahles wahrend der Ablenkung. Zugleich wird 
der beobachtete Vorgang dadurch vom Generator 
in der Phase synchronisiert. 

Um das Auftreten von Kavitationsblasen am Mi- 
krophon kontrollieren zu kénnen, wurde eine opti- 
sche Projektion der Mikrophonspitze in 10- bis 20- 
facher VergroBerang.senkrecht zur Zeichenebene in 
Abb. 1 vorgesehen. Zur Beleuchtung dienten ent- 
weder eine Bogenlampe oder eine Funkenstrecke. 

Das zum Fiillen des MeBgefaiBes benutzte Wasser 
konnte bis auf etwa 60% Luftgehalt untersattigt 
werden, indem es durch eine Diise in einen von 
einer Olpumpe_ evakuierten Glasballon gesaugt 
wurde. Nach langerem Stehen im Mefgefah nimmt 
der Luftgehalt dann langsam wieder zu, bis nach 
einigen Tagen Sattigung erreicht wird. Es wurde 
ausschlieBlich Leitungswasser verwendet. Besondere 
Vorkehrungen zur Reinigung durch Filtern oder 
ahnliches wurden nicht getroffen. 


1.2. Das Sondenmikrophon 


Wegen der beim Kavitationsproze} zu erwartenden 
sehr kurzen Impulse beim Zusammenfallen der Hohl- 
rdume mute ein Mikrophon mit médglichst groBer 
Bandbreite benutzt werden. Zu diesem Zweck wurde 
ein nach dem Magnetostriktionsprinzip arbeitendes 
Sondenmikrophon, das schon vor einigen Jahren in 
einer einfacheren Form ausprobiert worden war [1], 
weiterentwickelt und geeicht. Es besteht aus einem 
2,50 m langen Nickeldraht (Thermodraht DIN 43 712) 
von 1,0 bzw. 0,5 mm Durchmesser. Ein Ende befindet 
sich in der Fliissigkeit und ist durch einen Kunststoff- 
schlauch so abgeschirmt, daB nur die Endflache dem 
Schall ausgesetzt wird. Die Schallisolation wird durch 
ein zwischen Abschirmung und Draht befindliches Luft- 
polster erreicht, das durch einen spiralf6rmig hinein- 
gewickelten Wollfaden eine gleichmaBige Dicke erhalt. 
Am unteren Ende wird der Abschirmschlauch durch 
Lack oder Kitt luftdicht mit dem Nickeldraht verbun- 
den. Die durch den Schalldruck im Wasser im Nickel- 
draht erregten Longitudinalwellen sind infolge der 
Magnetostriktion mit der Magnetisierung gekoppelt, so 
da8 man an einer kleinen um den Draht gelegten In- 
duktionsspule eine induzierte Spannung erhalt, die pro- 
portional der zeitlichen Anderung der Magnetisierung 
ist. Da die Magnetisierung in Phase mit der Deforma- 
tion im Draht und diese wiederum in Phase mit dem 
an der Endflache wirkenden Druck ist, wird die Mikro- 
phonspannung proportional der zeitlichen Anderung 
des Druckes im Schallfeld. 

Um eine gleichmaBige Druck-Empfindlichkeit des 
Mikrophons in einem breiten Frequenzband zu errei- 
chen, sind drei Gesichtspunkte zu beachten: 

a) Reflexionen der mechanischen Welle am Draht- 
ende miissen vermieden werden, damit sich keine 
stehenden Wellen bilden kénnen. Dies wurde durch 
Anbringen einer in der Dicke keilférmig zunehmenden 
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Schicht von Absorptionsmaterial (,,Antidrdhnmittel“ 
oder Wachs-Kollophonium-Mischung) am freien Ende 
des Nickeldrahtes erreicht. Bei hohen Frequenzen wer- 
den die Wellen aufierdem durch die Absorption im 
Draht zunehmend gedampft. Abb. 2a zeigt eine Mef- 


ip oe 
Nickeldraht E 
8 
| 
[P 
Abnahme = p—-—--—— 1. 
| 
(a) MefS{anordnung, 
Anregung (b) und (c) Impulsbilder 


auf dem Oszillogra- 
phen ohne und mit 
reflexionsfreiem Ab- 


schluf. 


Abb. 2. Messung der Daémpfung von Longitudinalwel- 


reflexionsfreier 
Abschlu8 


len auf dem Nickeldraht und der Reflexion am 


Ende. 


anordnung zur Bestimmung der Absorption im Draht 
und des Reflexionsfaktors am Ende. Durch eine Er- 
regungsspule werden im Draht Schwingungsimpulse 
(Impulsdauer 0,2 ms; Folgefrequenz 250 Hz) ange- 
regt, die durch den Draht laufen und am Ende reflek- 
tiert werden. Abb. 2b zeigt die vom Mikrophon regi- 
strierten primadren und am Ende reflektierten Schwin- 
gungsimpulse bei 200 kHz, aus denen unter der An- 
nahme einer vollstaéndigen Reflexion am Drahtende die 
Dampfung im Nickeldraht abgeschitzt werden kann 
(Abb. 3a). Abb. 2c zeigt den Fall des durch Absorp- 
tionsmaterial reflexionsfrei gemachten Abschlusses. Bei 
200 kHz ist hier innerhalb der (nicht sehr hohen) 
MeBgenauigkeit keine Reflexion mehr wahrnehmbar. 


}=-—75cm 
brew. 150cm 


10° kHz 


Abb. 3. (a) Dampfung von Longitudinalwellen auf 
dem Nickeldraht, 
(b) Amplitudenreflexionsfaktor R an dem mit 
Absorptionsmaterial yersehenen  Draht- 
ende. 
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Jedoch wurde bei der tiefen Anregungsfrequenz yon 
14,6 kHz (Wellenlange auf dem Draht 35 cm!) nur ein 
Reflexionsfaktor von 20 bis 25% erreicht (Abb. 3b). 
Die hierdurch bedingten stehenden Wellen wurden 
in Kauf genommen, weil das Mikrophon nicht zu un- 
handlich werden sollte. 


b) Wegen des induktiven Effekts der Mikrophon- 
spule steigt die Druckempfindlichkeit mit der Frequenz 
an. Dies wurde mit der unter a) beschriebenen Mef- 
anordnung bestatigt (Abb. 4a, Kurve 1). Zur Entzer- 
rung wird ein integrierender Verstérker (kapazitiver 
Anodenwiderstand) hinter das Mikrophon geschaltet 


(Abb. 4a, Kurve 2). 
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Abb. 4a. 1: Frequenzkurve bei Ubertragung von 


Schwingungsimpulsen iiber den Nickel- 
draht, Anregung durch eine Hilfsspule, 
Ordinate willkiirlich. 

2: Frequenzkurve des Breitbandverstarkers 
mit Integrationsstufe, Ordinate uw: 107°. 

3: Akustische Eichung, Ordinate Volt/atm., 
A Mikrophonvergleich, © Schalldruck- 
waage. 


_-Nickeldraht, 0.5mm ¢ 


ee 
Papierrohre 


Abb. 4b. Konstruktion der Induktionsspule des magne- 
tostriktiven Mikrophons. 


c) Die Breite der Mikrophonspule mu8 klein gegen 
die kleinste zu messende Wellenlange der mechanischen 
Welle auf dem Draht sein, um einen Spalteffekt zu 
vermeiden. Bei 1 MHz betragt die Wellenlinge nur 
noch 5mm. Da die effektive Spulenbreite kleiner als 
die halbe Wellenlange bleiben muf, wurde die Spulen- 
wicklung etwa 1,2 mm breit ausgefiihrt (Abb. 4b). Der 
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Spulentrager besteht nur aus einer diinnen Zellophan- 
papierhiille, damit die Windungen (500 Windungen 
von 50u Durchmesser) den Draht méglichst eng um- 
schlieBen, mit seitlichen Stiitzen aus starkem Plexiglas, 
yon denen eine zur Befestigung der ganzen Spule an 
einer Aluminiumstange dient, welche zugleich zur Hal- 
terung des ganzen Nickeldrahtes benutzt wird. In 
Abb. 4a, Kurve 1, wird oberhalb 800 kHz der Abfall 
infolge des Spalteffekts deutlich. 


Nachdem durch die yorstehend beschriebenen Mab- 
nahmen eine gleichmafige Empfindlichkeit des Mikro- 
phons mit einer oberen Grenzfrequenz von nahezu 
1 MHz bei magnetischer Anregung sichergestellt war, 
wurde noch eine grobe akustische Absoluteichung vor- 
genommen (Abb. 4a, Kurve 3), und zwar bei den vier 
Frequenzen 15 kHz, 175 kHz, 365 kHz, 500 kHz durch 
Vergleich mit einem nach der Reziprozitatsmethode ge- 
eichten Turmalinmikrophon und im Bereich von 300 
kHz bis 1 MHz durch Vergleich mit einer Schalldruck- 
waage. Infolge der sehr unterschiedlichen Schallfelder 
im ersten Fall und der unterschiedlichen Struktur des 
Schallstrahles bei verschiedenen Frequenzen bei der 
zweiten Methode zeigen die Messungen starke Streu- 
ungen. Da es sehr fraglich erscheint, ob die akustische 
Eichung im freien Schallfeld fiir die Anzeige von star- 
ken Kavitationsst6Ben, die sich in kleinen Bereichen 
unmittelbar an der Endflache des Mikrophondrahtes 
selbst abspielen, iberhaupt mafgebend ist, wurde diese 
grobe Eichung als ausreichend betrachtet. 


Der Durchmesser des Mikrophondrahtes mu so 
klein gehalten werden, dai keine Schwingungsform 
mit Querunterteilung und damit eine Dispersion ein- 
treten kann. Daher wurde zunachst ein Durchmesser 
von 1mm gewéahlt. AuBerdem soll ein Sondenmikro- 
phon das Schallfeld méglichst wenig stéren. Spater 
ergab sich jedoch bei der Messung der Kavitations- 
impulse noch ein weiterer Gesichtspunkt, der die Wahl 
eines noch diinneren Drahtes ratsam erscheinen lief. 
Das Zusammenfallen der beim Kavitationsvorgang ge- 
bildeten Hohlraume findet im allgemeinen in sehr eng 
begrenzten Bereichen von vielleicht nur 10-3? oder 
10-4cm Durchmesser statt, so dah erst eine Art 
Transformation auf den viel groferen Drahtquerschnitt 
und damit eine Verminderung der durch den Impuls 
erzeugten Deformation eintritt, wahrend das anregende 
Schallfeld auf den gesamten Querschnitt des Draht- 
endes wirken kann. Infolgedessen waren die Impulse 
so schwach, da man sie auf der Sinuskurvye der An- 
regungsfrequenz kaum wahrnehmen konnte. Als dar- 
aufhin ein Draht von 0,5 mm Durchmesser verwendet 
wurde, traten die Impulse viel stérker heraus. Mit die- 
sem Drahtdurchmesser erfolgten samtliche im folgen- 
den wiedergegebenen Messungen. 

Die beim magnetostriktiven Wandler erforderliche 
magnetische Polarisation kénnte durch das remanente 
Feld geliefert werden, wie es beim Senderaggregat der 
Fall ist. Da sich der diinne Nickeldraht jedoch hinsicht- 
lich der Remanenz als sehr inhomogen erwies und man 
somit beim Verschieben der Induktionsspule lings des 
Drahtes unterschiedliche Polarisationen erhalt, was eine 
Eichung praktisch unméglich machen wiirde, wurde 
unmittelbar unter der Induktionspule parallel zum 
Draht ein stabformiger Permanentmagnet (8 mm Durch- 
messer, 25 mm lang) befestigt (Abb. 4b), der so stark 
ist, da in der Spule eine méglichst hohe Spannung 
induziert wird, da aber bei den héchsten zu erwar- 
tenden Amplituden noch keine Verzerrungen infolge 
magnetischer Sattigungserscheinungen auftreten k6n- 
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nen. Eine Beriihrung zwischen Stabmagnet und Nickel- 
draht wurde vermieden, tm keine Reflexionsstelle fiir 
die mechanische Welle zu schaffen. 


2. Theorie 
2.1. Kavitationstheorie 


Um die spatere Diskussion der MeBergebnisse zu 
erleichtern, sollen hier kurz einige Resultate von 
theoretischen Uberlegungen zum KavitationsprozeB 
angefiihrt werden, und zwar hauptsachlich unter 
Benutzung der Arbeiten von Notrinck und Neppiras 
[2], [3], in denen die Ansitze von Lord Rayieidén 
erweitert und speziellsauf den Fall der Schwingungs- 
kavitation in gashaltigem Medium angewandt wer- 
den. Die Autoren setzen die Differentialgleichung 
einer gashaltigen Keimblase in einem zeitlich ver- 
anderlichen Schallfeld p=py sin @¢ an und werten 
diese mit Hilfe einer elektronischen Rechenmaschine 
aus. Als Ergebnis erhalt man Beziehungen zwischen 


Schalldruck und BlasengroBe in Form der Abb. 5a. 


Schalldruck 
———} 


4 6 atm 8 
Schalldruck ——=- 


Abb. 5. (a) Zeitliche Anderung des Radius einer Gas- 
blase in einem starken sinusformigen 

Schallfeld (nach Noxrinck und Nepprras), 

(b) Zeit von der Entstehung bis zum Zusam- 
menfallen einer Kavitationsblase in Ab- 

hangigkeit vom Druck im Schallfeld (nach 


Notrinck und Nepprras). 


Die Blase wachst vom Anfangsradius Ry ab wahrend 
der Unterdruckperiode, wird nach Einsetzen des 
Uberdrucks gebremst, erreicht nach der Zeit ¢,, ihre 
maximale GroBe R,, und fallt dann in der Zeit t zu- 
sammen, allerdings, abweichend vom _klassischen 
Rayxeicuschen Ansatz, unter abnehmendem auferen 
Druck. AnschlieBend sind Nachschwingungen der 


gasgefillten Keimblase in ihrer Eigenfrequenz zu 
erwarten, 


4 
if 
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Von den zahlreichen Ergebnissen der Rechnungen 
von Noxrinck und Neppiras seien hier einige an- 
gefuhrt: 

Der maximale Blasenradius R,, ist fast unabhan- 
gig von der GréBe der Keimblase Ry. Zum Beispiel 
andert sich R,, bei einer Variation von Ry zwischen 
0,8 u und 80 uw, d.h. im Verhaltnis 1 : 100, nur um 

-etwa 20%. R,, wird wesentlich durch die Amplitude 

der anregenden Schallschwingung bestimmt. Da nun 
die Zeit t,, +t vom Anwachsen der Blase vom Ra- 
dius Ry ab bis zum Zusammenfallen wesentlich vom 
erreichten Radius R,, abhangt, wird also die Phase 
der beim Zusammenfallen entstehenden Druck- 
impulse in bezug auf die anregende Schallschwin- 
gung ebenfalls hauptsachlich von der Schalldruck- 
amplitude bestimmt werden (Abb. 5b). Selbstver- 
standlich hangt R,, auch von der Dauer der Unter- 
druckperiode im Schallfeld, doh. von der Schall- 
frequenz ab. Diese Abhangigkeit interessiert aber 
hier nicht, da in der vorliegenden Arbeit nur eine 
feste Anregungsfrequenz benutzt wurde. Nach Abb. 5 
ist zu erwarten, dai die Schallimpulse beim Zusam- 
menfallen der Hohlraume bzw. die dadurch ausge- 
Josten Blaseneigenschwingungen erst auftreten, wenn 
im Schallfeld das Druckmaximum erreicht oder tiber- 
schritten ist. 

Eine gewisse Einschrankung der Ergebnisse von 
Nottrnck und Neppiras bedeutet die Annahme einer 
ruhenden Keimblase zu Beginn des Kavitationspro- 
zesses. Wahrscheinlich wird die Keimblase noch in 
Nachschwingungen begriffen sein, wenn im Schall- 
feld der Unterdruck einsetzt. Dies kann zur Folge 
haben, daB der Kavitationsablauf in aufeinander- 
folgenden Perioden voneinander abweicht oder dai 
die maximalen BlasengréBen R,, und die Druck- 
impulse kleiner ausfallen. Trotzdem darf man die 
Ergebnisse wohl als eine gute Naherung ansehen. 


2.2. Zum Kavitationsspektrum 


In einer friiheren Arbeit wurden schon Kavita- 
tionsspektren gemessen [4], die allerdings nur qua- 
litativ zu bewerten sind, da die Frequenzkurve des 
benutzten Mikrophons nicht genau genug bekannt 
war. Dabei trat neben einem Linienspektrum aus 
Harmonischen und auch Subharmonischen der An- 
regungsfrequenz ein kontinuierliches Spektrum in 
Erscheinung, das auf statistische Schwankungen der 
Phase oder Amplitude der Vorgange hindeutet. Die- 
ser Rauschpegel ist oft als Kriterium fiir den Einsatz 
der Kavitation benutzt worden [4], [5], wobei man 
teilweise nur beim Auftreten eines kontinuierlichen 
Spektrums von ,,echter Kavitation** spricht. Dazu 
wire zu bemerken, dai} nach den Ergebnissen dieser 


1 Unter ,echter Kavitation“ wird in [4] die durch 
nur mit Wasserdampf gefiillte Blasen hervorgerufene 
Kavitation verstanden. 
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Arbeit der Kavitationsvorgang in einem stationaren 
Schallfeld quasi-periodisch ablauft, wenn man ihn 
an einer Stelle lokalisiert erzeugt und miBt. Es ware 
im Prinzip vorstellbar, daB unter extrem sauberen 
Bedingungen die Kavitationsvorgange so periodisch 
und in gleicher Starke auftreten, daB ein Rauschen 
neben dem Stérpegel der MeBapparatur nicht mehr 
wahrnehmbar ware, d.h. es ware Kavitation vor- 
handen, ohne daf ein kontinuierliches Spektrum ge- 
messen werden kann. Es erscheint daher zweckma- 
Biger, als Kavitations-Kriterium das Auftreten von 
Spektrallinien bei im Vergleich zur Anregungsfre- 
quenz sehr hohen Frequenzen zu benutzen, besser 
noch das Auftreten von Linien oder Rauschen, weil 
die Linienenergie im Falle starker  statistischer 
Schwankungen bei hohen Frequenzen zugunsten des 
Rauschpegels stark reduziert sein kann (siehe unten). 
Auf jeden Fall ist bei periodischer Anregung eines 
nichtlinearen Vorganges das Linienspektrum das Pri- 
mare; erst durch statistische Schwankung kann ein 
Teil der Linienenergie in Rauschenergie umgewan- 
delt werden. 

Bei den in dieser Arbeit berichteten Untersuchun- 
gen wurde die Kavitation im allgemeinen auf einen 
sehr kleinen Bereich an der Spitze des Sondenmikro- 
phons beschrankt. Der in den gemessenen Spektren 
auftretende Rauschuntergrund kann nur durch ge- 
ringe statistische Schwankungen der Amplituden 
und Phasen der Zeitvorginge um die Mittellage her- 
vorgerufen werden. Ein solcher Fall ist auf nach- 
richtentechnischem Gebiet von Macrariane [6] in 
einer theoretischen Arbeit behandelt worden, deren 
Ergebnisse hier kurz zitiert seien: Bei der stati- 
stischen Amplitudenschwankung einer 
periodisch wiederkehrenden Zeitfunktion ergibt sich 
auBer einem Linienspektrum, dessen Enveloppe das 
Spektrum des Einzelvorgangs G(w) ist, ein Rausch- 
spektrum mit der Energiedichte (Abb. 6a) 


oy? | G(w) P 


Ep(w) = rT 


(1) 
(6, mittlere quadratische Schwankung der Ampli- 
tude, 7 Periode der Wiederholung des Vorgangs). 

Bei statistischer Schwankung der 
Phase um die periodische Lage erhalt man bei 
Annahme einer Gaufiverteilung fiir die Energiever- 
teilung im Spektrum den Ausdruck 

r 


E(a) = Vie> re + (2) 
4+.@ 797%. oa 6(@—2mnx) ete) [2 


worin 67 die mittlere quadratische Schwankung 
der Phase um den periodischen Wert nT und 
0(@—2manx/T) die Diracsche 6-Funktion bedeutet. 
Der zweite Summand in der Klammer stellt also 
ein Linienspektrum dar, wahrend der erste ein kon- 
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tinuierliches Spektrum ergibt. Dieses Rauschspek- 
trum hat bei » =0 den Wert Null und steigt dann 
bei kleinen @ mit ?, bis es schlieBlich, sobald das 
erste Glied in der Klammer Eins ist, proportional 
zu G(w)? wird. Zugleich nimmt das Linienspektrum 
mit wachsender Frequenz ab (Abb. 6b), bis bei 
hohen Frequenzen fast nur noch statistische Energie 
vorhanden ist. 


\ l6lw)l? 


mittl. Energie der Linien—> 
—a 
+ 


wal 


clGlwhi? 


Energiedichte—> 


® 


Bei der Kavitation werden natirlich Amplitude 
und Phase schwanken, jedoch im allgemeinen nicht 
unabhangig voneinander, da ja eine besonders 
groBe Kavitationsblase zur Entstehung und zum 
Zusammenfallen eine langere Zeit benotigen wird 
als eine kleinere und zugleich einen besonders star- 
ken Impuls erzeugen wird, so das mit einer starken 
Phasenabweichung auch eine starke Amplituden- 
schwankung gekoppelt sein muB. Um den EinflufB 
einer solchen Kopplung abschatzen zu konnen, wurde 
unter Benutzung der Methode von Macrartane das 
Spektrum fiir den Fall einer linearen Kopplung in 
der Form 


An,=A+Ctn, O,=COp 


berechnet (A, Amplitude des n-ten Impulses und 
o, mittlere quadratische Amplitudenschwankung, 
A mittlere Amplitude, ¢, Abweichung und or mitt- 
lere quadratische Schwankung des Impulses von der 
periodischen Lage aus). 

Aus diesem Ansatz ergibt sich die spektrale Ver- 
teilung der Energie in der Form (vgl. Gl. (1) und 
(2)) (3) 


F (os) 02 COVE 5 j1- (1 4.042 072 2) e230" 
Prat 
(1+ 042 op? 2) e- F*2% §( 2.0m) eet 


Das bedeutet im wesentlichen eine noch starkere 
Beschneidung des Rauschspektrums bei tiefen Fre- 
quenzen und eine entsprechend geringere Abnahme 
der Linienenergie in diesem Bereich sowie eine et- 
was langsamere Abnahme der Linien nach hohen 
Frequenzen zu. 
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AuSer der Amplituden- und Phasenschwankung 
tritt natiirlich eine Schwankung in der Impulsdauer 
ein. Die Berechnung des Spektrums fiir diesen Fall 
wird sehr kompliziert und das Ergebnis sehr un- 
iibersichtlich, da mit der Impulsbreite der die Form 
des Spektrums wesentlich bestimmende Parameter 
verandert wird. Daher soll auf eine Wiedergabe des 
Resultats der naherungsweise durchgefiihrten Rech- 
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Abb. 6. Verlauf des Energiespektrums 
eines willkiirlich Lemma 
nen Zeitvorgangs 


(a) bei zeitlich  statistischer 
Amplitudenschwankung, 
(b) bei zeitlich  statistischer 


Schwankung der Phase; 
oben Linienspektrum, un- 
ten Rauschspektrum. 


nung verzichtet werden, zumal der Einfluf sich ge- 
gentiber dem der anderen Schwankungen als gering 
erwiesen hat. Im wesentlichen ergibt sich ein Rausch- 
spektrum, das an der Stelle der Einbriiche im nor- 
malen Impulsspektrum Maxima hat und so die Ein- 
briiche ,,ausfillt*, was ja zu erwarten ist. 


3. MeBergebnisse 


3.1. Zum Kavitationseinsatz 


Bei einer Untersattigung des Wassers bis auf etwa 
60% des Sattigungs-Luftgehaltes trat in einem Wech- 
selfeld von 3 bis 4 Atmospharen, oft allerdings erst 
nach einigen Minuten, spontan Kavitation am Mikro- 
phon auf. Dieser Schwellenwert stimmt etwa mit 
den Angaben von Srrassperc [7] und Brake [8] 
iiberein. Zur Erleichterung des Kavitationseinsatzes 
wurde eine Vorrichtung zur Erzeugung kiinstlicher 
gasgefiillter Keimblasen durch Elektrolyse vorgese- 
hake (Abb. 1). Diese gestattet entweder die Einschal- 
tung eines schwachen Gleichstroms (etwa 0,5 mA) 
von beliebiger Dauer oder die Entladung eines auf 
einige kV aufgeladenen Kondensators von 2 bis 
5 uF iiber die Mikrophonspitze zu einer Gegen- 
elektrode. Dabei wird im zweiten Fall eine groBe 
Zahl von Gasblasen, deren Gréfe je nach Dauer und 
Starke des Entladungsstromes zwischen wenigen 
und etwa 100 liegen kann, an der Endflache des 
Mikrophons gebildet und groBtenteils etwa 1 mm in 
das Wasser hineingeschleudert. Man erhalt so eine 
Wolke von Keimblasen, deren GréBe dem Augen- 
schein nach relativ homogen ist. Diese Methode 
wurde jedoch nicht weiter ausgebaut. Eine weitere 
Moglichkeit zur Erleichterung des Kavitationsein- 
satzes wurde kurz untersucht: Eine sehr starke Kon- 
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densatorentladung, die zu Stromen von einigen Am- 
pere im Mikrophondraht fihrt, wird durch eine 
KurzschluBfunkenstrecke nach 10°4s abgebrochen, 
so da kein auslaufender Gleichstrom-,,Schwanz“ 
vorhanden ist, auf den die Bildung von Blasen bei 
der vorher erwahnten normalen Entladung zuriick- 
zufiihren ist. In der optischen Projektion konnten 
auch keine Blasen festgestellt werden; diese miiBten 
sonst Durchmesser unter etwa 1 uw haben. Da der 
Widerstand der Wasserstrecke im wesentlichen an 
dem kleinsten Querschnitt der Mikrophonoberflache 
konzentriert ist, wird dort beim Stromdurchgang 
eine Erwarmung zu erwarten sein, womit sich wahr- 
_scheinlich der spontane Kavitationseinsatz, der auf 
diese Weise erzielt wurde, erklaren laBt. 


3.2. Kavitation bei schwacher Anregung 


Unter schwacher Anregung® sollen Schalldruck- 
amplituden bis zu etwa 1,5 Atmospharen verstanden 
werden. Damit lassen sich bei einem Luftgehalt von 
= 70% des Sattigungsgehaltes Kavitationsvorginge 
in Form von ,,StoBschwingungen“ gasgefiillter Bla- 
sen anregen. Das soll heifen: Eine luftgefiillte Keim- 
blase wird zu nichtlinearen Schwingungen angeregt, 
wobei sie sich wahrend der Unterdruckphase der 
Anregungsschwingung nur wenig vergrofert, so da 
sich beim anschlieBenden Zusammenfallen kein star- 
ker und scharfer Impuls ausbildet, sondern nur die 
Kigenschwingung der Keimblase ,,angestoBen“ wird, 
die dann als ,,StoBoszillation“ gedampft ausschwingt. 
Diese schwachen Kavitationssto{e enthalten keine 
sehr hohen Frequenzen, da sich nur Eigenschwin- 
gungen bis etwa zum 6- bis 8-fachen der Anregungs- 
frequenz beobachten lieBen. Kleinere Keimblasen 
lassen sich nur durch die bei starker Kavitation, 
d.h. bei hoheren Anregungsschalldrucken, erzeug- 
ten scharferen Impulse Zu Eigenschwingungen an- 
regen (s. Abschnitt 3.3.). 

Bei der Analyse dieser Stofschwingungen ergab 
sich eine sehr interessante Erscheinung: Die schwin- 
genden Luftblasen stabilisieren sich vorwiegend in 
bestimmten Gréfen und schwingen bei den entspre- 
chenden Pulsationseigenfrequenzen, die sich mit dem 
Luftgehalt des Wassers und zum Teil auch mit der 
Anregungsamplitude nur in diskreten Springen an- 
dern, so, als ware nur bei bestimmten charakteristi- 
schen Eigenfrequenzen eine stabile Blasenschwin- 
gung méglich. Dabei beobachtete man zugleich, wie 
sich von der schwingenden Blase stoBweise Wolken 
von winzigen Blaschen ablésen. Jn nahezu luftgesat- 
tigtem Wasser bilden sich auch oft an bestimmten 
Stellen im freien Medium Kavitationsprozesse aus 
(von Brake ,,trapped bubbles“ genannt). Die Blasen 
schieBen dabei ruckartig hin und her, als wirden 
sie von einer Rakete getrieben, und lassen jedesmal 
einen Schweif winziger Blaschen hinter sich (Abb. 
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17a zeigt in der Mitte links einen solchen Effekt). 
Ihre GroBe bleibt dabei konstant. Es wurden auch, 
wie schon von Brake erwahnt, Blasen beobachtet, 
die mit einer Frequenz von einigen Hertz auf einem 
Kreis vom Radius 1 bis 3 cm umlaufen, wobei die 
Kreisbahn, genau betrachtet, aus vielen kleinen ge- 
raden Stiicken zusammengesetzt ist, die von der 
Blase ruckartig unter Zuricklassung einer Wolke 
kleiner Blaschen durchlaufen wird. Eine solche Kreis- 
bewegung wird oft eine halbe Stunde véllig gleich- 
maf ig aufrechterhalten, ohne daf sich die GroBe der 
Blase, (die, wohlbemerkt, meist kleiner als die Re- 
sonanzgrofe ist) andert. 


Als Erklarungsméglichkeit fiir diese beobachtete 
Stabilisierung wurde eine als ,,Faraday crispation“ 
bezeichnete Erscheinung herangezogen, bei der durch 
Anregung sehr starker Kapillarwellen auf ebenen 
oder auch Kugelflachen entweder kleine Luftblaschen 
abgeschniirt und in die Fliissigkeit geschleudert oder 
Wassertropfchen in den Gasraum zerspriiht werden. 
Uber diese Erscheinung speziell bei Luftblasen in 
Wasser und Wassertropfchen ist eine kurze Vor- 
tragszusammenfassung von Wittarp [9] bekannt. 
Witiarp beobachtete besonders starke ,,crispations“ 
bei bestimmten Kombinationen von Blasenradius 
und Anregungsfrequenz. Da nun seit langerem be- 
kannt ist, daB eine stark nichtlinear schwingende 
Luftblase ihren Luftinhalt vergro8ert, ein Vorgang, 
der von Brake als ,,gleichgerichtete Diffusion“ be- 
zeichnet wird (und dessen Ursachen noch nicht ganz 
klar sind [10]), ware also eine Moglichkeit fiir 
eine VergroBerung der Blase gegeben. Andererseits 
kann durch starke Oberflachenschwingungen bei der 
»Faraday crispation™ der Luftinhalt durch Abschni- 
ren kleiner Blaschen herabgesetzt, also die Blase 
verkleinert werden, Eine Stabilisation bei bestimm- 
ten Blasenradien ware also méglich, sobald die An- 
regung einer Kapillarschwingung an einem bestimm- 
ten Radius gebunden ist. In diesem Falle konnte die 
Gaszufuhr durch ,,gleichgerichtete Diffusion“ kom- 
pensiert werden durch die Abtrennung von Blaschen 
bei der ,,Faraday crispation™. 

Witiarp auferte vor einiger Zeit die Vermutung, 
da® sich die fiir eine ,,Faraday crispation™ giinsti- 
gen Kombinationen von Blasenradius und Anre- 
gungsfrequenz wahrscheinlich durch das Zusammen- 
fallen der Grundpulsation nullter Ordnung mit ir- 
gendwelchen Oberschwingungen hoherer Ordnung 
ergaben. Dieser Gedanke fiihrt auch tatsachlich zu 
der Moglichkeit einer Kopplung der Oberschwin- 
gungen an bestimmte Blasenradien und damit zu 
einer Stabilisationsmoglichkeit. Fiihrt man namlich 
in die Differentialgleichung der Kapillarschwingun- 
gen eine auf die gesamte Oberflache wirkende sinus- 
formige Beschleunigung ein, so erhalt man eine 
Mathieusche Differentialgleichung. Damit ist eine 
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Méglichkeit gegeben, im Falle der Blasenschwingung 
aus der Grundpulsation Energie in die Oberflachen- 


wellen zu tbertragen und diese anzuregen. Da nun‘ 


aber nur bei bestimmten Radien gleichzeitig die 
Schwingung nullter Ordnung und eine Oberschwin- 
gung irgendeiner Ordnung existieren kénnen, gibt 
es ,,bevorzugte“ Radien bei der Anregung von Ka- 
pillarwellen. Bei schwacher Anregung tritt bevorzugt 
eine subharmonische Oberschwingung mit der hal- 
ben Frequenz der Schwingung nullter Ordnung auf. 
Bei starkerer Anregung durch die Grundpulsation 
bildet jedoch auch die Oberschwingung von gleicher 
Frequenz wie die anregende Schwingung nullter 
Ordnung eine Losung der Mathieuschen Differential- 
gleichung [11]. Man hat jetzt also die Blasenradien 
zu ermitteln, fiir welche die Schwingung nullter 
Ordnung mit irgendeiner der zahlreichen Schwin- 
gungen hédherer Ordnung zusammenfallen kann. Am 
Mikrophon wird dabei natirlich immer nur die 
Grundpulsation angezeigt, da die Schwingungen 


hoherer Ordnung praktisch nicht abstrahlen. 
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Abb. 7. Grundschwingung nullter Ordnung und Ober- 
schwingungen verschiedener Ordnung einer 


Gasblase in Wasser als Funktion des rezipro- 
ken Radius. 


In Abb. 7 sind die Schwingungen hoherer Ord- 
nung (f,) zusammen mit der Grundschwingung 
nullter Ordnung (fo) als Funktion des reziproken 
Radius dargestellt. Die Schnittpunkte von f, und f, 
ergeben die Bedingung fiir das gleichzeitige Auf- 
treten beider Wellentypen von gleicher Frequenz, 
die Schnittpunkte von f)/2 und f, fiir das Zusam- 
menfallen von fy mit einer Oberflachenschwingung 
der Frequenz f)/2. In der Tabelle I sind diese 
»Eigenwerte* zusammengestellt. 

Die Abb. 8 und 9 enthalten zwei verschiedene ge- 
messene ,,StoBschwingungen“, die jeweils bei ein- 
maliger Strahlablenkung auf dem Oszillographen 
photographiert wurden *, und zwar einmal mit vol- 


2 Samtliche in‘ der Arbeit wiedergegebenen Oszillo- 
gramme zeigen nach oben Uberdruck und nach unten 
Unterdruck an. 
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Tabelle I 
| _berechnete 
_ Eigenfrequenzen ; 
ovbile Eige | fiir gleich- Anregungs- 
poara le 18en- peitige Anregung bedingungen 
ses bac der nullten und | 
| m-ten Ordnung 
| fo =fn: |\fo =2fn| Luft- | Schall- 
| | gehalt | druck 
Abb. kHz | kHz n/| kHz n) % atm 
| 
1560 2 
16b 400...450 | 470 3] | GOs.) ROS Iiat eos eae 
i 390 2 | 
12c 200...220 | 210 “4 ea CAN: 3 
81ers LL 3270.22.80 
11215 | 0,5...1,0 
9 Gb i 0a sOnne es 6 | 80...90 1,0 
| 52,5 4| 
*) | 43,6 7 
29.7 8 | 
19b 25...30 28” ileus 2, Le 0,5 
10 22 21,5 9 | 
10 16 16,6 6) 
15,8 10 


*) In diesem Bereich wurden Schwingungen beobachtet, — 
deren Frequenz sich schwer bestimmen lieB. 


ler Bandbreite einschlieBlich der Anregungsfrequenz — 
von 14,6kHz, in den Abb. 8c’ und 9a’, 9b’ aufer- 
dem noch einmal tiber ein Hochpaffilter von 40 kHz — 
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Abb. 8. Blasenschwingungen bei schwacher Kayitation. 
Eigenfrequenz 110 bis 120kHz, Anregungs- 
frequenz 14,6kHz, Schalldruckamplitude (a) 
<~0,5atm, (b) ~0,7atm, (c) ~1,0atm, 
Luftgehalt des Wassers 70 bis 80%; Spektrum 
(csp) entspricht Oszillogramm (c); Spektrum 
(d) : Oberwellen yor Kavitationseinsatz, Schall- 
druck 2 bis 3 atm., Bandbreite bei (a), (b) 
und (c) 10 kHz bis 1 MHz, bei (c’) iiber Hoch- 
paB 40 kHz. 
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Abb. 9. Blasenschwingungen bei schwacher Kavitation. 
Figenfrequenz 65 bis 70kHz; Anregungsfre- 
quenz 14,6kHz, Schalldruckamplitude: (a) 
und (a’) 0,8 bis 1,0 atm, (b) und (b’) 1,0 bis 
1,5 atm, Luftgehalt des Wassers 80 bis 90%; 
das Spektrum (bsp) entspricht dem Oszillo- 
gramm b); Bandbreite bei (a) und (b) 10 kHz 
bis 1 MHz; bei (a’) und (b’) tiber Hochpah 
40 kHz. 


unterer Grenzfrequenz, wobei im wesentlichen die 
Anregungsfrequenz weggeschnitten wird. In den 
letztgenannten Oszillogrammen kann man die Ver- 
anderung der Blasengrofe wahrend der Kavitations- 
schwingung verfolgen; z. B. sieht man in den Abb. 
9a’ und 9b’, daB kurz vor dem Einsetzen einer 
neuen StoBschwingung die Periodendauer schnell 
abnimmt, d.h. also, da die Blase schnell kleiner 
wird, wahrend die Perioden davor relativ grof sind, 
was auf einen grofen Blasenradius hinweist. Die 
ersten Perioden bei sehr grofer Amplitude sind 
wahrscheinlich infolge der Nichtlinearitat vergro- 
Bert. — Der Nachschwingvorgang der Abb. 8 ent- 
halt etwa sieben bis acht Perioden auf eine Periode 
der Anregungsfrequenz von 14,6 kHz; das entspricht 
etwa 110 bis 120 kHz. Aus Abb. 9 ergibt sich etwa 
eine Frequenz von 65 bis 70 kHz. Bei denselben 
Frequenzen zeigt sich in den Spektren (Abb. 8 und 
9) ein Kontinuum, das auf statistische Schwankun- 
gen beim Ablauf der Nachschwingungen hindeutet. 
Die Form entspricht etwa der ,,Glockenkurve“ einer 
gedimpften Schwingung. Dieses Rauschspektrum 
gibt die Frequenzlage der Eigenschwingungen im 
allgemeinen deutlicher wieder als das Linienspek- 
trum, 

Vergleicht man die gemessenen bevorzugten Fre- 
quenzen mit den berechneten in Tabelle I, so ergibt 
sich fiir 65 bis 70 kHz eine Stabilisationsméglich- 
keit durch das Zusammenfallen der Schwingung null- 
ter Ordnung mit der von sechster Ordnung bei 


gleicher Frequenz. Fiir die Schwingung bei 110 bis 
120 kHz ergeben sich sogar zwei Moglichkeiten: 
Das Zusammenfallen der nullten und fiinften Ord- 
nung bei gleicher Frequenz oder der nullten Ord- 
nung mit der dritten Ordnung von der halben 
Grundfrequenz. Interessant ist, da sich bei der 
Frequenz von 110 bis 120 kHz die Nachschwingung 
schon bei sehr geringen Schalldrucken anregen lief 
(vgl. Abb. 8a), im Gegensatz zu der bei 65 bis 
70 kHz, was auf eine Stabilisierung durch die leich- 
ter anzuregende Subharmonische, mithin auf die 
Frequenz von 117 kHz hindeutet. An dieser Stelle 
mag auch erwahnt werden, daB im Abschnitt 3.3. 
bei groBeren Schalldrucken auch die Anregung von 
stabilen Schwingungen bei 200 bis 220 kHz und bei 
etwa 450 kHz beschrieben wird, die durch das Zu- 
sammenwirken der Schwingung nullter Ordnung mit 
denen der vierten bzw. siebenten Ordnung gleicher 
Frequenz erklart werden konnten. 
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Abb. 10. Spektrum von Blasenschwingungen in luft- 
gesattigtem Wasser. Schalldruck etwa 1 atm. 


In Abb. 10 ist ein Spektrum bei Kavitation in 
Leitungswasser bei schwacher. Anregung wiederge- 
geben, wobei direkt neben der Linie der Anregungs- 
frequenz von 14,6 kHz noch eine Spitze im kon- 
tinuierlichen Spektrum bei etwa 16 kHz auftritt; 
dort liegen auch etwa zwei mogliche ,,stabile“ Fre- 
quenzen (15,8 kHz und 16,6 kHz, siehe Tabelle I). 
Dasselbe gilt fiir die Spitze zwischen der 1. und 2. 
Spektrallinie (etwa 22 kHz). 


Abb. 11. Blase an der Mikro- 
phonspitze (Kavitation ent- 
sprechend Abb. 8c). Blasen- 
durchmesser etwa 65, Be- 


4 
leuchtung durch Bogenlampe 
mit 1/100 s. Imm 


In der Abb. 11 ist noch eine Abbildung der Mi- 
krophonspitze in etwa 20facher Vergr6Berung wie- 
dergegeben, die zugleich mit den Messungen der 
Abb. 8 gemacht wurde. Abb. 11 ist mit einer langen 
Belichtungszeit von etwa 1/100 s aufgenommen und 
zeigt eine stabil schwingende Blase, deren Durch- 
messer von 60 bis 65 .u etwas groBer ist als der 
einer 110- bis 120-kHz-Blase, was wahrscheinlich 
auf die sehr starken Schwingungsamplituden zuriick- 
zufiihren ist. 
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In der Tabelle I sind die Luftsattigungen und 
Anregungsamplituden angegeben, bei denen die ent- 
sprechenden Nachschwingungsfrequenzen bevorzugt 
auftreten. Diese Werte sind nur als grobe Anhalts- 
punkte zu werten, da die Kinstellung bestimmter 
Blasengrofen in sehr komplizierter Weise von Luft- 
gehalt und Schallamplitude abhangen. Bei der Er- 
héhung der Schallamplitude beobachtet man oft 
einen Ubergang von einer ,,stabilen“ Blasengrofe 
zur niachst gréBeren stabilen Blase, wahrscheinlich 
infolge verstarkter ,,gleichgerichteter Diffusion“, 
falls der Luftgehalt gro® genug ist. Bei zu geringem 
Luftgehalt scheint auch das Umgekehrte méglich zu 
sein. Vielleicht wird in diesem Fall durch sehr 
starke Oberflachenschwingungen soviel Luft abgege- 
ben, da die Blase unter Umstanden eine geringere 
Gro8e erreichen kann, um sich dort zu stabilisieren, 
wenn es tiberhaupt zu einer Stabilisierung kommt. 
Denn um stabile Eigenschwingungen zu erreichen, 
ist anscheinend ein ganz bestimmtes Verhaltnis von 
Luftgehalt und Anregungsamplitude notwendig. Dies 
ergab sich bei dem Versuch, mit Hilfe kiinstlicher 
Gaszufuhr durch Elektrolyse in untersattigtem Was- 
ser die verschiedenen moglichen Eigenresonanzen 
nacheinander herzustellen. Es gelang dabei, durch 
passende Abstimmung von Gaszutubr und Schall- 
feldamplitude auBer den schon erwaéhnten Schwin- 
gungen bei 110 und 67 kHz weitere stabile Schwin- 
gungen bei 40 bis 50 kHz und bei etwa 30 kHz 
einzustellen. 

Die Dampfung der Blaseneigenschwingungen hielt 
sich in der tblichen Grofienordnung der in verschie- 
denen Arbeiten, z. B. [12], an schwingenden Blasen 
gemessenen Werte, kann jedoch zum Teil bis zu 50% 
hoher liegen, was sich vielleicht durch Energieabgabe 
an die Oberschwingungen erklaren lieBe. 


3.3. Starke Kavitation 


3.3.1. Untersadttigtes Wasser 
70 bis 80% Luftgehalt 


von etwa 


Bei starker Schallanregung, d.h. bei Amplituden 
von etwa 1,5 Atmospharen bis maximal 3 bis 4 
Atmospharen, werden die Kavitationsvorgange sehr 
kompliziert und kénnen in verschiedenen Formen 
ablaufen. Im folgenden wird unter Abschnitt a) auf 
einige interessante ,,Zwischenphasen* eingegangen, 
die nicht immer auftreten (siehe die Abb. 12, 13 
und 14), wahrend unter b) der normalerweise ein- 


tretende Ablauf behandelt wird (Abb. 15 ff.). 


a) Eine Erhéhung der Schalldruckamplitude iiber 
0,5 bis 1,0 atm fiihrt zuerst oft zu einer Anregung hé- 
herer Resonanzen zugleich mit den im Abschnitt 3.2. 
beschriebenen ,,StoBschwingungen“ (Abb. 12a), die 
auch im Spektrum hervortreten kénnen. Zum Beispiel 
ist in Abb. 13b auf er einer Anhebung des gesamten 
Pegels bei hohen Frequenzen deutlich ein ,,Formant- 
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Abb. 12. Starke Kavitation. eer bei 
(a) 1,5 bis 2,0 atm, bei (b), (c), (d) etwa 
3 atm, Luftgehalt des Wassers 70 bis 80%; 
das Spektrum (cgp) entspricht dem Oszillo- 
gramm c, Eigenfrequenzen: (b), (c) 200 bis 
220 kHz, (d) 110 bis 120 kHz. 


bereich“ des Rauschpegels um 200 kHz wahrnehmbar. 
Zugleich beobachtete man in der optischen Projektion 
in unregelmafigen Abstaénden ein Zerplatzen der sta- 
bilen Blase, wobei jedesmal ein Strahl von kleinen 
Blaschen von der Mikrophonoberflache fortgeschleudert 
wird, und zwar zum Teil senkrecht zum Schallstrahl. 
Erst nach etwa 0,5 bis 1 cm Laufweg werden die Blis- 
chen vom Strahlungsdruck erfaBt und in Richtung des 
Schallstrahles mitgefiihrt, wo sie sich, je nach Luft- 
gehalt des Wassers, entweder bald auflésen oder sich 
an bestimmten Stellen des Schallfeldes zu einer gro- 
Beren Blase vereinigen. Nach dem Zerplatzen einer 
Blase bildet sich, meist an derselben Stelle, eine neue, 
die innerhalb yon 1/10 bis 1/100s bis auf die gleiche 
GroBe anwiachst, die die zerplatzte Blase hatte. Bei 
einer weiteren. Steigerung des Schalldrucks tiber etwa 
1,5 atm sieht man nur noch zerplatzende Blasen am 
Mikrophon, oft mehrere gleichzeitig, die unter Umstan- 
den bei zwei verschiedenen GréSenstufen zerplatzen 
kénnen. Das Anwachsen und Zerplatzen erfolgt jetzt 
so schnell, da8 man mit dem Auge nicht mehr folgen 
kann. Man sieht nur das Auftauchen und Verschwin- 
den von schwarzen Punkten gleicher Gréfe. Abb. 14a 
zeigt die Bahn der abgeschleuderten Bruchstiicke zer- 
platzender Blasen in einer Streulichtaufnahme bei 1 s 
Belichtungsdauer. Im oberen Teil ist die Mikrophon- 
spitze als Silhouette erkennbar. In Abb. 14b sind drei 
willkiirlich herausgegriffene Phasen mit Funkenbeleuch- 
tung aufgenommen worden. Die Blasen auf dem ersten 
Bild (Pfeile) haben etwa ,normale“ GréBe, die auf 
dem zweiten sind stark vergroBert und befinden sich 
anscheinend in starker Bewegung, und im dritten Bild 
sind sie offensichtlich zerplatzt. In den sehr unregel- 
maBigen Oszillogrammen treten meistens wenige, oft 
nur drei oder vier Schwingungen héherer Frequenz, 


z. B. 200 kHz, stark hervor (Abb. 12b). Die Verande- 
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Abb. 13. Kayitationsspektren bei steigendem Anre- 
gungsschalldruck. Druckamplitude bei (a) 
1,0atm, (b) 2,0atm, (c) und (d) 3 bis 
4atm; Luftgehalt des Wassers bei (a), (b), 
(c) 80%; bei (d) 60 bis 70%. 


Abb. 14. (a) Streulichtaufnahme der vom Mikrophon 
fortgeschleuderten Blischen; Belichtung 
ls, Druckamplitude etwa 2 atm. 

(b) Kavitationsblasen in drei willkiirlichen 
Phasen, Luftgehalt 80%, Druckamplitude 
2 bis 3 atm, Funkenbeleuchtung 2 bis 
B>107-*s. 
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rung des Spektrums beim Ubergang von schwacher 
Anregung zu starkerer Anregung zeigt Abb. 13. Wah- 
rend bei Abb. 13a nur eine starke Schwingung bei 
110 bis 120 kHz hervortritt, wird bei Abb. 13b der 
Pegel bei hohen Frequenzen angehoben, und es tritt 
zugleich ein zusatzliches Maximum bei 200 kHz auf. 
Bei Abb. 13c und d zeigt sich ein starkes Kontinuum 
bei hohen Frequenzen mit einem flachen Maximum 
bei 200 kHz (groBerer Luftgehalt) oder sogar bei 400 
bis 500 kHz (geringerer Luftgehalt). Das Linienspek- 
trum fallt nach hohen Frequenzen sehr stark ab, so 
da dort fast nur noch statistische Energie vorhanden 
ist, was auf besonders starke Schwankungen bei Kavi- 
tationsvorgaingen dieser Art schliefen laBt. (In fast 
allen Spektren sind bei starker Kavitation die Linien 
bis 50 kHz, zum Teil auch bis 100 kHz, etwas iiber- 
steuert, um bei der beschrankten Dynamik des Ana- 
lysators das Rauschspektrum noch mit erfassen zu 
konnen!) 


Gelegentlich kommt auch der Fall vor, da sich Bla- 
sen trotz starker Anregung (bis zu 2,5 oder 3 atm) 
stabilisieren und nicht zerplatzen. Sie fiihren dann 
Schwingungen mit extrem groBen Amplituden aus. 
Zwei Beispiele sind in den Abb. 12c¢ und d fiir Fre- 
quenzen von 200 bis 210kHz? (Spektrum cg, in 
Abb. 12) bzw. 110 bis 120 kHz dargestellt. Die Auf- 
nahmen enthalten auch die Anregungsfrequenz von 
14,6 kHz (Abb. 12c bei einmaliger Strahlablenkung, 
Abb. 12d Zeitaufnahme bei periodischer Ablenkung). 
Wenn man annimmt, da’ es sich bei den am Mikro- 
phon beobachteten, in den Umrissen etwas unscharfen, 
aber vollig runden Blasen um einzelne, bei der im 
Oszillogramm angezeigten Eigenfrequenz schwingende 
Blasen handelt, so erhalt man fiinf- bis zehnmal zu 
grofe Blasenradien; d. h. der Maximalradius Rp miiBte 
das fiinf- bis zehnfache des Ruheradius R, betragen. 
Theoretisch waren nach den Rechnungen von Nottinex 
und Neprrras derartige Verhiltnisse von R»/Ry mog- 
lich. Dann lieBe sich allerdings die beobachtete gleich- 
mafige Eigenfrequenz kaum verstehen. 


b) Sehr oft treten die unter a) beschriebenen 
Kavitationserscheinungen gar nicht auf. Dann beob- 
achtet man bei einer Erhéhung der Schallamplitude 
lediglich, wie die Blaseneigenschwingungen zuneh- 
mend spitzer werden, wie z.B. in den Abb. l5a 
und b, und wie sich schlieBlich bei noch starkerer 
Anregung ein einzelner starker Impuls ausbildet 
(Abb. 15 ¢ und c’). Die Spektren bei starker Kavi- 
tation sind nur sehr schwer zu analysieren, weil bei 
den in dieser Arbeit vorliegenden Verhaltnissen hin- 
sichtlich Luftgehalt des Wassers und Schalldruck- 
amplitude die entstehenden Kavitationsimpulse ge- 
rade eine solche Dauer haben, dafi ihr Spektrum 
im MeBbereich schon Nullstellen und Maxima zei- 
gen kann. Dieses tiberlagert sich also dem durch 
eventuell auftretende Blasenschwingungen erzeugten 
Spektrum. Bei starker Kavitation konnen, wie es 
scheint, mehrere Blasenschwingungen gleichzeitig 
auftreten, die zuweilen miteinander abwechseln. Da- 
zu kommt die Modifikation der Spektren durch sta- 
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Abb. 15. (a), (b), (c) Ubergang von Blasenschwingun- 
gen in einen starken Impuls bei Steigerung 
des Anregungsdrucks; (csp) Beispiel eines 


Spektrums beim Auftreten starker Impulse. 
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Abb. 16. 


Kinige spezielle Spektren bei starker Kavita- 
tion ((a) ‘und (c)), (b) Oszillogramm des 
zeitlichen Differentialquotienten des Druckes; 
Frequenz der angeregten Eigenschwingungen 
etwa 450 kHz. 


tistische Schwankungen der verschiedenen Parameter. 
Aus diesen Griinden laBt sich ein solches Spektrum 
nur schwer genau erklaren, selbst unter Zuhilfe- 
nahme des Oszillogramms, aus dem sich beim gleich- 
zeitigen Auftreten mehrerer Eigenschwingungen von 
Blasen deren Frequenzen meist auch nicht genau 
bestimmen lassen. Daher sollen im folgenden nur 
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einige Beispiele aus der groBen Zahl der gemesse- 
nen Spektren wiedergegeben werden. 

Das Spektrum in Abb. 15 zeigt im Gegensatz zu 
denen der Abb. 13, die unter a) beschrieben wur- 
den, bei etwa 200 kHz einen tiefen Einbruch. Aus 
dem Impulsoszillogramm Abb. 15 c’ kann man ent- 
nehmen, daf} die Impulsdauer t~7 us betragt. Wie — 
im Abschnitt 4.2. spater gezeigt wird, kann man die 
Impulsform durch einen halben Kosinus annahern, 
d.h. aber, da die Nullstellen des Spektrums bei 
den Frequenzen 1,5/t, 2,5/t, 3,5/t usw. liegen, in 
unserem Falle die erste bei etwa 210 kHz, die nach- 
ste bei 350 kHz. Bei diesen Frequenzen liegen aber 
gerade Minima im gemessenen Spektrum, so dab 
man dessen Verlauf wohl zum grofen Teil durch die 
Impulsform erklaéren muf. 

Andererseits treten aber in den Oszillogrammen, 
z. B. in den Abb. 15 und c’, neben den Impulsen 
deutlich Nachschwingungen auf, die sich auch im 
Spektrum sehr oft, z. B. durch ungewohnlich hohe 
Nebenmaxima oder andere vom Impulsspektrum 
abweichende Formen bemerkbar machen (Abb. 21 
Cgp, die bei stark untersattigtem Wasser ohne merk- 
liche Blasenschwingungen aufgenommen wurde, 
kommt der ,,Normalform“ der Impulsspektren sehr 
nahe). In den Abb. 16a und c sind zwei Beispiele 
angefihrt, bei denen der Rauschpegel bei hoheren 
Frequenzen ganz besonders hoch liegt. Der Abb. 16a 
entspricht das in Abb. 16b wiedergegebene Oszillo- 
gramm (hier ist der zeitliche Differentialquotient 
des Druckes aufgenommen worden, um die héheren 
Frequenzen anzuheben, damit sie sich besser aus 
dem Verstarkerrauschen herausheben), in dem deut- 
lich eine Schwingung von etwa 450 kHz erkennbar 
ist, die auch im Spektrum Abb. 16a besonders her- 
vortritt. In diesem Frequenzgebiet besteht auch die 
Moglichkeit einer Stabilisierung (siehe Tabelle I). 

Am Mikrophon beobachtet man bei sehr starker 
Kavitation wieder das schon erwahnte Zerplatzen 
von Blasen, wobei jetzt die ,,Bruchstiicke“ mit sehr 
groBer Geschwindigkeit geradlinig auf bestimmten 
Bahnen fortgeschleudert werden. Abb. 17a zeigt 
eine Streulichtaufnahme dieser Blasenbahnen, die 


(a) (b) 

Abb. 17. Streulichtaufnahme der vom Mikrophon fort- 
geschleuderten Blaschen bei sehr starker Ka- 
vitation. Luftgehalt des Wassers bei 
(a) 80 bis 90%, (b) 60%. 
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2 bis 4 mm lang sein kénnen, und interessanterweise 
ihre Richtung, die senkrecht zum Schallstrahl liegen 
kann, beibehalten. Wahrscheinlich haben sich an be- 
stimmten Stellen der Mikrophonflache durch die 
Kavitationseinwirkungen Locher gebildet, welche 
die fortgeschleuderten kleinen Blaschen in bestimmte 
Richtungen lenken. Abb. 17 a wurde in starker luft- 
haltigem Wasser von mindestens 80% Luftgehalt 
aufgenommen. In starker untersattigtem Wasser von 
nur 60 bis 70% Luftgehalt fliegen die Blaschen nicht 
so weit (Abb. 17b), wahrscheinlich weil sie sich 
schnell auflosen oder weil sich bei der ,,harteren“ 
Kavitation kleinere Blaschen bilden. 

c) Bei mindestens 80% Luftgehalt oder mehr und 
maximaler Schalldruckamplitude von etwa 4 atm 
bilden sich zuweilen sehr viel groBere Kavitations- 
blasen aus, die dann beim Zusammenfallen beson- 
ders kurze Impulse erzeugen (Abb. 18), deren Dauer 
oft nur 3 us betragt. Derselbe Effekt tritt ein, wenn 
man bei starker Kavitation in weniger lufthaltigem 
Wasser. durch Elektrolyse zusatzlich Gas erzeugt. 
Abb. 18 ag, zeigt ein dabei aufgenommenes Spek- 
trum, bei dem erst bei etwa 450 kHz ein Einbruch 
zu erkennen ist. Die UnregelmaBbigkeit kommt durch 
die kiinstliche Gaszufuhr. 
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Abb. 18. Zwei Beispiele besonders kurzzeitiger Kavita- 
tionsimpulse. Das Spektrum entspricht etwa 
den Oszillogrammen. Die Oszil- 
logramme wurden iiber einen 40- 
kHz-Hochpa8 aufgenommen. 


Abb. 19. Druckverlauf bei Kavitation in 
gesattigtem Wasser. 
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3.3.2. Kavitation in gesattigtem oder 
tibersattigtem Wasser 


Die Kavitationsvorgange in gesattigtem oder nahe- 
zu gesattigtem Wasser laufen ahnlich ab, wie im 
Abschnitt 3.3.1.,c) beschrieben, solange sich keine 
Blasen bilden, die bei der Anregungsfrequenz von 
14,6 kHz in Resonanz kommen, d.h. die Radien 
von etwa 0,22 mm erreichen, was bei sehr starker 
Anregung selbst bei etwas untersattigtem Wasser 
moglich ist. AuBer am Mikrophon treten an vielen 
Stellen im freien Medium ,,Kavitationsherde“ auf. 
Da die Abstrahlung der groSen Blasen schon recht 
betrachtlich ist, kann das Mikrophon auch durch 
diese in einiger Entfernung stattfindenden Prozesse 
angeregt werden, sofern sich nicht am Mikrophon 
unmittelbar eine groBe Blase bildet, was praktisch 
fast immer der Fall ist, die dann den Einflu8 des 
Schallfeldes weitgehend abschirmen kann. 


Wie schon erwahnt, bilden sich noch starke Ka- 
vitationsstoBe aus, solange die Blasen unterhalb der 
Resonanzgrofe fiir 14,6 kHz bleiben. Abb. 19 a zeigt 
ein Oszillogramm eines scharfen Impulses, dessen 
Spektrum sich von dem in Abb. 18 nur durch starke 
Schwingungsvorgange bei sehr tiefen Frequenzen 
unterscheidet. Auch der Einbruch bei etwa 450 kHz 
liegt an der gleichen Stelle. In Abb. 20 ist eine 
Streulichtaufnahme der 
gegeben, auf der sich ein ,,Kreislauf“ ausgebildet 
hat. Etwas schrag nach links werden in einem Strahl 


Mikrophonspitze wieder- 


kleine Blaschen weggeschleudert, die dann vom 
Strahlungsdruck nach rechts mitgenommen werden, 
sich dabei langsam vergrofern und aufsteigen und 
dann schlieBlich wohl durch eine sich ausbildende 
Strémung, zum groBten Teil wieder zum Kavitations- 
zentrum an der Mikrophonspitze zuriickgetrieben 
werden. 

Sobald die Blasen die Resonanzgrofe fiir 14,6 kHz 
erreichen, treten kaum mehr scharfe Impulse auf 
(Abb. 19 cc). Die sich an der Mikrophonspitze bil- 


denden grofien Resonanzblasen scheinen das Schall- 
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Abb. 20. Streulichtaufnahme bei starker Kavitation in 
nahezu luftgesattigtem Wasser mit einem Bla- 
sen-,,Kreislauf“, Schalldruck etwa 3 atm. 


feld weitgehend abzuschirmen, denn man beobachtet 
auf dem Oszillographen Schwingungsformen von 
der Art der in den Abb. 19b* und c gezeigten, die 
fiir stark nichtlinear schwingende Blasen charakte- 
ristisch sind. Das zugehorige Spektrum ist in Abb. 
19 cg, wiedergegeben. Grofe Blasen lassen sich an- 
scheinend besonders leicht zu subharmonischen 
Schwingungen in bezug auf die Anregungsfrequenz 
anregen. Abb. 19d zeigt ein Spektrogramm bei 
schwacherer Anregung in gesattigtem Wasser. Es 
treten deutlich Spektrallinien bei 1/4 und 1/2 und 
bei Vielfachen der Anregungsfrequenz auf. 


Uber die Ursache fiir das Auftreten yon Subharmo- 
nischen sind schon sehr viele Uberlegungen angestellt 
worden. Darauf soll hier nicht naher eingegangen wer- 
den, jedoch sollen einige Beobachtungen berichtet wer- 
den, die vielleicht auf eine weitere Deutungsméglichkeit 
hinweisen: 

Subharmonische werden nicht nur in der Umgebung 
der Grundfrequenz, sondern zuweilen bis zur 30. bis 
40. Harmonischen hinauf beobachtet, meist als ,,halb- 
zahlige* Harmonische, die manchmal nur 10 dB unter 
den ganzzahligen liegen. Diese hauptsachlich auftreten- 
den halbzahligen Vielfachen der Grundfrequenz weisen 
oft eine von den normalen Harmonischen abweichende 
»Enveloppe“ auf, deren Maximum bzw. Minimum dort 
liegt, wo bei starker Kavitation ein kontinuierliches 
Spektrum lige. Zum Beispiel tritt an die Stelle des 
scharfen Kontinuums bei 110 kHz im ‘Spektrum der 
Abb. 8 beim Herabsetzen des Schallpegels ein Linien- 
spektrum yon ,,halbzahligen“ Harmonischen. Sehr oft 
liegen schon auf dem Kontinuum solche Zwischenlinien 
(Abb. 13a und b). Auf Grund dieser Beobachtungen 
k6nnte man yvermuten, dafi sich bei der Schwingungs- 
kavitation irgendwelche der resultierenden Effekte, Im- 
pulse oder Blasenschwingungen nicht in jeder Periode 
gleich ausbilden, sondern z. B. in aufeinanderfolgenden 
Perioden abwechselnd starker und schwacher auftreten. 
Das bedeutet praktisch eine periodische Amplituden- 
und zugleich wahrscheinlich Phasenmodulation mit der 
halben normalen Folgefrequenz und ergibt Linien zwi- 


4 Abb. 19b entspricht sehr gut dem in einer kiirz- 
lich fertiggestellten Arbeit [131] mit optischen Metho- 
den gemessenen Verlauf des Blasenradius iiber eine 
Kavitationsperiode. 
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schen den normalen Harmonischen. Eine Méglichkeit 
fiir das alternierende Verhalten lieBe sich vielleicht fin- 
den, wenn man in den theoretischen Uberlegungen ein- 
mal die bei Kavitation in gashaltigem Wasser meist 
vorliegende ,,Randbedingung™ beriicksichtigen wiirde, 
da8 zu Beginn der Unterdruckperiode im Schallfeld 
die Eigenschwingungen der beim vorhergehenden Zu- 
sammenfallen angeregten Blasen noch nicht abgeklun- 
gen sind, was gewif den Ablauf des neuen Kavitations- 
prozesses beeinflussen k6nnte. 


3.3.3. Kavitation in stark untersattig- 
tem Wasser (Luftgehalt < 60%) 


In Wasser-von etwa 60% Luftgehalt bilden sich 
bei der hochstmoglichen Schalldruckamplitude von 
3 bis 4 atm scharfe Impulse aus, die im allgemeinen 
keine oder nur sehr schwache Blasenschwingungen 
zur Folge haben (Abb. 21a, b, c). Das Spektrum 
Abb. 21 cg, scheint ausschlieBlich durch die Impuls- 
form bestimmt zu sein und laBt sich am besten mit 
dem Spektrum eines Kosinusimpulses von etwa 
6 bis 7 us Dauer in Einklang bringen. 
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Abb. 21. Druckverlauf und Spektrum bei Kavitation in 
stark untersattigtem Wasser. Luftgehalt 60%; 
Schalldruck bei (a), (b) 3 atm, bei (c) 3 bis 
4 atm, das Spektrum (csp) entspricht etwa 


(c). 
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4. Phasenlage, Schwankung und Form der 
Kavitationsimpulse 


4.1, Phase der Impulse in bezug auf die Anregungs- 
schwingung 


Fiir die Bestimmung der Phase der Impulse in 
bezug auf die Phase der anregenden Schwingung 
von 14,6 kHz ergab sich eine Schwierigkeit daraus, 
daB fir 14,6 kHz auf dem Nickeldraht des Sonden- 
mikrophons stehende Wellen auftreten, weil der 
Reflexionsfaktor am Ende nur 20 bis 25% betragt. 
Man kann sich leicht tiberlegen, daB dadurch die 
Lage der Impulse auf der Sinuskurve von 14,6 kHz 
bis zu +10°, bezogen auf 360° fiir die Grund- 
periode, verschoben werden kann. Dieser Fehler 
wurde dadurch herabgesetzt, da8 die Abnahmespule 
des Mikrophons méglichst genau im Abstand einer 
Wellenlange vom Drahtende befestigt wurde, was 
sich durch einen Phasenvergleich des mit einem 
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Kristallmikrophon unmittelbar neben der Spitze des 
Sondenmikrophons gemessenen Schallfelddrucks mit 
dem yom Sondenmikrophon tiber den Integrations- 
verstarker angezeigten Druck leicht einstellen lief. 
Aufferdem wurden die Kavitationseffekte mit einem 
Kristallmikrophon aufgenommen, wobei sich inner- 
halb von etwa 5% der Grundperiode Ubereinstim- 
mung in der Lage der Impulse ergab. Diese Ge- 
nauigkeit dirfte ausreichend sein, zumal wegen der 
meist etwas verzerrten Kurvenform der Sinuskurve 
die Ausmessung der Impulslage noch ungenauer 
werden kann. 

Aus Abb. 21b entnimmt man, bezogen auf den 
Beginn der Unterdruckperiode, eine Zeit ¢,,+7 fiir 
das Anwachsen und Zusammenfallen des Kavitations- 
hohlraumes von etwa 57 us. (Die letzten 20% der 
Oszillogramme sind etwas komprimiert; dieser Teil 
darf nicht ausgewertet werden;) Aus Abb. 6 ent- 
nimmt man einen von Norrinck und Neprrras [3] 
berechneten Wert von 60 us bei 3 bis 4 atm. Diese 
Ubereinstimmung kann als befriedigend angesehen 
werden. 


Um auf bequeme Art die Anhangigkeit der Impuls- 


lage von der Schallamplitude zu bestimmen, wurde 
der ,,Nachteil“, daB die Senderamplitude durch den 
1-kW-Kraftverstarker mit 100-Hz-,,Brumm“ modu- 
liert ist, ausgenutzt, indem die Brummkompensation 
ganz ausgeschaltet und auf diese Weise eine 100-Hz- 
Modulation der Schalldruckamplituden im Verhalt- 
nis 14/10 erzielt wurde. Unter diesen Verhaltnissen 
ergaben sich die Oszillogramme in den Abb. 22a, 
b, c. Bei periodischer Ablenkung mit 20 bis 30 kHz 
wurde der Verlauf des Schalldrucks (Abb. 22a) und 
des zeitlichen Differentialquotienten des Druckes 
(Abb. 22b) bei starker Kavitation in nahezu ge- 
sittigtem Wasser photographiert, wahrend Abb. 22 c 
den Druckyerlauf bei Synchronisation mit der Netz- 
frequenz zeigt. Die Hohe der Impulse schwankt sehr 
stark mit 100 Hz (nach Notrinck und Nepriras mit 
der vierten Potenz des Schalldrucks, das ware hier 
im Verhaltnis (1,4)4~4). In den Abb. 22a und b 
sieht man deutlich die Verschiebung der Impulse 
zwischen zwei Positionen, die etwa 7 bis 8 us aus- 
einanderliegen. Entnimmt man wieder aus der Abb. 6 
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on 


die Anderung von ¢,,+t mit der Druckamplitude, 
z. B. fir 4 und 2,8 atm, so erhalt man eine Impuls- 
verschiebung von 6,5 bis 7 us. Wenn man beriick- 
sichtigt, dafS die Kurve in Abb. 6 fiir einen sehr 
kleinen Keimradius (3,2 1) berechnet ist und die 
Kurve fiir gashaltigeres Wasser mit groBeren Blasen 
trotz der geringen Abhangigkeit vom Keimradius 
etwas steiler verlaufen wurde, so erhalt man einen 
etwas groferen Wert, der dem experimentell be- 
stimmten noch naher kommt. 


4.2. Impulsdauer und Impulsform 


Aus den Oszillogrammen ist zu entnehmen, dah 
die Impulsdauer unter den vorliegenden Verhaltnis- 
sen von etwa 7 bis zu 2 bis 3 us betragen kann. Die 
Impulse sind besonders kurz und scharf, wenn der 
Luftgehalt gering oder wenn der maximale Blasen- 
radius besonders gro ist. 

Die Impulsform scheint zwischen einem oben ab- 
gerundeten Dreieck und einem halben Kosinus zu 
liegen (Abb. 22d und e) und Jaft sich vielleicht am 
besten durch die Kosinusform annahern. Die vor- 


Abb. 22. 
(a), (c) Druck und (b) zeitlicher Differential- 
quotient des Druckes bei Kavitation in luft- 
gesattigtem Wasser bei 100-Hz-Modulation 
des Senders. (a), (b) Kippfrequenz 10 kHz, 
(c) 50 kHz. (d), (e) Druckverlauf iiber Hoch- 
paB 40 kHz bei 25kHz Kippfrequenz; (f), 
(g) zeitlicher Differentialquotient des Druk- 
kes, (f) bei 25 kHz, (g) bei 1,5 kHz Kipp- 


frequenz. 


dere Flanke ist oft, besonders bei starker Kavitation, 
etwas steiler, was sich auch aus dem zeitlichen Dif- 
ferentialquotienten (Abb. 22 f, g) ergibt, der im po- 
sitiven und negativen Impuls die Steilheiten der 
vorderen und hinteren Flanke anzeigt. Schwiachere 
(Abb. 22 f) 


trisch zu sein. 


Impulse scheinen annahernd symme- 
Uber die wirklich im Kavitationsimpuls auftre- 
tenden Drucke lat sich aus den Messungen nichts 
entnehmen, da, wie im ersten Teil der Arbeit schon 
einmal erwahnt wurde, eine Drucktransformation 
auf den vollen Querschnitt des Nickeldrahtes statt- 
findet, bei der die Deformationsamplitude stark her- 
abgesetzt wird. Wenn man angeben konnte, auf 
welchem Querschnitt die DruckstéBe beim Zusam- 
menfallen eines Hohlraumes an der Mikrophonspitze 
wirken, koénnte man mit Hilfe der gemessenen Im- 
pulse, die bis zu 8 atm betragen, die wirklichen 
Kavitationsdrucke abschatzen. Das Transformations- 
verhaltnis ware z. B. bei Annahme eines Einwirkungs- 
querschnitts von 50 u etwa 100, d.h. man kame auf 
wirkliche Kavitationsdrucke von nahezu 1000 atm. 
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4.3. Mittlere quadratische Schwankung der Kavita- 
tionsvorgdnge um ihre periodische Lage 


Mit Hilfe der in Abschnitt 2.2. angegebenen Formel 
(3) fiir die Energieverteilung im Spektrum bei sta- 
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malem Wasser bei anderen als in den in Tabelle I an- 
gefiihrten Frequenzen Blasen stabilisiert werden kén- 
nen. Daher soll hier auf eine ausfiihrliche Wiedergabe 
der Experimente bei veranderter Oberflachenspannung 
verzichtet werden. 


Abb. 23. Spektren bei schwacher Kavi- 
tation in Fliissigkeiten mit 
herabgesetzter Oberflachen- 
spannung. Druckamplitude 
bei (a) etwa 1 atm, (b) etwa 


1,5 atm, (c) 1 bis 2 atm (ver- 
gleiche (c) mit Abb. 13). 


tistischer Schwankung yon Amplitude und Phase last 
sich eine Abschaétzung der mittleren quadratischen 
Schwankung or um die periodische Lage durchfihren, 
wenn man bei einer bestimmten Frequenz das Verhalt- 
nis von mittlerer Linienenergie zu Rauschenergie aus 
den gemessenen Spektren entnimmt. Diese Abschatzung 
ist von der speziellen Form des Zeitvorgangs unabhan- 
gig. (Die direkte Messung ist wegen der schon er- 
wahnten schwachen Amplitudenmodulation des Schall- 
feldes mit 100 Hz nur schwer moglich.) Die Rechnung 
soll hier nicht im einzelnen wiedergegeben werden; es 
sei nur bemerkt, da sich ein relativer Fehler bei der 
Bestimmung des Energieverhaltnisses nur zu 1/50 als 
relativer Fehler in o7 auswirkt. 

Bei starker Kavitation ergeben sich maximal Werte 
von 2 bis 3:10°®s fir die mittlere quadratische 
Schwankung der Kavitationsvorginge, und zwar fiir 
die Abweichungen nach einer Seite. Die gesamte 
Schwankungsbreite der Impulse ergibt sich dann als 
207 ~4 bis 6:10 ®s = 4 bis 6us, d.h. etwa 6 bis 
9% der Anregungsperiode. Dieser Wert kénnte gréfen- 
ordnungsmaig bei den beobachteten Kavitationsvor- 
gangen moglich sein. Daraus folgt weiter, dafi der Ein- 
satz der Kavitationsimpulse relativ periodisch erfolgt. 


4.4, EinfluB der Oberflichenspannung auf die 
Blasenschwingungen 


Zum Schluf8 soll noch erwahnt werden, da zur Prii- 
fung der Hypothese der Blasenstabilisation durch Ober- 
flachenschwingungen auch einige Messungen in Fliissig- 
keiten mit herabgesetzter Oberflachenspannung (Sa- 
poninlésung 0,4 g/l; Nitrobenzol) erfolgten. Bei nied- 
rigeren Oberflachenspannungen liegen die Oberflichen- 
resonanzen bei tieferen Frequenzen (proportional der 
Wurzel aus der Oberflaichenspannung), und die sich 
aus dem Zusammenfallen mit der Schwingung nullter 
Ordnung ergebenden stabilen Frequenzen erhéhen sich 
proportional der Oberflichenspannung. Es wurden tat- 
sachlich solche Verschiebungen in der berechneten Gré- 
Benordnung beobachtet; z. B. traten Blasenschwingun- 
gen bei etwa 90 und 140 bis 150 kHz auf (Abb. 23), 
als die Oberflachenspannung auf 57 dyn/cm _herabge- 
setzt wurde. Bei diesen Frequenzen, besonders bei 
150 kHz, waren sonst keine stabilen Blaseneigenschwin- 
gungen beobachtet worden. Trotzdem miissen diese Ex- 
perimente mit Vorbehalt aufgenommen werden, solange 
der ,,Stabilisationsmechanismus* nicht véllig geklirt ist 
und damit sichergestellt ist, da& nicht auch in nor- 


@)” 50 100 200 kHz ©) 10 50100 200 kHz © 10 50100 200 400 600 kHz 


Herrn Professor Dr. E. Meyer mochte ich an die- 
ser Stelle fiir die Anregung zu dieser Untersuchung 
und fur zahlreiche forderliche Diskussionen meinen 


herzlichen Dank zum Ausdruck bringen. 
(Eingegangen am 20. Juli 1956.) 
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ageow bJECTIVE LOUDNESS SCALE 
by D. W. Rosryson 


Communication from the National Physical Laboratory 


Summary 


Discrepancies between published determinations of the sound intensity/loudness relation 
are not satisfactorily accounted for by differences in the test procedures. Statistical exami- 
nation suggests that the variations, both between different observers and between the 
various inyestigations, may be interpreted as random errors in the subjective estimation 
of sensation ratios. The function relating loudness to loudness level is derived from the 
data. For practical purposes a simple approximation to it is proposed (2-fold change 
equivalent to a 10-phon step). Further work is necessary to establish the transformation 
function for use in computing the loudness of complex sounds from their spectra. It is 
probably similar to but may be distinct from the loudness function. 


Sommaire 


Les écarts entre les relations intensité sonore/intensité subjective qui ont été publiées 
ne sexpliquent pas d’une maniére satisfaisante par la différence des modes opératoires 
employés dans les déterminations. Un examen statistique permet de penser que les variations 
constatées, soit d’un observateur 4 |’autre, soit d’une détermination 4 l’autre, peuvent étre 
dues 4 des erreurs aléatoires dans l’estimation subjective des rapports de sensation. On 
établit, 4 partir des résultats obtenus, la fonction reliant l’intensité auditive au niveau de 
force sonore. On en propose une approximation simple (une multiplication par 2 équivant 
a une différence de 10 phones) qui convient pour la pratique. D’autres études seront 
nécessaires pour établir la fonction transformation a utiliser pour calculer l’intensité sub- 
jective de sons complexes a partir de leur spectre. Cette fonction est probablement ana- 
logue, mais peut étre distincte de la fonction intensité subjective ordinaire. 


Zusammenfassung 


Die Unterschiede, die zwischen den bisher veréffentlichten Untersuchungen iiber die 
Beziehungen yon Intensitéit und Lautheit auftreten, konnen nicht befriedigend durch die 
unterschiedlichen Versuchsbedingungen erklirt werden. Eine statistische Uberlegung lat 
yermuten, dafi die Abweichungen auf der natiirlichen Streuung der Urteile tiber die Ver- 
haltnisse von Sinnesempfindungen beruhen. Die Funktion zwischen Lautheit und Lautstarke 
wird aus den gegebenen Daten ermittelt. Fiir praktische Anwendungen wird eine einfache 
Naherung gegeben (Verdopplung entspricht einem 10-phon-Schritt). Erst auf Grund 
weiterer Arbeiten wird es moglich sein, die Funktion aufzustellen, die es gestattet, die 
Lautheit von Schallereignissen aus ihren Spektren zu berechnen. Diese Funktion ist wahr- 


scheinlich ahnlich, aber doch verschieden von der Lautheitsfunktion. 


1. Introduction 


This paper is concerned primarily with the re- 
lation between changes in the loudness of sounds 
and the changes in the physical intensity producing 
them, as determined by methods of direct subjective 
judgment. From the practical standpoint, interest 
in this relation arises from a demand to express the 
loudness of sounds in the units of a scale taking 
account of the characteristics of hearing, and free 
from the arbitrariness of equivalent loudness scales, 
such as the phon scale. It also has important theo- 
retical consequences. It has been suggested [1] for 
example that indirect methods, for instance the 
comparison of monaural and binaural sensations, 
or the loudness of tone-pairs in relation to the 
loudness of the components, may lead to a scale of 
loudness essentially the same as that derived by 


direct judgments of ratio. The physiological hypo- 
thesis underlying this idea is the arithmetical ad- 
ditivity of component loudnesses, and on it rests a 
number of proposed methods for the calculation of 
loudness from objective spectrum analysis. The 
identification of the direct scale with one having 
the additivity property cannot, however, be assumed 
but must rest on experimental confirmation. A com- 
parison of the scales of loudness obtained by direct 
and indirect methods is made in a later section. 


Methods of direct estimation of relative loudness 
have been employed in various forms by at least 
12 authors, and it is with the results of these and 
their apparent conflict that this paper is mainly con- 
cerned. The tests fall into two principal classes, 
ratio methods and absolute loudness estimations, 
though each type of test may be carried out by a 
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variety of procedures. Two principal reasons for the 
divergent results have been put forward, namely, 


the influence of test procedure on the process of’ 


loudness judgment, and the failure of observers to 
correlate their sensations accurately with the scale 
of numbers, Garner [2] has described an experi- 
ment under extreme conditions in which different 
observers were persuaded to give results at the 
centres of various non-overlapping ranges of the 
stimuli, but it is clear that the constraints required 
to bring this result above were incompatible with 
the aim of deducing the loudness scale, and would 
have been rejected by investigators having this aim 
in view. Some new experiments have been carried 
out, in which substantial alterations in test proce- 
dure were deliberately introduced, corresponding to 
the sort of variations met with in the 12 investiga- 
tions mentioned, rather than to “restrictive prac- 
tices’ such as were the object of Garner’s study. 
The experiments failed to show any substantial in- 
fluence of the test procedure. In another context 
Garner [3] has suggested that the scaling of loud- 
ness by a given observer may be consistent and yet 
differ from another observer’s equally consistent 
results due to a failure to correlate sensation with 
number, and has shown that if the individual loud- 
ness scales are erected without the intervention of 
numerical judgment (this is done by the determina- 


tion of equal loudness-intervals and e qual loud- 


ness-ratios) the results for different observers appear 
more consistent with each other. Garner’s results, 
however, lack conviction, owing partly to the rather 
untypical results obtained by his group of observers, 
and also to some extent, to the sensitivity of his 
method to very small changes in the estimated bi- 
section level of loudness intervals (see Appendix 1). 
The disparity between Garner’s results and others’ 
is further aggravated when the deduced loudness 
scale units are reintroduced'. Interpretation of 
Garner’s results also presents some difficulty, for 
two sounds may appear to his observers to stand in 
the loudness relation 2:1, yet he is able to infer 
that the “true” relation is 1.7: 1. Can it be said 
that either of these figures is more valid than the 
other? 


2. The experimental data and its treatment 


With a few exceptions, the test procedures adopt- 
ed by different workers differed in material ways, 


1 The general consensus of results is that a reduc- 
tion of about 11 phons corresponds to a halving of the 
loudness sensation for sounds of moderate intensity. 
Garner, however, obtained a value of about 17 by di- 
rect ratio estimations, but a still higher value, about 
21, as the indirect result of the bisection and equal- 
ratio technique. 
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and in order to make a proper comparison, the 
various data have first to be presented in commen- 
surate units. 


Perhaps the most obvious difference between the 
investigations is the variety of sounds used as test 
material. In view of the divergence between the 
equivalent loudness relations for such sounds as 
low-pitched pure tones, white noise and speech, it 
cannot, be expected that loudness scales based on 
decibel intervals will be the same for all. Stevens 
[4] has proposed for practical purposes that the 
scale of loudness for sounds of 50 dB or above may 
be uniquely related to the sound pressure level pro- 
vided a “spectrum parameter’, characteristic of the 
sound to be measured, is added to the sound pres- 
sure level. However, this proposal is certainly not 
rigorous in the case, for example, of pure tones, 
owing to the nature of the equal loudness contours. 
It has previously been shown by the author [5] and 
is confirmed by experiments described in this paper 
that the scale relation between loudness and loud- 
ness level under fixed conditions of test is indepen- 
dent of the sound material, whereas the relation 
between loudness and sound pressure level is not. 
This result is a logical consequence of the experi- 
mental fact that loudness is a single-valued and 
measurable attribute of sounds. Fortunately, the re- 
lations between loudness level and sound pressure 
level are known with considerable accuracy for most 
of the test sounds used. 


Some investigations present results for various 
test ratios, and with one exception the consistency 
within each investigation is satisfactory in this re- 
spect. In order to combine results for different ra- 
tios, use has been made of the fact that the varia- 
tion of the phon interval for a given ratio is a 
relatively slow-moving function of the absolute level, 
in other words, the relation between loudness S and 
loudness level P may be expressed, for moderate 
ranges of level, by the equation S=exp(aP+)), 
where a and b are constants. The results of tests. 
can thus be normalised without appreciable error 
in terms of an equivalent constant ratio, which is 
taken for convenience and without loss of generality 
to be the two-fold ratio. Thus, if a starting level 
of P phons requires an increment or decrement of 
AP phons to satisfy the criterion of r-fold loudness 
ratio, the equivalent two-fold interval is taken to be 
AP/logsr. 

Experimenters have found that the results for 
different test ratios are consistent in each of the 
domains r>1 and r<1 but not in both. This is 
clearly demonstrated by some of the data shown in 
Fig. 1, where, for convenience of simultaneous pre- 
sentation, incremental values 4,P are associated 
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Fig. 1. Experimental data. 


(a) Lairp, Taytor and Witie 
(b) Gricer and Firestone 

(c) Ham and Parkinson 

(d) Rscuevern and Rasinovircn 


(f) Garner 
(g) Porrack 
(h) Rosrnson 


with the centre of the range, namely P + ~AsP , 
and similarly decremental values Aj/2? with the 
level P—FAypP, the subscripts denoting the loud- 
ness ratio and P the initial level. 


It is found that the original data from two of the 
investigations become substantially more uniform 


(e) Cxurcuer, Kine and Davies 


(j) Povurron and Srevens 
(k) Quirrzscu 
(1) Srevens, Rocers and HerrnsteIn 


when an allowance is made for a well-known effect 
associated with the order of presentation of sound 
stimuli. Briefly the effect is that, of a pair of sounds 
heard consecutively, the second appears relatively 
louder than it would if the sounds were presented 
in the other order. The “order effect’’ manifests it- 
self directly in equal-loudness judgments and can 
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be eliminated by presenting the comparison sounds 
in random order. In loudness ratio estimations, on 
the other hand, the sounds are necessarily presented 
in the order indicated by the terms of the ratio and 
in consequence the operation of the “order effect” 
can only be indirectly inferred. Experimental data 
indicating the magnitude of the effect are presented 
later. 


Fig. 1 shows some of the stages of reduction of 
the data to a consistent basis. The convention has 
been adopted of dotted lines connecting points re- 
presenting fractional judgments and solid lines those 
representing multiple judgments. No significance 
attaches to the use of straight line segments joining 
successive points in a series except as an aid to 
clarity. For descriptions of the test procedure re- 
ference may be made to the original papers. The 
abscissa scale is to be interpreted as loudness level 
referred to the 1000 c/s tone. For the free-field de- 
terminations (iii) and (viii) below, it becomes the 
phon scale. 


(i) Larrp, Taytor and Witte [6], Fig. 1a 


Data for loudness halving of tones of 64, 256, 1024 
and 4096 c/s have been included, all except 1024 c/s 
having been corrected for equal-loudness contours ”. 


(i) Gercer and Firestone [8] 


Values for 1000 and 60c/s tones are given, the 
latter after allowing for the equal-loudness relation. 
The data for the 40-component complex tone are omit- 
ted owing to lack of information about the equivalent 
loudness. Fig. 1 b shows that the equivalent two-fold 
interval — Ajj2P or + A2P is systematically smaller the 
more remote the loudness ratio from equality; this is 
the only case in which different ratios show appreciable 
inconsistency. 


(iii) Ham and Parkinson [9] 


Figs. 1c; and 1 cs show the original data for Ham’s 
and Parkinson’s tests normalised to equivalent values 
for a 2-fold ratio. The Fig. 1 c, refers to their experi- 
ments on direct loudness magnitude estimation in terms 
of a standard assigned the value 100. Values for 350 
and 2500 c/s are included after correction for the equal 
loudness relations. Fig. 1 cy shows the same data after 
application of a correction for “order effect”, from 
which it is seen that the consistency of the data in each 
sequence is substantially improved. 


(iv) Rscwevkrn and Rasrnovircu [10] 


Fig. 1d contains the normalised data of Rscuevxin 
and Rasinovircu for ratios of 1/2, 2, 1/4 and 4 and illu- 
strates the consistency of fractional and multiple judg- 
ments within themselves but not with each other. 
> In this and similar cases, the pure-tone equal- 
loudness relations by Rosinson and Dapson [7] were 
used. 
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(v) Cxurcuer, Kine and Davies [11] 


Vol. 7 (1957). 


Fig. le shows their data for ratios of 1/2 and ~ 


1/4, after conversion of the original values for 800 e/s 
to the equivalent loudness at 1000 c/s. 


(vi) Garner [3], [12] 


Fig. 1 £ shows the results of experiments by Garner — 
in 1952 and 1954 on loudness halving, the mean of © 


which is regarded as typifying Garner’s results. 


(vii) Potrack [13], [14], [15] 


Potrack has given loudness ratio data for speech, 
white noise and 1000 c/s tone, together with the rele- 
vant equal-loudness relations. The latter have been 
used to reduce the values to equivalent loudness and 
the transformed values are shown in Fig. 1 g; and 1 gs 
for doubling and halving respectively. In spite of some 
discrepancy between the results for the different 
sounds, the averages of the four determinations are 
remarkably smooth and exhibit the characteristic in- 
consistency between doubling and halving. 


(viii) Rosrnson [14] 


Fig. 1h shows the author’s data (1953) together 
with some values obtained for the purpose of the pre- 
sent study. Values for the 50 c/s and 83 c/s tones, the 
engine noises and the white noises are included after 
corrections for the known equal-loudness relations. The 
multiple/fractional disparity is again apparent. 


(ix) Srevens [16] 


Fig. 1 i, shows the original data of Srrvens obtained — 


by the method of direct loudness estimation. The uni- 
formity of these data appears to be improved, as in the 
case of Ham and Parkinson (iii), by the application of 
a correction for the order effect, notably in the removal 
of sharp discontinuities between the first and succeed- 
ing points in the descending sequences, as may be seen 


in Fig. lig. The values obtained by Srrvens’ method — 


of “no-fixed-standard” are included. 


(x) Povutron and Srevens [17] 


They give values only for a white. noise. In the ab- 
sence of specific loudness data for the noise used by 
them, the conversion to equivalent loudness has been 
made using average values from three previous de- 


terminations, namely Porrack [13], Brrrrar [18] and — 


Rosinson [5]. Fig. 1 j shows the data of Pourron and 
Stevens ‘after conversion to equivalent loudness at 


1000 c/s. This investigation again provides a good — 


example of the fractional/multiple anomaly. 


(xi) Qurerzscu [19] 


The determination by Quierzscu was apparently in- 
cidental to another investigation and very few details 


are given. So far as it is possible to read Fig. 7 of 


Quierzscu’s paper, the results of his doubling and haly- 
ing tests are shown in Fig. 1 k. 


(xii) Srevens, Rocers and Herrnstern [20] 


These experiments were based on attempts to repeat 
the values of Larrp, Taytor and Witxe by reproducing 
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as far as possible the same conditions of test. As may 
be seen in Fig. 1a the results of Larrp, Taytor and 
Witie are somewhat out of line. Srevens, Rogers and 
HerrnstEIn were partially successful in repeating the 
experiment of Lairp, Taytor and Wut, the data 
shown in the upper curve of Fig. 11 being obtained 
when the conditions of Lairp, Taytor and Witte were 
closely copied. Under other conditions, Srnvens, Rocrrs 
and Herrnstrin obtained the other two sets of results 
‘seen in Fig.1]. Apart from a scale factor, it would 
appear that the three sets are consistent in trend and 
for present purposes the average of the determinations, 
shown by the filled-in squares, has been adopted. 


3. Preference for central levels of the 
variable stimulus 


The persistent occurrence of the fractional/mul- 
tiple anomaly has been noted by various authors 
but no convincing explanation has yet been offered. 
It operates in a manner consistent with the assump- 
tion that a true observation of loudness level change 
(i. e. a change differing only in sign for incremental 
and decremental values over the same range of sub- 
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jective loudness) is contaminated by a secondary 
effect inseparable from the usual test procedure and 
is consequently not directly observed. The secondary 
effect may be interpreted as a preference for central 
levels or as a subjective resistance to remote levels 
of the adjustable stimulus. It is interesting to note 
that in Stevens’ absolute loudness estimations em- 
ploying the unusual method of “no fixed standard”, 
the results are intermediate between those obtained 
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Fig. 2. Magnitude of the “centring correction”; curve a: 
decremental loudness judgments, curve b: incre- 
mental loudness judgments. 
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222 


in incremental and decremental situations relative 
to standard reference levels, as illustrated in Fig. 1 ig. 
Examination of the remaining data shows that the 
level of the adjustable stimulus satisfying the cri- 
terion of test may be regarded as being pulled to- 
wards a fixed value of 85 phons by an amount 
which is roughly proportional to the departure 
therefrom, as shown in Fig. 2. The fractional/mul- 
tiple anomaly may be described as a “centring er- 
ror’. It is not the same as the “order effect” and 
apparently operates independently of it. The “order 
effect” tends to inflate all incremental values of AP 
and to reduce decremental values, whereas the “cen- 
tring effect”? operates selectively, depending on the 
absolute level. 


The corrections shown in Fig. 2 have been ap- 
plied throughout the data described, in order to 
present the values in a unified manner. This pro- 
cedure does not, of course, explain away the effect, 
but permits a rational amalgamation of those de- 
terminations which contain values for either halving 
only, doubling only or a mixture. The various data, 
freed from the “centring error”, are presented in 
Figs. 3a—1. 


4. Functional dependence of the increment 
on the loudness level 


The data assembled in Fig. 3 enable some general 
deductions to be made about the shape of the func- 
tion connecting P and AP. The investigations cover- 
ing a wide intensity range, especially Garner, Pot- 
Lack and Rogrnson, suggest that a maximum value 
of AP occurs at around 55 phons, and there is some 
evidence of a minimum at around 90 phons. Other 
investigations covering a more limited range, notably 
Rscuevkin and Rasrnovircu, Poutron and SrTevens, 
and Lairp, Taytor and Witte, show a marked con- 
formity with these trends. Substantial differences in 
the magnitude of AP, however, persist. As will be 
seen later, these would be expected to appear in the 
form of variations in the ordinate values of the 
component parts of Fig. 3 while the functional de- 
pendence on loudness level would be similar through- 
out. In order to illustrate this, the values have been 
fitted in each case by smooth curves (shown as 
broken lines). These curves are all alike except for 
a multiplying factor on the respective ordinate sca- 
les, and a small adjustment in the position of the zero 
on the abscissa scale. The choice of the multiplying 
factors, the horizontal shifts, and the dependence on 
loudness level were governed by the fit of the 
curves to the data points shown, using a process of 
successive approximation. The curves represent the 
data very well in some cases, e. g., Larrp, Taytor 
and Wiiie, Rscuevxin and Rasrnovircu, Portack, 
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Rosrnson, and Stevens, and in the others, except 
for Cuurcuer, Kine and Davis, they are as good 
a fit as the internal consistency of the individual 
investigations would demand. It is interesting to 
note that the fitting of the curves led naturally to 
horizontal scale shifts consistent with the mode of 
listening ?. The shape of the curves is well defined 
over the central portion but the terminal portions — 
are probably less reliable owing to the relative 
scarcity of data below 20 and above 110 phons. 

The multiplying constants range from 1 to 2.7, 
the data of Larrp, Taytor and WIE requiring the 
largest factor. There appear, however, to be grounds 
[20] for excluding the values of Lamp, Taytor 
and Wie, and if this is done the variations 
amongst the remaining 11 investigations are re- 
duced to a range of 2 to | in the multiplying fac- 
tors. The cause of the variations within this range 
and the estimation of the most probable value are 
considered in later sections. 


5. Some new experimental determinations 


The various experiments described below were 
carried out with a view to providing information 
on the effects of alterations in the test procedure. 


5.1. Influence of silent interval on the observed 


value of AP 


Method: C.S.*4 with randomised levels of the second 
tone; free-field listening; halving judgments for levels 
of the fixed tone 110, 70 and 40 phons and doubling 
judgments for levels of 90, 60, and 30 phons, a sepa- 
rate test session being devoted to each of the 6 refe- 
rence levels. The period of silence between the reference 
and test tones was changed at random during the 
course of each test, generally avoiding consecutive 
pairs having the same time sequence. The mark period 
for both tones, however, was always 1.5s. A clickless 
switch of 0.15s time constant was used, and space ~ 
periods were measured by the time interval between 
commencements of decay and build-up. 

The results, which are shown in Fig. 4, were un- 
expected. It was found that values of — AiP for halv- 
ing decreased as the time interval increased, while 
values of A2P increased, but to a lesser extent. The 


5 Relative to the position of the curves fitting the 
data of Ham and Parxinson and Rosinson, for which 
the abscissa scale is phons, the average shift for Pout- 
ton and Stevens, Gercer and Firestone, and Cuurcuer, 
Kine and Davies was 3 phons to the right and that for 
Larrp, Taytor and Wirie, Rscuevkin and Rasryoviren, 
Garner, Poxtack, and Stevens, Rocers and HerrnsTein 
was 5 phons to the right. These groups used binaural 
and monaural earphone listening respectively and the 
shifts are in the direction consistent with the fact that 
values in phons are less than the loudness levels for 
earphone listening. 


4 Constant stimulus method. See reference [5] for 
a fuller description of the test procedure. 
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ness ratio tests. 


average variations were significant in a statistical sense 
and were well-marked and progressive with the time 
interval for many of the individual subjects. The ob- 
servers generally remarked that they were unaware of 
any shift of judgment, though some found that the 
shortest and longest periods (0.3 and 3.0) tended to 
upset the rhythm of their judgments. Consideration of 
the trend of the results and a search for the possible 
origin of the observed variations prompted a further 
series of experiments (5.2). 


5.2. Variations of the order effect with silent inter- 
val and loudness level 


As has been seen, the order effect in loudness ratio 
tests may only be demonstrated indirectly, but inferen- 
ces of its occurrence may be drawn from the manner 
in which the order effect varies with the silent interval 
and with the loudness level. 

Direct estimation of the effect as a function of loud- 
ness level for the constant space period of 1.0 is af- 
forded by numerous data obtained by Rosryson and 
Dapson [7] in the course of an investigation of the 
equal-loudness relations for pure tones. The average 
values at intervals of 10 phons in the loudness level 
are shown in Fig. 5a. In allowing for the order effect 
in relation to data by Ham and Parktnson, and Stevens, 
various trial values were assumed. Those leading to 
the optimum improyement in uniformity are shown by 
the solid line in Fig. 5 a, which can be seen to follow 
the directly determined values fairly closely. 

Some further equal-loudness tests by the constant- 
stimulus method were carried out to determine whether 
the order effect depended on the silent interval. The 
tones compared were (i) 1500 c/s at 75 dB and (ii) 
1000 c/s at 80 dB, the reference tone being a 1000 c/s 
tone in each case. Listening was in free-field, the mark 
period of each tone was 1.0 and the space period was 
set at random to 0.3, 1.0 or 3.0s. The level of the re- 
ference tone and the order of presentation of the test 
and reference tones were also varied at random from 
judgment to judgment. The results are shown in Fig. 
5b, supplemented by some earlier data by the author 
for tones of 100 and 1000 c/s. It will be noted that the 
magnitude of the order effect diminishes as the silent 
period increases. It is greatest in the heterophonic com- 
parison (1500/1000 c/s) but is still appreciable in the 
comparison of the 1000 c/s tone with itself. Neither of 
these situations corresponds exactly to the conditions 
of a loudness ratio test where the sounds, though of 
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Fig. 5. (a) Variation of “order effect’ with loudness 
level in constant stimulus equal-loudness 
tests with 1.0 second silent interval. 
— curve used in applying correction 

to loudness-scaling data of Ham 
and Parkinson [9] and Srevens 
[16]; 

¢ Values due to Rosinson and Dap- 
son, based on equal-loudness tests 
with 1.0 second silent interval. 

(b) Variation of “order effect’? with silent in- 
terval. 

¢ 1000/1500 c/s, 28 subjects, 75 dB, 
x 1000/1000 c/s, 23 subjects, 80 dB, 
* 1000/100 c/s, 31 subjects, 85 dB. 


the same character, are distinct in loudness. It seems 
reasonable, however, by comparing Figs. 4 and 5b to 
conclude that here are two manifestations of the same 
effect. Thus, the variations in the loudness ratio tests 
observed with varying time interval are probably ex- 
plained by the order effect, i.e. by an apparent change 
in the loudness of the second sound, and there is no 
evidence of a change in the ratio judgment itself. 


5.3. Comparison of the methods of adjustment and 
constant stimuli 


Reasons for preferring C.S. technique to the “Method 
of adjustment” have been argued in another place [5] 
and up till now all the author’s data were obtained 
in this way. The method of adjustment (M.A.) has, 
however, been widely employed in various forms by 
other workers and it seemed desirable to demonstrate 
that the procedures are essentially equivalent even if 
not equally rapid or reliable. In the C.S. method, the 
observer has no control over the presentation of the 
stimuli. M.A., on the other hand, gives the observer 
control over the intensity of one of the tones by means 
of a knob covering an adequate range. Control of the 
order and durations of the stimuli may alo be given 
to the observer (M.A. 1) or reserved to the experimen- 
ter (M.A. 2). 

The method of adjustment, however, introduces a 
variable which does not occur in the randomised C.S. 
procedure, namely the gain/rotation characteristics of 
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the observer’s knob. A series of M.A. tests was there- 
fore carried out by 26 observers, their task being to 


determine the half-loudness level of a free-field 1000 , 


c/s reference tone of 70 phons. Three knobs were pro- 
vided having widely differing characteristics but each 
covering a total range of 40 dB: 


a) loudness roughly proportional to rotation, 
b) uniform decibel steps, 
c) a linear series rheostat. 


Each observer made three judgments with each of the 
three knobs in cyclic order, first by the M.A. 1 pro- 
cedure, and then on a separate occasion by the M.A. 2 
procedure. The values obtained are shown in Table if 
each entry being the median of 78 observations. No 
appreciable shift in the half-loudness level was found 
even in the case of knob c) which was very ill-condi- 
tioned in the neighbourhood of the typical setting. 


Table I 


Interval for halving 


Knob (median value) 
MALI | M.A. 2 
ooo 
a) 10.1 phon | 10.7 phon 
b) | 10.7 phon 11.4 phon 
c) | 10.2 phon 10.5 phon 


Experiments were also carried out to provide a di- 
rect comparison with the results obtained by M.A. and 
C.S. procedure. Three tasks in the order below were 
given to 26 observers, and these were carried out at 
a single session. 

a) “doubling” a free-field 1000 c/s reference tone of 
60 phons by the M.A.1 procedure. The observers 
had a 3-position clickless switch giving the reference 
tone, the adjustable tone and silence respectively, 
and a knob giving loudness roughly proportional 
to rotation. Three adjustments were made for ran- 
dom “bias” settings of an independent attenuator 
set by the experimenter. 

b) as above but using the M.A.2 procedure with a 
fixed sequence of presentation giving the reference 
and test tones alternately for 1.5s separated by 
silent periods of 1s. 

c) the same judgment carried out by the C.S. method 
with randomised levels. This was a control test, the 
same observers having already carried out the same 
determination in tests unconnected with M.A. 

In a second test session, a similar procedure was 
followed for “halving” a 1000 c/s tone of 70 phons, 
except that on this occasion only one adjustment was 
made under a) and b). 

The results, shown in Table II, require careful inter- 
pretation, and should be taken in conjunction with 
those of the halving tests with the three knobs des- 
cribed above. If the latter results are averaged for the 
three knobs, values of 10.4 and 10.9 are obtained by 
M.A. 1 and M.A. 2 respectively, whereas the correspond- 
ing values in the present tests (Table II) are 9.9 and 
8.9. Analysis of the M.A. results disclosed a small but 
constant trend. The first in any set of similar adjust- 
ments gave a smaller value of —AjjP or a larger 
value of 42P than succeeding adjustments. 

This is shown in Table II in which the results are 
divided up in order of performance. It would appear 
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Table II 
Method ___ Median value of A P 
cas Doubling Halving 
M.A.1 | 11.9 phon* 9.9 phon* 
M.A. 2 | .14.7 phon* | 8.9 phon* 
C.S. | 14.4 phon 11.1 phon 
Value expected by 
C.S. from previous | 14.5 phon 11.9 phon 
data by same team | 


* Mean of three adjustments + Single adjustment 


Ste 
~ 


that-single adjustments should not be relied upon in 
M.A. The values 9.9 phons and 8.9 phons in Table Il 
are therefore probably subject to an upward correction 
of the order 1 phon. Taking this into account, there 
does not appear to be any systematic difference be- 
tween the C.S. and M.A. results. The author is, how- 
ever, confirmed in his preference for the C.S. method. 


Table III 
poe [7 nda aes Drift 
| adjustment | adjustment | adjustment 
Doubl- | 
ing | 13.9 13:1 —1l 


Ist Average of 9th 
adjustment | 7adjustm. | adjustment 
Halv- | ; oF} papery Ey , 
ing 9.9 10.7 11.1 1.2 


5.4. Ratio tests with low frequency tones 


The following experiments were carried out to ob- 
tain confirmation that the changes of level for a given 
loudness ratio are independent of the type of sound 
when expressed in equivalent loudness units but not 
when expressed in decibels. Three tests with pure tones 
were carried out by the C.S. method, by 26 observers, 
at separate test sessions at intervals of some months. 
The time sequence used was: mark 2.0s, space 1.05, 
in each case, and listening was in free-field. Low fre- 
quency pure tones were chosen as the test sounds be- 
cause the slopes of their equivalent loudness/sound 
pressure level relations depart substantially from unity, 
and are moreover known with considerable precision 


[7]. 


Table IV 
1 Peat 3 Le Od, 5 6 
; 8 ‘| splot | Cha 
2 | standard Loudness of s. —- Shae 
requency donee level of oatnnlne interval 
of sound Ratio | (@B rel | standard | (™ for 2- 
(c/s) | 9.0002 tone paper fold 
| dynfem’ | (hom) | “(qpy | tao ae 
3.3 | x% | 738 | 706 | 85 | 102 
1000 | xe | 19.6. 1 708 31.6) “PRS 
83.3 | x2 55.3 | 41.0 | 198) | eee 
1000 | x 2 39.6 39.6 | 17.2 | W712 
50 | xy% 74.6 | 55.0 7.6 10.8 
1000 Ue a 1 2:6 52.6 11.6 11.6 
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The results are shown in Table IV. Alternate entries 
are data for the 1000 c/s tone obtained in 1953 by the 
same method, at roughly corresponding loudness levels, 
for a group consisting mainly of the same observers. 
Comparison of pairs of entries in column 5 for the 
1000 c/s tone and the low frequency tones shows a 
considerable discrepancy in the decibel intervals, which 
is almost eliminated by the conversion to phon inter- 
vals, as shown in column 6. 


5.5, The comparison of the loudness of sounds se- 
parated by a long time interval 


In practical applications of the loudness scale, such 
as noise abatement, very long times may elapse between 
the sounds to be compared, and it is of interest to 
determine whether the scale relation determined with 
brief time lapses of the order of seconds maintains its 
validity over greater periods. In the present experiment 
28 observers listened individually to’ a free-field re- 
ference tone of 1000 c/s at 90 phons for as long as 
they desired. Some five hours later they again listened 
to the 1000 c/s tone at a modified level and were asked 
to judge its loudness in relation to that heard before 
(14 of the observers heard 100 phons on the second 
occasion and the remainder 80 phons). Estimates in 
terms of a standard of 100 ranged from 90 to 200 for 
the first group (only one observer judging the change 
in the wrong direction), and from 95 down to 25 for 
the other group with no inconsistencies. Interpreting 
the combined results, the median observation was equi- 
valent to a value of 15.1 phons for a 2-fold loudness 
change, with a standard deviation of 4.2 phons. It is 
encouraging to find that this experiment leads to a 
value comparable with those obtained under more 
usual laboratory conditions. 


5.6. The employment of unpractised observers 


Qurerzscu [19] has argued that the only loudness 
ratio which can really be comprehended by listeners 
is 2:1, and that tests for other ratios are made in 
terms of doubling or halving with appropriate adjust- 
ments. It does not seem likely to the author that ob- 
servers untrained in acoustics or mathematics would 
in fact proceed in this tortuous way if set a ratio task 
other than 2: 1, but would use a more direct correla- 
tion with the cardinal numbers. 

A group of unpractised subjects without previous 
experience of subjective tests or psycho-acoustics were 
set the tasks of (i) “multiplying” 1000 c/s tone by 10, 
(ii) “doubling” and (iii) “halving”, in that order, in 
a single test session by C.S. procedure. The results are 
shown in Table V and suggest that initiation into loud- 
ness ratio testing by way of a more “difficult” ratio 
does not materially affect the result. The consistency 
of the 2-fold and 10-fold values should be noted. 


Table V 


Phon Interval 4P 


Level of —_ = a 
standard Mean value Value expected 
tone for 10 | by interpolation 
(phon) | unpractised of author’s 
| observers | 1953 data 
| } / 
45 emietor ah 403: | 46.9 
60 x 2 16.6 16.0 
70 | x% 9.7 | 12.0 
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None of the variations of test procedure intro- 
duced in the experiments just described brought 
about changes in the judgments of observers com 
parable in magnitude with the overall discrepancies 
between different observers or different investiga- 
tions. To be conclusive, a very large investigation 
employing a fresh group of observers with each 
variant of the test procedure would be required, but 
from the results already obtained it seems more 
likely that an explanation of the differences is to be 
sought in another direction. 


6. Statistical considerations 


6.1. The problem of weighting 


It has been shown by Garner [12] and Rosin- 
son [5] that the variance in loudness ratio tests is 
dominated by the component due to differences 
between observers, while the error variance as esti- 
mated from repeated results is comparatively small, 
at least, under constant experimental conditions. The 
author has further demonstrated [6] an extremely 
high “concordance coefficient’ between different 
observers’ data as a function of intensity level, a 
result which indicates that there is general unifor- 
mity as to the form of the functional dependence 
of AP on P though not necessarily as to the abso- 
lute values of AP. Both these observations suggest 
that it is correct to weight the investigations in pro- 
portion to the number of observers taking part. 
Comparatively little is added by the repetition of 
tests, nor does the varying density of the data along 
the intensity scale make much difference. 


~~ 


The question of averaging the data arises at three 
stages, beginning with the conduct of individual 
tests, next, in the representation of the data for the 
various groups, and finally in relation to the data 
as a whole. 


6.2. The problem of averaging 


Of the various factors which may affect the result 
of an individual test, one, the centring error, has 
already been allowed for. Another is the order effect. 
It has been shown how an allowance for this aids 
the presentation of some of the data. In other cases, 
however, its effect on the data would be equal and 
opposite in respect of incremental or decremental 
changes of loudness. It is equivalent to a translation 
of the correction curves for the centring effect, de- 
creasing somewhat the loudness level at which they 
would intersect. The result, so far as the fit of the 
curves to the data points of Fig. 3 is concerned, would 
be negligible. Another factor concerns the spacing of 
the levels of the variable stimulus. As has been 
shown, strong variations in the “feel” of the control 
in MA make little difference to the result obtained. 
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Equally, in C.S. tests with random levels, it is im- 
material what levels are used provided an adequate 


range is covered in random steps sufficiently nu- ’ 


merous. In determining the transitional level from 
the test record, however, it is in principle necessary 
to take account of the scale in which the judgments 
are recorded. The range of uncertainty in practice, 
however, is only of the order 3 dB for a 2-fold test, 
and frequently less. Consequently very little error 
will result from taking the centre of the uncertainty 
range in whatever units it is expressed. Allowing 
for these factors, the conversions to a common ratio, 
to equivalent loudness units and to-a common zero 
for the scale of loudness level, it appears that one 
has gone as far as possible towards the elimination 
of systematic effects from the original data. 


Various statistics have been used to specify the 
“average” values for a group. Ham and Parkinson, 
Rscurevxin and Rasrnovircu, GARNER and PoLiack 
use the arithmetic mean of decibel values, while 
Cuurcuer, Kine and Davies, Gercer and Firestone, 
Rosinson, Poutron and Stevens, and Stevens, Ro- 
cers and HerrnsTEIn give both mean and median. 
For loudness estimation type tests, Ham and Par- 
KINSON give the arithmetic mean of loudness values, 
while Srrvens gives only the median. The data of 
QurerzscH are treated in an unspecified manner 
while Larrp, Taytor and Witte alone give original 
data. In view of the relatively large dispersions, the 
units in which the values are averaged becomes im- 
portant. For example, distributions of AP values are 
usually skewed towards the high values. If the values 
were ie an in linear sound pressure units the 
distribuffons would be entirely different. In terms 
of sound power units they would be different again. 
Since objective scales of the strength of sounds take 
no account of the characteristics of hearing, the use 
of any one is as arbitrary as the use of any other. 
It seems clear therefore that the averaging of group 
data should take account of the characteristics of 
hearing, even if somewhat conventionally, in order 
to define a unique value. One such proposal is em- 
bodied in a method recommended for some pur- 
poses by Stevens [21]. This consists of an iterative 
process by which the observed values are converted 
into loudness units using a trial function, averaged 
arithmetically, and the resulting average used to 
estimate a second approximation to the function, 
and so on. The procedure, however, is_ strictly 
applicable only to repeated results by a single ob- 
server, since it assumes a common loudness function 
for all the observers, a proposition which is not ex- 
perimentally confirmed and is, of course, inherently 
unobservable as regards absolute scale numbers. 
Moreover, for a reason which will appear later, 
Stevens’ method lacks generality and this restricts 
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its usefulness. If applied, for example, to identical 
sets of AP values, one representing “doubling” and 
the other “halving” intervals, Stevens’ procedure 
leads to two distinct average values. This would imply 
that the group values for “doubling” and “halving” 
were not consistent although the values for every one 
of the individuals were. It is useful to consider what 
general principles should guide the choice of averag- 
ing method, or even, in fact, whether what is re- 
quired is an averaging process at all. The use of 
averages is properly applied to a set of values con- 
sisting of a common “true”’ value contaminated with 
variable errors, and the averaging process is aimed 
at separating the true value from the error. If, how- 
ever, one has a set of numbers of which it is known 
that the errors are less than the deviations from a 
central value, a different situation exists. Averaging 
such a set of inherently different numbers is neces- 
sarily an arbitrary process; a more logical proce- 
dure is to nominate a typical member of the set, for 
example the modal or the median value. The con- 
clusions already arrived at in this paper might sug- 
gest that this is the situation in regard to loudness 
scales for different observers, but another, interpre- 
tation is possible, namely that unpractised observers 
start with an equal probability of arriving at the 
same judgment but owing to the imprecision of the 
subjective estimation of ratios this common judg- 
ment is subject to a random error. Judgments once 
arrived at, however, seem to be adhered to by the 
individuals, and thereafter may be regarded as a 
set of inherently different numbers in the sense dis- 
cussed above. The median value is evidently to be 
preferred on any interpretation as this has the pro- 
perty of invariance under all monotonic transforma- 
tions of the variate. A knowledge of the median 
value alone, however, is insufficient to permit an 
interpretation of the causes of the variance. 

Similar considerations apply to the results of the 
investigations considered as a whole, but the deter- 
mination of the grand median ideally requires a 
knowledge of all the individual values, and these are 
not available. Statistical information is lacking or 
scanty in some cases. Fortunately Geiger and Fire- 
stong, Larrp, Taytor and Wiiie, Garner, and 
Povutton and Srevens give sufficient data, which, 
together with the author’s original values, enable a 
reasonable estimate of the shape of the frequency 
distribution occurring in loudness ratio tests to be 
made, from which some interesting conclusions 
emerge. 


Table VI shows five determinations chosen as an 
example of the procedure. It fortunately happens 
that these five sets are not only comparable with 
each other as to ratio and absolute loudness level, 
but are also typical of their respective investigations 
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Table VI 
Experimenter | Betis Source Reference Number Statistical data 
| of sound level of observers available 
Larrp, Taytor and WILLE | xy 1024 c/s 60 dB | 10 Original values 
GEIGER and Firestone . xy, 1000 ¢/s 55 dB | 44 Mean, median, max. 
and min. values, 
| average deviation 
ROBINSON | xh 1000 ¢/s 53 phon 25 Original values 
GARNER x, 1000 e/s 60 dB 18 Mean, spread, 
standard deviation 
Povtton and Stevens xh White 50 dB 36 Mean, median, 
noise standard deviation 


(for example, the median of the group of GrrcER 
and Firesrone is 0.94 times the mean value, while 
the average value of this factor for the whole of 
their determinations is 0.92. The corresponding 
figures for Poutron and Srrvens are 0.88 compared 
to 0.90, and those for Rozsryson 0.99. compared to 
0.98). Assuming the distributions to be unimodal 
it was possible to reconstruct the histograms for 
Povutton and Srevens, Garner, and Gercer and 
Firestone in a fairly determinate manner, while 
those for Ropinson, and Larrp, Taytor and WILLE 
are available from the original data. The five distri- 
butions are shown in the torm of smooth curves in 
Fig. 6, the areas under them being proportional to 
the weight of the observations. The distribution of 
the sum has the shape of the dotted curve, and 
shows a fairly large positive skewness. 


“Aq? —~ 


Fig. 6. Reconstructed frequency distributions for “halv- 
ing” experiments at approx. 60 phons. 


Curve a: Geicer and Firestone [8], n=44, 
Curve b: Povrron and Srevens [7], n=36, 
Curve c: Rosinson [14], n=18, 
Curve d: Garner [3], [12], n=18, 
Curve e: Laren, Taytorand Witte [6], n=10, 


Curve f: Pooled data with scale factor adjust- 
ed for >)” =50. 


7. Transformation of the variate scale 


A non-normal distribution of AP values would 
be expected on general grounds. For example, what- 


ever variations may occur between observers, one 
would never obtain both negative and positive va- 
lues in the same experiment, for this could only 
result from a gross confusion of augmentation and 
diminution of loudness relative to a standard. The 
value AP =0 is clearly a limit point. Another limit 
is set by the threshold of hearing or pain depend- 
ing on the direction of the experiment but, except 
for tests at rather extreme reference levels, these 
constitute a practical rather than an inherent re- 
striction of the possible range of AP. The existence 
of an end point of the distribution at JP =0 sug- 
gests the logarithmic transformation x—=log| AP}. 
This relation satisfies the desired criterion of in- 
variance with reversal of the sign of AP and as will 
be seen, it permits of a simple interpretation in 
terms of random errors. Fig. 7 shows the distribu- 
tion curves for the data of Table VI transformed to 
the x-scale and it may be seen that their total is ap- 
proximated, except for the extreme high values, by 
a normal-error curve. By calculating y? it was found 


Relative frequency —= 


0 04 08 1.2 1.6 
X= log |AP1 —= 


Fig. 7. Experimental frequency distribution curyes trans- 
formed into the scale x=log| AP |. 
—-—w Distribution of the pooled data (with 
arbitrary scale reduction), 
eta. t Norma! error distribution; 


Curve a: Geicer and Firestone [8], 
Curve b: Pourron and Srevens [17], 
Curve c: Rosinson [14], 

Curve d: Garner [3], [12], 


Curve e: Larrp, Taytor and Wue [6]. 
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that the data of Grrcrr and Firestonr, Poutron and 
Stevens, and Ropinson could be regarded with rea- 
sonable probability as samplings from the normal- 
error distribution. The probability that Garner’s 
data also represented a sampling from this distri- 
bution could not be excluded (the significance level 
was between 0.02 and 0.01) though it would be 
fair to say that it is not a likely random sample. 
Thus, in spite of the apparently large disparities 
between the determinations, they are not as exclu- 
sive as they look at first sight. Moreover this argu- 
ment is supported by examining the individual di- 
stributions. The Poutron and Srevens distribution 
is positively skewed and lies generally to the left 
whereas that of Garner is negatively skewed and 
lies to the right. This is what would be expected if, 
in sampling from a normal distribution, it so hap- 
pened that the sample mean departed one way or 
the other from the grand mean. It would, however, 
indicate some lack of randomness in the samples, 
and at least one way in which a bias may have en- 
tered is suggested by the following consideration. 
The experimental data do not extend beyond AP 
values of about 28. At the other end a value as 
small as 2.3 is quoted (Geiger and Firestone). A 
corresponding value at the high end would be very 
large in the scale of phon intervals and if obtained 
it would probably have been ruled out as anoma- 
lous. Rscuevkin and Rasrvoyircu actually separated 
their results into those for “normal observers” and 
those for two with “reduced sensitivity” and reject- 
ed the latter. Fig. 8 shows the normal-error distri- 
bution in the x-scale transformed back into loudness 
level intervals JP and shows the tail extending into 
the region of Rscurvxrn’s and Rasrnovirtcn’s rejected 
values (28.5). Thus, these may have been admis- 
sible after all. The author has met with the same 
experience, one or two values exceeding 30 phons 
having been observed on occasions. 


Fig. 8 also illustrates the same distribution trans- 
formed into loudness units (with the arbitrary as- 
sumption that 10 phons corresponds to a 2-fold 
ratio of loudness). Two distributions result depend- 
ing on the sign of AP. In the case of decremental 
values the distribution assumes a form not unlike 
the original Gaussian curve and illustrates the case 
of successful application of Srrvens’ averaging pro- 
cess. However, as he and Poutron remark, when 4P 
is positive the distribution becomes more skewed 
than before. 

It is logical to seek an interpretation of the trans- 
formation x=log!| AP|, which converts the experi- 
mental data to the normal distribution, by enquiring 
into the nature of the “errors” which are so distri- 
buted. To the extent that the logarithm of the loud- 
ness is a linear function of the loudness level, 
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the transformation is formally equivalent to 
«=log|Inr| where r is the task ratio of the sub- 
jective judgment. For if one considers a ratio test 
with r as the ratio and Sp, as the loudness of the 
reference sound, then 


InS)=aP+b, 
In(r So) =a(P+ AP) +6 
whence 
|Inr|=a|AP| or log |Inr|=log | AP | + constant. 


_E| | Positions of the mean values onthe APsscale 


(i Common 
f, / median value 


Relative frequency —= 


2.0 x=log \API—-> 
phon 30.) Pi —= 


Fig. 8. Transformations of a variate, normally distri- 
buted in the z-scale. 

Normal-error distribution with mean 1.0 

and standard deviation 0.2 in the x-scale, 

Same distribution transformed into the 

scale of AP, 

Same distribution transformed into a 

scale of loudness numbers S with log S 

=0.0301P+K, for dP<0 and 

AP > 0 respectively. 


Scales are arbitrarily aligned at the common 
median value corresponding to 2-fold ratio. 


The exponential relation is not exactly true nor 
is the experimental distribution exactly normal in 
terms of log AP, but neither of these approxima- 
tions necessarily excludes the possibility that 
log|Inr| is distributed according to the normal- 
error law, and it is interesting to consider the im- 
plications of such a relation. 


Let us assume that the error in subjective assess- 
ments of loudness ratio is mainly confined to er- 
rors in the introspective judgments of the numbers 
assigned to the ratio, whilst the physiological fae- 
tors are fairly constant for different observers. Apart 
from secondary effects it is experimentally observed 
that values of AP for individual observers differ 
only in sign for a given task ratio and its recipro- 
cal, and the values AP for different values of r are 
equally consistent with each other whether one con- 
siders two ratios larger than unity, say, or their 
respective reciprocals. The term |Inr| takes these 
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factors into account since its value is unchanged if 
1/r is written in place of r. The discrimination of 
loudness ratio which an observer can make is not, 
however, directly measured by a quantity of this 
type as can be seen, for example, by considering the 
pairsh ot ratios (121); (1:2) and (1:10), 
(1:20) which are equidistant as measured by 
|Inr|. Subjectively these pairs of ratios are not 
equidistant. It is undoubtedly easier to distinguish 
in one’s mind between pairs of stimuli satisfying 
the criteria of “equality” and “doubling” than it is 
to distinguish between the remoter ratios even 
though one of these is twice the other. The number 
“20” would have to be much further away from 
“10” for the 4th ratio to be as distinct from the 3rd 
as doubling is from equality. It seems more likely, 
in fact, that sensation ratios are equidistant when 
they stand in the same power relationship, i. e. when 
they are separated by constant values of log |Inr|. 
Measured in this way, the ratios (1:2), (1:4) are 
as distinct from each other as the ratios (1:10), 
(1: 100) and much more so than the ratios 
(1:10), (1:20) of the previous example. If this 
interpretation is correct one concludes that the ex- 
perimentally observed normal distribution in the 
x-scale is due to the operation of random error in 
the judgment of subjective ratios. Supporting evi- 
dence is provided by a considerable amount of data 
on the relative dispersion in tests at different task 
ratios, notably by Ham and Parkinson, Geicer and 
Firestone, and Rosiyson. The argument rests on the 
observed close proportionality between the values 
AP and the standard deviations in AP-units if a 
small and constant component of the observed va- 
riance® values is first subtracted, representing pre- 
sumably the instrumental errors, inherent differen- 
ces in the hearing mechanism of different listeners, 
and the finite value of the differential intensity 
threshold. 


8. Assessment of the most probable form 
of the loudness function 


The foregoing arguments suggest that the experi- 
mental distributions of 4P values are samples of a 
variate normally distributed in the x-scale. A value 
typical of the respective investigations has already 
been determined, namely the multiplying constant re- 
quired to produce the best fit of the broken lines in 
Fig. 3. Under the logarithmic transformation this 
quantity is distributed in the same way as log | AP | 
differing only by an additive constant. The distribu- 
tions already analysed had an r.m.s. standard devia- 
tion of 0.18 in the z-scale, and this value was tenta- 


5 Corresponding to a standard deviation of the order 
2 phons. 
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tively attributed to the remainder of the data. 12 dis- 
tribution curves were then drawn, with areas propor- 
tional to the weight of the observations and median 
values estimated from the logarithm of the multi- 
plying constants with an adjustment in cases where 
the curve-fitting was to mean values. The summed 
distribution closely resembled the normal error form 
with a standard deviation of 0.2, this value con- 
forming with the calculated variance (compounding 
the variance within investigations, 0.187, with that 
between the median values of the investigations, 
namely 0.10?, the overall standard deviation should 
be (0.18? + 0.10?) "= 0.21). It is important to note 
that of the total variance, about 80% occurs within 
the investigations, a fact which puts the apparent 
discrepancies between investigations in a more 
favourable light. 

Over much of the intensity range there are more 
than 250 observations contributing to the total so that 
fairly narrow limits can be set on the median value 
of x in spite of the spread of the original data. A 
final assessment of the data, expressed in terms of 
P and AP, is shown in Fig. 9. Confidence limits are 
shown for the 0.001 significance level, justifying a 
high degree of reliance on the assertion that the 
function lies within the band. Owing to the relative 
scarcity of data at the highest and lowest levels and 
the fact that judgments in these regions are attended 
by special difficulties, the limits are omitted towards 
the ends of the curve. 
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se 
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3s 
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8s 
£2 
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0) 4 120 phon 

Loudness level at centre of 2-fold range —= 
Fig.9. Interval for half or double loudness expressed 


in terms of equivalent loudness in phons. 
———— Value arrived at in this paper, with 
99.9% confidence limits, 
Srevens’ proposed approximation, 
——-— Cnuourcuer 5th power law, S=P*- 1078, 
-xX-X— Former American Standard Z. 24. 2 
(1942) (from the table by Getcer), 
—-—-+— Kvwavss’ formula 
S=1(10~527+1)-28- 10-8, 
Gates’ function. 


Some other forms of the 4P-function which have 
been in vogue at different times are also shown on 
Fig. 9 for comparison. Fig. 10 shows the function 
as it would appear if referred to logarithmic units 
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of the sound pressure of a wide band white noise. 
It is interesting to note that in this form, which is 
perhaps rather less arbitrary than in terms of the 
equivalent loudness levels for the 1000 c/s pure tone, 
the function is somewhat simpler. 


= 
Oo 


(Sal 


Increment or decrement for 
2-fold ratio of loudness 


120 


40 80 
Sound pressure level (dB rel. 0.0002 dyn /em2) —= 


Fig. 10. Interval for half or double loudness expressed 
in terms of sound pressure levels for a wide- 
band white noise. 
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Fig. 11. Relation between loudness and loudness level; 
this paper, 
—-—-— proposal by Srevens. 


The deduced loudness function is shown in Fig. 11 
in the conventional form and the corresponding 
values are given in tabular form in Appendix 2. It 
may be remarked that Srevens’ proposal [16] to 
adopt a constant value of 10 for the change of level 
accompanying a 2-fold loudness change leads to a 
loudness function differing only by +2.5 phons 
from the values arrived at in this paper, at least 
over the greater part of the range. It would, of 
course, be a remarkable coincidence if the true form 
of the AP-function were so simply connected with 
the phon scale in view of the arbitrary definition of 
the latter in terms of a particular pure tone. The 
data presented in this paper indicate fairly defini- 
tely that a 10-phon reduction is subjectively larger 
at a level of 90 phons than at 60 phons. For engi- 
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neering purposes, however, it may well be that the 
advantages of a simple logarithmic relation having 
a convenient mnemonic property should prevail 
against a more complicated expression, If this 
course is adopted it should be recognised as a con- 
vention rather than as a fundamental relation. It 
seems possible that theoreticians and research wor- 
kers may find significance in the shape of the AP-— 
function in relation to the mechanism of hearing. In 
this connection it is interesting to compare the func- 
tion already arrived at with cognate values obtained 
by various indirect methods. 


9. Comparison with the loudness function 
determined by indirect methods 


Fig. 12 shows data obtained by Frercner and 
Munson [1] by the monaural/binaural method and 
by the combination-tone method using two equally 
loud component tones, together with that of Pot- 
tack [13], [15] and of Caussk and Cuavasse [22] 
by the monaural/binaural method for a variety of 
sounds. There is evidently substantial disagreement 
even in this field. Most of the values are lower than 
those corresponding to the AP-function obtained 
subjectively, and it appears improbable on this evi- 
dence that the functions are one and the same; on 
the other hand, the chain curve on Fig. 12, showing 
the unweighted average of these indirect deter- 
minations, is of the same general form as the AP- 
function of Fig. 9. 


oe a 


Increment —~ 


10 30 50 10 90 
Loudness level ——= 
Fig. 12. Values obtained by indirect loudness-sealing 
methods. 
Two-component tone method: ; 
Frercuer and Monson, 125/1000 c/s and 
1000/2000 ¢/s tones. 
Monaural/binaural method: 
* CavssE and Cuavasse, 500, 1000 and 1500 | 
c/s tones, 
Potrack, 1000 c/s tone (1951), 
Pottack, white noise (1951), | 
Pottack, speech (1951), 
Fiercuer and Munson, 125, 1000 and 4000 
c/s tones, 
© Porrack, 1000 c/s tone (1955), 
—-— Unweighted average of seven determi- 
nations. 


xdpo 


Comparison of the ordinate values of these two 
average curves suggests that binaural loudness is 
about 1.6 times the monaural loudness for the same 
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stimulation, instead of double as has sometimes been 
supposed. A factor less than 2 is not, perhaps, un- 
reasonable bearing in mind the possibility of mutual 
inhibition of neural responses. The value of 1.6, 
however, rests on the rather variable data presented 
in Fig. 12. There is one objective factor which 
would tend to produce low values of the loudness 
level increment in the binaural/monaural method, 
namely the inequalities of left- and right-ear acuity 
of the subjects. If, for any individual, this difference 
is large, it will result in a false value of the observed 
increment of binaural over monaural loudness. In 
practice, experimenters usually equate the ears by 
introducing compensating attenuation equal to the 
threshold difference. This may not be exactly equi- 
valent to equalisation at supraliminal levels, but a 
study of the distribution of the difference of left 
and right ear thresholds for normal subjects has 
been made and indicates that even with totally un- 
equalised ears the error from this cause averaged 
over a group of normal subjects would be much 
less than the departure of the monaural/binaural 
results from the subjective values of AP. 


One of the curves on Fig. 9 is derived from the 
function proposed by Gates [23] for use in loud- 
ness computation. Garss’ function was arrived at 
heuristically by comparing the subjective loudness 
of various complex sounds with the values computed 
from a scale assumed only to have the additivity 
property. It is not necessarily to be expected that 
his function will have the ratio property, and in fact 
the equivalent values of AP are considerably smal- 
ler than those obtained subjectively. Nevertheless, 
Gares’ function, like the results of the indirect scal- 
ing methods, has some resemblance to the 4P-func- 
tion, which indicates that the summation of loudness 
is closely connected with, though perhaps not im- 
mediately deducible from the ratio scale of loudness. 
Srrvens [24] also has recently questioned the vali- 
dity of the principle of arithmetical loudness additi- 
vity and has proposed a modified law of addition. 


10. Conclusions 


(i) Twelve previous determinations of the ratio 
loudness scale are compared after reducing the ori- 
ginal data to commensurate terms, the necessary 
steps being: a) conversion of intensity scales into 
loudness levels, b) the reduction of values for tests 
with different loudness ratios in terms of the equi- 
valent values for a common ratio, c) the application 
of a correction for an “order effect’ which tends 
to enhance the subjective loudness of the last-heard 
of a pair of stimuli, d) the application of a common 
correction for a “centring effect” which operates in 
such a way as to inhibit judgments which would place 
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the adjustable stimulus towards the extremities of 
the intensity range. 

(ii) The modified data are well fitted in every 
case but one by similar curves relating the change 
in loudness level for a given loudness ratio to the 
loudness level. The curves have the same functional 
dependence on the loudness level but vary in their 
ordinate scales. In form they show a maximum at 
59 phons and a minimum at 90 phons. 


(iii) The discrepancies between the results of dif- 
ferent investigations are not believed to be explain- 
ed by differences in the test procedures. Several 
experiments introducing such changes gave results 
that differed very little. A statistical study indicates 
rather that the larger part of the variance may be 
attributed to error in the introspective judgment of 
numbers expressing the loudness ratio. 


(iv) The various data are averaged on the basis 
of a logarithmic transformation suggested by the 
statistical analysis. An assessment of the most prob- 
able value is made, the limits on which are fairly 
narrow owing to the large amount of data. It is 
considered that a simplified approximation to the 
result obtained would best meet practical require- 
ments, and this would take the form that a 2-fold 
loudness change is equivalent to an interval of 
10 phons in the loudness level. 


(v) Results of various indirect determinations of 
the loudness scale show a trend similar to the sub- 
jective result. Binaural loudness is estimated to be 
1.6 times the monaural loudness for the same stimu- 
lation. This suggests that the principle of arith- 
metical additivity of component loudnesses on 
which existing methods of loudness computation 
from spectrum analyses are founded may require 
some modification. Further work is in progress at 
the Laboratory which it is hoped will lead to the 
derivation of a loudness function suitable for this 
purpose. 
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Appendix 1 


Influence of experimental error 
in the equal ratio/equal interval scaling method 


The essence of the equal ratio/equal interval pro- 
cedure is the determination of loudness levels which 
divide a given range into equal arithmetic and geo- 
metric intervals of the sensation respectively. Since 
the derivation of a loudness scale depends on the 
relative positions of these levels, the result is 
rather sensitive to small errors. The error will be 
least in the case of the bisecting and geometric 
mean levels and somewhat greater in the case of 
subdivided intervals. If we assume, for simplicity, 
the approximate relation 


InS=aP+b 


in which P is the loudness level, S the loudness, and 
a and b constants, then for a pair of loudness levels, 


P, and P,=P)+4P, 


we have 
So=exp(aPy)+)), 
S,=exp(aP,+adAP +b) =exp(aAP)-exp(aP) +b). 
In the nature of the function assumed, the geo- 
metric mean of the sensation interval, 

VS S; =exp(4 a AP)- exp(a Py +) 
corresponds to the loudness level Py +$4P. The 
bisector of the sensation interval is 


}(Sy-+S,) = $exp(a Py +b) [1 + exp(@ AP) ] 
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and the corresponding loudness level is Py +a AP 
where @ is given by 


exp(a AP) —2exp(aa AP) +1=0. 


If A,P is the loudness level increment correspond- 
ing to a sensation ratio of 2: 1 we have 


a A,P =In2, 


and the equation may be re-expressed in terms of 
the derived increment 4,P in the form 


ZATIAE ot rece) aoe] One 


Fig. 13 shows the solution of this equation in 4,P 
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Fig. 13. Determination of loudness level increment for 
2-fold loudness ratio from observed bisection 
levels; 


—-— Level corresponding to geometric means 
of loudness intervals. 


Table VII 


Relative loudness 


Relative loudness Relative loudness 


Loudness | Loudness Loudness 
level rgeraued — = level Se a ee ee ea level SEEEEEEAEEEY (suetpeanen caro wee ee 
(phon) |» This | Logarithmic (phon) This | Logarithmic (phon) This Logarithmic 
paper | approximat. / paper approximat. paper | approximat. 


28 0.523 0.435 
30 0.604 0.500 
32 0.696 | 0.574 
34 0.796 | 0.660 
36 | 0.909 0.758 
38 | 1.08 0.871 
40 | SedRh7 1.00 
42 1.32 1.15 
44 _ 1.50 1.32 
46 1.69 1.52 
48 1.90 1.74 
50 2.14 2.00 
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for assumed values of AP of 10, 20, 30, 40 phons 
and for values of a running from 0.5 to 1. It may 
be seen, for example, that an error of only 0.5 phons 
in the position of the bisecting level of a 40 phon 
range would change the deduced value of 4,P from 
16.5 to 20 phons. Uncertainties of larger magnitude 
than this have been reported by Garner [3] de- 
pending on the time sequence of presentation of the 
stimuli in equisection tests. Moreover, the propor- 
tionate error would be still greater the smaller the 
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loudness level range, and greater for the 1/4, or = [7] 
3/4 or other subdivided intervals than for the bisec- (3) 
tion point. This factor might be sufficient to account 
for the discrepancy between the subjective and “cor- [9] 
rected” ratios in GARNER’s tests. 
[10] 
Appendix 2 11] 
Table of loudness against loudness level 12] 
Table VIL gives the relative loudness values 153. 
deduced in this paper against the loudness level and [13] 
corresponds to the full-line curve of Fig. 11. It was a 
obtained by numerical integration of the reciprocal [16] 
of the AP function of Fig. 9. Values of the loudness [17] 
computed from the logarithmic approximation, by 
which a change of 10 phons corresponds to a con- [18] 
stant loudness ratio of 2:1, are also shown for 113. 
[19] 
comparison. The unit of loudness defined as the [20] 
sone is the loudness of a sound of loudness level 
40 phons and this convention is adopted for the [21] 
logarithmic approximation. Loudness values for [22] 
os [1942], 405. 
the curvilinear relation have been adjusted by an [23] 
arbitrary multiplying factor for best fit to the log- (discussion). 
arithmic relation. [24] 
(Received 22nd June, 1956.) 807. 
Corrections 


R. Cuocnorre, La sensibilité auditive différentielle d’intensité 
en présence d’un son contralatéral de méme fréquence. (Acustica 


7 [1957], No. 2, pp. 75 —83). 


Les tableaux B et C de la page 79 doivent étre rectifiés comme suit: 


Tableau B: Les colonnes 2 4 4 sont relatives aux essais 4 10 dB au lieu de 40 dB. Les 
colonnes 8 a 10 sont relatives aux essais 4 40 dB au lieu de 20 dB. 

Tableau C: Les colonnes 2 a 4 sont relatives aux essais 4 10 dB au lieu de 40 dB. Les 
colonnes 8 4 10 sont relatives aux essais 4 40 dB au lieu de 10 dB. 
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Summary 

The N-terminal acoustical junctions under discussion consist of a device with one input 
duct and N—1 non-reflecting output ducts, all of which have the same cross-sectional di- 
mensions. A particular group of these junctions with the simplest.possible impedance 
matrix is considered. a 

A general phase relation between the characteristic reflection and transmission factors 
of waves travelling across such loss-free acoustical junctions is derived, and its usefulness 
is demonstrated experimentally for two cases in which V=3. 

In practice the required conditions of symmetry can be fulfilled for values of V=2,3 
and 4 but only approximately when N >4. The method can be applied to certain 2 V- 
terminal electrical circuits and transmission lines, where instances exist in which V > 4. 


Sommaire 

Les jonctions acoustiques 4 VV terminaisons étudiées ici se composent d’un dispositif 
ayant un conduit d’entrée et N—1 conduits de sortie non réfléchissants, ayant tous les 
mémes dimensions de section. On considére un groupe particulier de ces jonctions ayant 
la matrice impédance la plus simple possible. 

On établit une relation générale de phase entre les facteurs caractéristiques de réflexion 
et de transmission d’ondes traversant de telles jonctions acoustiques sans pertes, et on 
démontre expérimentalement leur validité pratique pour deux cas ou N=3. 

Les conditions de symétrie nécessaires peuvent étre satisfaites pour V=2, 3 et 4, mais 
d’une fagon approchée seulement si NV > 4. La méthode peut étre appliquée a certains 
circuits électriques 4 2 N terminaisons et a des lignes de transmission, des cas se pré- 


sentant dans lesquels V > 4. 


Zusammenfassung 


Die hier behandelten /V-fachen akustischen Verzweigungen bestehen aus je einem Ein- 
gangsrohr und V—1 nicht-reflektierenden Ausgangsrohren, die alle die gleichen Quer- 
schnittsabmessungen haben. Es wird eine besondere Gruppe dieser Verzweigungen mit der 
einfachst méglichen Widerstandsmatrix behandelt. 

Eine allgemeine Phasenbeziehung zwischen den charakteristischen Reflexions- und Uber- 
tragungsfaktoren ftir Wellen, die sich entlang solcher verlustfreier akustischer Verzweigun- 
gen ausbreiten, wird hergeleitet und ihre Niitzlichkeit fiir zwei Falle mit V =3 experimen- 


tell dargestellt. 


In der Praxis kénnen die erforderlichen Symmetriebedingungen fiir die Werte V=2, 3 
und 4 erfillt werden; aber nur angendhert, wenn NV >4. Die Methode ist anwendbar 
auf gewisse 2 N-polige elektrische Schaltungen und Ubertragungsleitungen, bei denen auch 


Beispiele mit V > 4 vorkommen. 


1. Introduction 


Little experimental or theoretical work has been 
reported on sound propagation in duct systems which 
are terminated by more than one (non-reflecting) 
output duct; in this paper such an arrangement will 
be called a multi-terminal acoustical junction. 

The purpose of this paper is to show that a ge- 
neral phase relation exists between the characteristic 
reflection and transmission factors of a specific 
group of multi-terminal acoustical junctions. This 
phase theorem which will be derived below and 
demonstrated for two cases with three ducts, is ap- 
plicable in analogous fields, e. g. in the filter theory 
of certain periodic structures of multi-terminal elec- 
trical circuits. 


2. Electrical analogy for a multi-terminal 
acoustical junction 


The N-terminal acoustical junctions which are 
considered in this paper have a number of specific 
features. There is an input duct and V—1 output 
ducts, all of which have hard and rigid walls and 
the same cross-sectional dimensions. Only one mode 
of wave is propagated along the various ducts, e. g. 
the principal mode or plane waves. Such one-dimen- 
sional cases of wave propagation can be treated by 
electrical analogy, e. g. by regarding the input duct 
as corresponding to an electrical transmission line 
and the acoustical junction as corresponding to a 
multi-terminal electrical circuit. Sound pressure and 
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volume velocity (that is, particle velocity multiplied 33’, ..., NN’ on the various output lines by the 
by the cross-sectional area) at the various reference corresponding co-ordinates 7, %3,..., %y as shown 
planes may be regarded as analogous to electrical _ in Fig. La. 


voltage and current respectively at the correspond- 
ing circuit terminals and this analogy will be used 
in this paper. 

Each reference plane of an acoustical duct cor- 
responds, by analogy, to a two-terminal electrical 
circuit and a 2 \-terminal electrical circuit is there- 
fore analogous to an /V-terminal acoustical junction. 
Such a circuit is shown schematically in Fig. 1 a, the 
electrical network between the various pairs of ter- 
minals being omitted. 


~ Thput transmission line 


Fig. 1. (a) 2 N-terminal electrical circuit terminated by 
non-reflecting transmission lines with the 
characteristic impedance Z,. 

Incident waye = 4 E, exp (—j22y/A). 
Outgoing wave = E, exp (+j2%2m/A); 
BDicet Dna. <-> 

(b) 2 N-terminal electrical circuit equivalent to 
(a) with a generated voltage Ey and supply 
and load resistances of value Zo. 


The loss-free transmission lines connected to the 
N pairs of terminals in Fig. 1a have the characte- 
ristic impedance Z, for waves of wavelength 4 being 
propagated along them. The transmission line lead- 
ing to the terminals 1 1’ is called the input line and 
the other V—1 lines at the terminals 2 2’, 33’,..., 
NN’ are called output lines which are non-reflecting. 
The distance from 1 1’ on the input line is spezified 
by the co-ordinate y and the distance from 22’, 


i 


If a wave represented by 
E. : : 
= exp (—j2ay/A) 


(where the time factor exp(jm¢) is included in E, 
and in the following £,,) is incident on the input 
line there will be, in general, a reflected wave at 
the junction forming a standing wave on this line 
and there will be progressive waves 


En exp ( a= j 25 Lm| A) 


with m=2, 3, ..., N on the various output lines. 
A piece of transmission line 1/2 in length at each 
pair of terminals has been included in the 2 /\-ter- 
minal circuit shown in Fig. 1a to ensure that the 
voltages Ey, Hy, ..., Ey at the various terminals 
are free of disturbances from other modes of waves. 
The circuit of Fig. 1a has been changed into the 
equivalent one of Fig. 1 b which is usually more con- 
venient for computing purposes. These two circuits 
are equivalent as far as the voltages and currents 
at the various terminals are concerned, and the cir- 
cuit of Fig. 1 b is obtained from that of Fig. 1 a by 
replacing the non-reflecting output lines by the loads 
Z) and the input line at the terminals 11’ by a 
generator with the generated voltage Hy (open cir- 
cuit voltage of the input line at 11’) and with the 
supply resistance Z) (cf. THEveNtIN’s theorem). 


3. The characteristic factors of multi-terminal 
junctions with the highest possible degree 
of symmetry 


The 2 N-terminal electrical circuit of Fig. 1 a or b 
can be defined in general by its linear relationships 
between the voltages and currents at the various 
terminals 

N 
oven ts (1) 


/ 
_ 


v=] 


TN ee ae he 


‘ 


Eu a 


where 


The coefficients of these V equations form a square 
matrix which is called the impedance matrix of the 
2 N-terminal circuit. 

For linear and passive 2 \-terminal circuits for 
which the reciprocity theorem holds, as in the pre- 
sent cases, it has been shown that the following 
relation is valid [1]: 


Lye = Lr (2) 


with Vz be. 


The impedance matrix of a completely asymmetrical 
linear and passive 2 V-terminals circuit, therefore, 
has N(N+1)/2 independent coefficients. A reduc- 
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tion of the number of independent coefficients occurs 
when the circuit possesses specific qualities of sym- 
metry. If, for instance, the open-circuit input im- 
pedance (the ratio of voltage and current at the 
input terminals with all the other N—1 pairs of 
terminals open circuited) is the same if any one of 
the V pairs of terminals is taken as input, it can be 
shown from eq. (1) that 


Vie anne (3) 


In addition to the property leading to -eq. (3), 
the symmetrical 2 V-terminal circuits with which we 
are concerned in this paper have the property that 
if any one of the NV pairs of terminals is taken as 
input, the output voltages (and currents) at the 
other V—1 pairs of output terminals have all the 
same value if all output terminals are equally load- 
ed, e.g. by Zy. This reduces the number of in- 
dependent coefficients of the impedance matrix to 
two, which may be called Z,,; and Z,,, and eqs. (1) 
are simplified to 


Ey=Ly1, + 2Zyglyt ... +Zyo1y 
Ey =Zyl,+Zylt .-. tZly (4) 
Ey=Zyo1,+Ziolo+ ... +241, In- 


A 2 N-terminal circuit with an impedance matrix 
as shown in eq. (4) will hereafter be called highly- 
symmetrical. 

The characteristic factors of a highly-symmetrical 
2 N-terminal circuit (or its corresponding /V-terminal 
acoustical junction) are specified from the wave 
patterns on the transmission lines of Fig. 1a as 
follows: the characteristic reflection factor 


Rvyy = Avy exp (j ayy) 


is the complex ratio of the voltages (or sound pres- 
sures in the acoustical case) of the reflected and of 
the incident waves at the input terminals 1 1’; the 
characteristic transmission factor 


T (xy) = Buy exp (j Bay) 


is the complex ratio of the voltages (or sound pres- 
sures) of the transmitted wave at the output ter- 
minals 2 2’ (or 33’... or NN’) and of the incident 
wave at 1 1’ if all output lines are non-reflecting. 


The characteristic factors can also be related to 
the voltages shown in Fig. 1 b (cf. [2]). 
From 


E, 
Ey = 52+ Ep, 


(5) 
where E£,/2 is the voltage of the incident wave and 
Ep that of the reflected wave at 1 1’, it follows that 
Ey 


Rwy) =2 E, 


1 (6) 
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and 
E E Ex 
IP ea i el 
(iy) =2 E, D E, 0 E, (7) 
With the impedance ratios 
“i> Z11/Zo and Z49= Z2/Zo (8) 
and the relations 
E,/ —I,=E,/ —~1,= ... =Ey/ —Iy=Zq 9) 
and 
| Ge ea FS (10) 


of Fig. 1b the characteristic factors can be found | 


in terms of z,, and z,) from eqs. (4), (6) and (7); 
(11) 
(2411) [4s +1 + (W — 2) 249] = (NED) aga 


"(244 + 1) [244-2 + UN — 2). 2) 


oy 242 


(12) 


4. Magnitude and phase relation of the characte- 
ristic factors of loss-free highly-symmetrical 
junctions 

A highly-symmetrical 2 \-terminal circuit which 
has no resistances (excluding the various load re- 
sistances) and therefore no consumption of energy 
within the network connecting the NV pairs of ter- 
minals is called loss-free, and is characterized by 
24, andz,y being purely imaginary [1]: 


Re(z,,)=0 and Re(zy.) =0 (13) 


(where Re(w) means the real part of w). Applying 
this to equations (11) and (12) it can be shown 
that there exist two fundamental relations between 
the magnitude and phase values of R,y) and Ty) 
for any loss-free highly-symmetrical 2 N-terminal 
circuit. The first relation is 


Aiyy + (N—1) Biyy=1. (14) 
This is fairly obvious from energy considerations, 
as Atyy is proportional to the ratio of reflected to 
incident energy at the input terminals, and B7y, is 
proportional to the ratio of transmitted energy at 
each output terminal to the incident energy at the 
input terminals. 
In addition, it yields: 
oe mg a expLj (a) — Bay) 
Zi + (N — 2) 244 249 — (N- 1) 242" — ay 
- Seer = ie 


(15) 


N—2 
eerie 
Since the first term on the right hand side is purely 
imaginary (eq. (13)), eq. (15) becomes 


es? 3 


(4, +1) 4 +1 + (N—2) 29] —(N—1) ae? | 
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bine, { N= 2 
Ree expla ~ Boo)! + Fae (10) 
(N 


) 


and this is equivalent to 


A N-2 
B, , (08 (20 Roto. =o (17) 
(N) 
or 
N—2 Bw 
aw) — By) =a + are cos 5 == ie enGl/-a)) 


This is the other relation between the phase values 
of loss-free highly-symmetrical 2 V-terminal circuits. 
Eq. (14) and (17) may be rewritten: 


Ry) = ee aN et) Beny x (18) 


e N-2 ; Bw) : 
Rep Buy) Hare cos 2 (l-W-1) Ba, 


which shows that R:y) can be expressed in terms of 
B yy and fy) or by the magnitude and phase values 
of Tw) . 

For loss-free and symmetrical four-terminal cir- 
cuits (or corresponding discontinuities in acoustical 
duct systems) where V=2, equations (14) and 
(17a) are simplifield to 

A24+B?=1 and a—f=tn2/2 (19) 
(leaving out the various subscripts 2 in brackets). 
These formulae have been derived and used pre- 
viously [3], [4]. 


5. Examples 


Instances of the type of acoustical junctions which 
are highly-symmetrical as defined above are shown 
in Fig. 2. There are two junctions (a) and (b) with 
N=3 and one (c) with N=4. Highly-symmetrical 


Fig. 2. Examples of highly-symmetrical multi-terminal 
acoustical junctions: (a) plane, 120° junction 
with V=3, (b) orthogonal junction with N=3, 
(ec) tetragonal junction with NV =4. 
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acoustical junctions with V>4 do not exist but the 
matrix conditions of eq. (4) can be approximated 
for certain junctions with V>4, which will not be 
elaborated here. The coupling system connecting the 
various ducts is represented by a hollow sphere in 
which the circular ducts are set so that there is com- 
plete symmetry and equivalency between any two 
of the ducts. In Fig. 2a the axes of the three ducts 
are in one plane and meet at the centre at an angle 
of 120°; in Fig. 2b they join as three orthogonal 
co-ordinate axes at the origin; in Fig. 2c the four 
axes are so arranged that points on each axis, equi- 
distant from the centre, form the corners of a tetra- 
hedron. These acoustical junctions are highly-sym- 
metrical if the input and output reference planes on 
the various ducts are chosen so that they are all at 
the same distance from the centre. Other instances 
of highly-symmetrical junctions can be given, e.g. 
where the cross-section of the ducts is not circular 
or where the mechanical coupling system differs from 
a hollow sphere. 

To check the theory of loss-free highly-symmetrical 
junctions given above, measurements of the charac- 
teristics factors were made with the type of junction 
illustrated in Fig. 2a. A drawing of it is shown in 
Fig. 3 and a photograph in Fig. 4. This highly- 
symmetrical junction consists of an input duct with 
a square inner cross-section of 7.6 cm by 7.6 cm and 
brass walls 0.95 cm thick, and two similar output 
ducts branching off to the right and left at 120°. In 
the centre of the junction is a cylindrical brass ob- 
stacle with a diameter of 5.1 cm (2 in.) and the same 


Fig.3. Sketch of plan and elevation of the acoustical 
junction used in the experiment (measurements 
are in cm). 
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height (7.6cm) as that of the duct walls. Holes 
(0.5 cm diameter) are drilled through the brass cy- 
linder to permit passage of the probe tube along the 
duct axes. The cylinder can be removed as required. 

For practical reasons the centre piece of the junc- 
tions was made of hardwood with a wall thickness of 


Fig.4. Photograph of the 
phone used in the 
moved). 


feast 
re 
Stag-wool ri 
T T 7 
| | | 
Fi Del |e Seay ena ete 
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Frequency f —= 


Sketch of non-reflecting terminal. 


Measured reflection factor plotted against 
the frequency of the non-reflecting termi- 
nals, I and II. 
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1.8 cm and was bolted to the input and output ducts 
where particular care was taken to ensure a good 
seal at the various joints. In Fig. 4 the top of the 
wooden part of the junction has been removed. Each 
output duct connected to the wooden part of the 
junction extends to a distance of 80cm from the 


ase 


acoustical junction and probe tube micro- 
experiments (top part of the junction re- 


centre, where it is terminated by a non-reflecting 
cone made of slag-wool of which one was used in 
previous experiments [4] and the dimensions of 
which are given in Fig. 5a. In Fig. 5b the reflection 
factor | P,/P;| (the magnitude of reflected over in- 
cident sound pressure amplitudes) of each of the two 
terminals is plotted against the frequency, which 
shows that in both cases the reflection factor is 
small (about 0.02) and similar in a wide frequency 
range. 


6. Technique and results of measurements 


For any highly-symmetrical junction the sound 
fields in the various output ducts must be the same 
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when a sound wave is generated in the input. This 
property was first checked with the junction describ- 
ed by probing the sound field along the axis of each 
output duct with a probe microphone through the 
holes in the walls, shown in Fig. 3 and 4 [4]. From 
measurements at a number of different frequencies 
it was found that the sound pressure magnitudes of 
the outgoing progressive waves were practically the 
same in both output ducts (measured differences 
<3.5 percent) and that places with equal phase of 
the sound pressure occurred at practically the same 
distance from the centre of the junction (measured 
differences <1.0mm). This result shows that any 
of the two output ducts can be taken for measuring 
the characteristic transmission factor as is required 
for a highly-symmetrical junction. 

The characteristic reflection and transmission fac- 
tors were then measured by probing the sound fields 
along the axis of the input and of one output duct, 
applying the method described elsewhere [3], [4]. 
The positions of the input and output planes of 
reference were arbitrarily chosen to be perpendicular 
to the axis of each duct and at a distance of 7.6 cm 
(width of the square cross-section) from the two 
corners of the corresponding duct. 

The measurements were made at five frequencies, 
400, 800, 1200, 1600 and 2000 c/s, with and without 
the cylindrical obstacle. As the measured junction 
belongs to the case V = 3 the characteristic reflection 
and transmission factors are called 
Ris) = Aig) exp(j%s))_ and Ts) = Bis) exp (7 fis))- 
The measured values of 4,3), Bis) and a3), Aig) are 
given in Table I and are plotted against the dimen- 
sionless factor O/x=2a/A in Fig. 6a and b, where 
a is the width of the input and output ducts. The 
measured junction (with and without the cylindrical 
obstacle) is approximately loss-free because of the 
hard and smooth wall surfaces. Therefore, the meas- 
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Fig.6. (a) Measured magnitude values A,a), Brg) . 
(b) Measured phase values aig), 63) . 
The values are plotted against the dimensionless 
factor O/x=2 a/A for the two cases without and 
with the cylindrical obstacle. 


ured characteristic factors should fulfil (or nearly 
fulfil) the relations of eqs. (14) and (17a) which, 
with VV =3, read 


Aj +2 Be 


al (20) 


and 


3) — Brg) =a + arc co fe : (21) 
The validity of these formulae was checked with the 
measured values given in Table I. At first, the left 
hand side of eq. (20) was computed and, as shown 
in Table I, it is close to the value of unity, proving 
that the junction is nearly loss-free. The last two 
columns in Table I give the measured and computed 


Table I 
| | ee tea | 
/ eee | Aw | Bey} Oa)/* | Bilt | Ag +2 Be \(@%3)- Bay) | (a¢a)— Bray) = 
e/s em . | | measured measured | computed 
Without obstacle: 
400 86.9 0.175 0.33 0.65 — 1.36 | — 0.38 0.96 — 0.98 — 1.00 
800 43.7 0.348 0.33 0.65 — 1,94 — 0.84 0.96 —1.10 — 1.06 
1200 297 | 0.512 0.33 0.64 033 — 1.24 0.93 + 0.91 + 0.92 
1600 PEG O:705.-'|, . 0,36 0.62 — 0.78 — 1.68 0.90 + 0.90 + 0.83 
2000 17.4 O85 { 638 0.62 — 1.32 — 0.10 0.91 — 1.22 — 1.20 
With obstacle: 
400 88.4 0.172 0.35 0.65 — 1.55 — 0.40 0.97 — 1.15 — 1.12 
800 43.1 0.352 0.48 0.58 — 0.20 — 0.87 0.90 + 0.67 + 0.67 
1200 29.2 0.521 0.61 0.52 — 0.67 — 1.28 0.91 + 0.61 + 0.64 
1600 a7 0.700 0.71 0.46 —1.12 — 1,70 0.93 +. 0.58 +. 0.60 
2000. 17.4 0.873 0.84 0.29 — 1.68 — 0.22 0.87 — 1.46 — 1.44 
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values of (03) —A(3))/a; the computed values are 
derived from A.) and Bi) using eq. (21). The 
measured and computed values agree well and de- 
monstrate the practicability of both the magnitude 
and phase relations derived above. 


7. Discussion 


The phase relation of eq. (17) or (17a) reveals 
a property of loss-free highly-symmetrical junctions 
which is significant in related practical and theoreti- 
cal problems. It has been used to derive eq. (18) 
which provides a general relation between the cha- 
racteristic reflection and transmission factors. This 
relation is useful in acoustical measuring technique 
for checking the accuracy obtained if both factors 
have been measured, or for simplifying the experi- 
mental procedure by measuring one factor only and 
computing the other. 

Although there are limitations in applying the 
theory in acoustical cases where V>4, it can be 
shown that no such restrictions apply to 2 N-terminal 
electrical circuits. Such electrical circuits with the 
required symmetry can be designed for any value! 
of NV. 

In previous work on sound propagation in acousti- 
cal four-terminal filters and ducts [4] the character- 
istic factors were shown to be accurately measurable 
and their particular usefulness was demonstrated in 
the filter theory of certain periodic structures which 


' One way of obtaining highly-symmetrical 2 /-ter- 
minal electrical circuits is by suitable combinations of 
a number of equal symmetrical four-terminal circuits. 
For example, this may be done by connecting the invut 
and output transmission lines in Fig. la (or corres- 
ponding generator and loads in Fig. 1b) with one side 
of each of N symmetrical four-terminals, and circuiting 
the remaining sides of all four-terminals in parallel. 
In another instance it mav be done by connecting 
N(N—1)/2 symmetrical four-terminal circuits in such 
a way that one is always circuited directly between 
any of the N pairs of terminals of Fig. la or b. 
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Die 4. Auflage stellt eine Neubearbeitung dieses Standard- 
werkes der Filtertheorie dar. Die bewahrte Kapiteleinteilung 
(I: Spulen, Kondensatoren, Schwingkreise, II: Zweipole, III: 
Grund-Siebketten, IV: Siebketten mit Grundgliedern und 
M-Gliedern, V: Bandpa8-Siebketten mit Hauptgliedern und 
Nebengliedern, VI: Einflu8 der Verluste) und ihr Inhalt 
wurde im wesentlichen beibehalten. Neu ist ein Abschnitt 
iiber Bandpa8&-Siebketten mit Hauptgliedern und Nebenglie- 
dern. Er enthalt anschauliche Ableitungen dieser Schaltun- 
gen mit einer Reihe von Abwandlungen, die viele verschie- 
dene Schaltungsméglichkeiten erlauben. Unter den Beispie- 
len sind besonders wirtschaftliche Schaltungen besprochen, 
die im Gegensatz zu normalen Bandpafketten weniger Spu- 


BUCHBESPRECHUNG 


ACUSTICA 
Vol. 7 (1957) 


include duct portions connected in cascade at the 
input and output of each section [5], [2]. Following 
the argument of these former results the phase re- 
lation (eq. (17 a)) and eq. (18) will also be helpful 
in the filter theory of similar n-sectional (single- or 
multi-connected) periodic structures of the multi- 
terminal junctions (or analogous electrical circuits) _ 
considered here. 

Loss-free multi-terminal junctions or electrical 
circuits with the highest possible degree of symmetry 
as spezified above are the exception rather than the 
rule in practical applications. However, the present 
method and results are of significance in the funda- 
mentals of acoustics and applicable in similar cases 
of one-dimensional wave propagation in other fields 
of physics. 


8. Conclusions 


Magnitude and phase relations between the cha- 
racteristic reflection and transmission factors of 
highly-symmetrical loss-free multi-terminal 
acoustical junctions have been derived. These re- 
lations constitute an extension to those given pre- 
viously [3] for symmetrical and loss-free discon- 
tinuities in acoustical ducts, and can be applied in 
other fields of analogous cases of one-dimensional 
wave propagation, e.g. on electrical transmission 


lines and in multi-terminal circuits. 
(Received June 15th, 1956.) 
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len als Kondensatoren enthalten. — Um Platz zu gewinnen, 
muften andere Stellen entsprechend gekiirzt werden. 

Eine weitere Neuerung besteht in der zusatzlichen Einftih- 
rung einer von T. Laurent vorgeschlagenen speziellen log- 
arithmischen Frequenzskala, die die Konstruktion von Damp- 
fungs- und Phasenwinkelkurven mit Hilfe von Schablonen 
erlaubt. Auf die Darstellung des Caverschen Verfahrens zur 
Berechnung von Siebschaltungen auf Grund einer vorgegebe- 
nen Betriebsdimpfungskurve ohne Benutzung der Wellen- 
widerstandsanpassung wurde bewuBt verzichtet. 

Die Erginzungen bedeuten also einen Ausbau der ein- 
fachen Filtertechnik, die auch von nicht besonders geschulten 
Mathematikern oder Filtertheoretikern angewendet werden 
kann. Sie kommen besonders ihrer praktischen Anwendung 
zu gute. Der Wert des Buches hat damit nicht nur als Lehr- 
buch fiir den Studenten zugenommen, sondern auch als An- 
leitung zur Konstruktion von einfachen Filtern fiir den aus- 
ubenden Physiker und Techniker. K. Tamm 
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KINEMATOGRAPHISCHE UNTERSUCHUNGEN 
DER SCHWINGUNGSKAVITATION 


von E. Munpry und W. Giru 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Mit einem bei 2,5 kHz angeregten Nickelschwinger wurden Kayvitationsblasen erzeugt 
und mit Hilfe eines Funkenkinematographen bei Bildfolgefrequenzen zwischen 12 500 und 
105 000s“! photographiert. In gashaltigem Wasser wurden bei nicht zu starker Schall- 
energie schwingende Gasblasen beobachtet (,,Gasblasenkavitation“), deren Schwingungs- 
formen als Uberlagerungen der nichtlinearen erzwungenen Blasenschwingung mit einer 
rasch abklingenden Blaseneigenschwingung gedeutet werden konnten. Bei starker Anregung 
wurde sowohl in gashaltigem als auch in entgastem Wasser an Stelle der schwingenden 
Gasblasen eine gro8e Zahl yon kleinen Dampfblaschen gefunden, die nach jeder Anregungs- 
periode wieder nahezu vollstandig zerfallen. Die bei dieser ,,Dampfblasenkavitation“ ent- 
stehenden starken Stofwellen wurden schlierenoptisch photographiert. 


Summary 


Cavitation bubbles are produced from a nickel transducer of frequency 2.5 kc/s and 
photographed at 12500 to 105 000 frames/sec. In aerated water and moderate sound fields, 
bubbles are observed to vibrate in a mode indicating a combination of non-linear forced 
vibration and sudden excitation at the natural frequency. In stronger fields both in aerated 
and degassed water a large number of tiny bubbles, which almost disappear during the 
time-period of one oscillation, are seen. The strong shock waves accompanying this pheno- 
menon are photographed in a schlieren system. 


Sommaire 


On a produit des bulles de cavitation au moyen d’un oscillateur a nickel fonctionnant 
a 2,5 kHz, et on a cinématographié par étincelles ces bulles, la fréquence de succession 
des images étant 12500 et 105 000 images par seconde. Dans le cas de l’eau contenant 
du gaz, on a observé, si l’énergie sonore n’est pas trop grande, des bulles gazeuses oscil- 
lantes («cavitation par bulles gazeuses») dont la forme d’oscillation peut étre considérée 
comme le résultat de la superposition d’une oscillation forcée non linéaire et d’une oscil- 
lation propre s’affaiblissant rapidement. Lorsque l’excitation est intense, on constate, que 
l'eau contienne du gaz ou soit dégazée, un grand nombre de petites bulles de vapeur au 
lieu de bulles de gaz oscillantes; aprés chaque période d’excitation, ces bulles s’écrasent 
presque complétement. On a photographié par striation optique les ondes de choc intenses 
qui se forment dans cette «cavitation par bulles de vapeur». 


1. Einleitung 


Eine haufige und meistens unerwiinschte Folge- 
erscheinung hoher Zugbeanspruchung von Fliissig- 
<eiten ist die Bildung dampf- oder gasgefiillter 
Blasen, die bei Nachlassen der Zugspannung je nach 
ihrem Gasgehalt mehr oder weniger weit wieder zu- 
sammenschlagen. Dieses Phanomen bezeichnet man 
als Kavitation. Die Entstehung solcher Kavitations- 
blasen findet an den Stellen statt, wo sich in Form 
<leiner ungeléster Gasblasen oder nicht vollkommen 
benetzter Festkdérperteilchen Keime befinden, wo also 
sine definierte Grenzflache zwischen Fliissigkeit und 
Gas bereits vor Einwirkung der Zugspannung be- 
steht. Je machdem, wie die erforderliche Zugspan- 
nung hervorgerufen worden ist, unterscheidet man 
verschiedene Kavitationsarten; die Strémungskavita- 


tion beispielsweise wird tiberall da beobachtet, wo 
durch Stromlinienverengung in Fliissigkeiten ein ent- 
sprechender Unterdruck (Bernoullischer Druck) ent- 
steht. Die Schwingungskavitation hingegen ist an die 
Existenz periodischer Druckschwankungen hoher In- 
tensitat gebunden. 

Die Kavitation und ihre Begleiterscheinungen, wie 
Zerst6rung der Oberflache von festen Kérpern [1]. 
chemische und physiko-chemische Prozesse [2], Lu- 
mineszenz [3], [4] und Schallabstrahlung sowie Ver- 
zerrung eines urspringlich vorhandenen Schallfeldes 
[5]... [8] sind bisher Gegenstand zahlreicher Unter- 
suchungen gewesen. Nicht selten wurde dabei ver- 
sucht, auf dem Umweg iiber diese Begleiterscheinun- 
gen Riickschliisse auf das Wesen der Kavitation selbst 
zu ziehen. Von den Arbeiten, die sich mit einer 
direkten optischen Beobachtung weniger oder einzel- 
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ner Kavitationsblasen befassen [9]... [12], be- 
ziehen sich nur zwei auf die Schwingungskavitation. 
So haben Kornretp und Suvorov [13] durch die An- 
fertigung einer groBen Zahl von Momentphotogra- 
phien einzelner Blasen AufschluB iiber den Kavita- 
tionsablauf zu gewinnen versucht. Neuerdings liegt 
eine Arbeit von E1iis [14] vor, in der kinemato- 
graphische Bildreihen von Kavitationsblasen gezeigt 
werden, die mit Hilfe eines in radialer Eigenresonanz 
angeregten Bariumtitanatzylinders erzeugt wurden. 
Jedoch sind die erzielten Ergebnisse nicht zu einer 
Deutung des Phanomens herangezogen worden. 

Der Sinn der vorliegenden Arbeit soll es sein, mit 
Hilfe kinematographischer Untersuchungen einen 
Beitrag zur Beschreibung und Erklarung der Schwin- 
gungskavitation zu liefern. 


2. Die Rayleighsche Beschreibung des Blasen- 
zerfalls 


Eine theoretische Behandlung des Problems haben 
vor allem Lord Rayxereu [15], Notrinck und Nep- 
prras [16], Presser [10] und Brake [19] vorgenom- 
men. An dieser Stelle soll lediglich die Beschreibung 
des titberaus anschaulichen Gedankenganges von Lord 
Ray.rieu das Verstandnis fiir das physikalische Ge- 
schehen bei der Kavitation erleichtern: 

Man stelle sich in einem unendlich ausgedehnten. 
homogenen, inkompressiblen fliissigen Medium einen 
kugelformigen Hohlraum vom Radius R vor. Inner- 
halb des Hohlraumes herrscht der Druck Null, auBer- 
halb in groSer Entfernung soll der konstante Druck 
P herrschen. Infolge dieses AufSendruckes wird sich 
die Fliissigkeit auf das Zentrum des Hohlraumes zu 
bewegen, wobei die genau in radialer Richtung lie- 
gende Geschwindigkeit der Flissigkeitsteilchen um 
so groBer sein wird, je naher sie sich am Zentrum 
des Hohlraumes befinden. Das ist eine direkte Folge 
der Inkompressibilitat. Der Hohlraum wird so mit 
wachsender Geschwindigkeit kleiner, bis schlieflich 
die aus allen Richtungen zusammenstromende Fliis- 
sigkeit heftig aufeinanderprallt. Dabei wiirde an der 
Stelle des Zusammenpralls ein unendlich hoher Druck 
entstehen. Die Voraussetzung der Inkompressibilitat 
ist jedoch jetzt nicht mehr erlaubt, vielmehr bewirkt 
die Kompression der Fliissigkeit, daB es zu einem 
endlichen Spitzendruck kommt, der zum Zentrum 
einer StoBwelle wird, welche sich mit Schallgeschwin- 
digkeit kugelf6rmig ausbreitet. 

Aus einer einfachen energetischen Betrachtung er- 
gibt sich die Zeitabhangigkeit des Blasenzerfalls zu 
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Daraus Jat sich die Zerfallskurve der Blase nume- 
risch ausrechnen. Die gesamte Zerfallszeit wird 
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Nimmt man jetzt an, da8 der Hohlraum nicht exakt 


R’! 2 


eee 


3@ 
aR 0.9 Row |/5 pie) 


evakuiert ist, sondern eine geringe Gasmenge ent- | 


halt, so wird sich der Gasdruck bei der Volumen- 


verringerung soweit steigern, bis schlieBlich der auf - 


diese Weise in der Blase entstandene Gegendruck die 
Bewegung der Fliissigkeit zunaichst abbremst und 
dann umkehrt. Es ergibt sich so eine stark nicht- 


lineare Schwingung der Blase um ihr Gleichgewichts- | 


volumen, die jedoch infolge der hohen Dampfung 
schon nach wenigen Perioden zum Stillstand kommt. 
Auff§erdem wird die Hohe des beim Zusammenfallen 
der Blase erreichten Druckes in der abgestrahlten 


StoBwelle durch die Federwirkung des Gases herab- | 


gesetzt. 


3. Beschreibung der Apparatur 


Zur Erzeugung der Kavitationsblasen dient ein — 


1m langes, in der Mitte eingespanntes Nickelrohr. 


Kine das Rohr umgreifende Luftspule mit zwei Wick- | 
lungen sorgt fiir die notige Vormagnetisierung und — 


die Anregung des Schwingers in seiner tiefsten Eigen- 
frequenz bei etwa 2,5kHz. Zur Vermeidung von 
Wirbelstromverlusten ist das Rohr langs aufgeschlitzt. 
In einer zweiten Spule nahe seinem unteren Ende 
wird infolge der Rohrschwingungen eine Wechsel- 
spannung induziert, die zur Kontrolle der Schwin- 
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gung hinsichtlich ihrer Sinusform auf die vertikalen — 


Ablenkplatten eines Kathodenstrahloszillographen — 
gegeben wird. Das mit einem Messingplattchen ab- — 
geschlossene Rohr taucht in eine wassergefiillte Kii- — 
vette, und die Kavitationsblasen erscheinen nun ge- | 
nau unter der Mitte dieses Messingplattchens. Die — 


zur Anregung des Schwingers erforderliche Wechsel- 


spannung wird aus einem Schwebungssummer be- | 


zogen und gelangt tiber einen 70-W-Leistungsverstar- 
ker und einen Abstimmkondensator an die Erreger- 


spule. Diese Art der Blasenerzeugung hat den Nach-— 


teil, daB die Blasen nicht im freien Medium ent- 
stehen, sondern an der Grenzflache zwischen Schwin- 


ger und Fliissigkeit; sie hat jedoch andererseits den 


fur eine optische Beobachtung erheblichen Vorzug, 
da’ der Ort der Blasenbildung festliegt. AuBerdem 
ist die Tatsache von Vorteil, da nur einzelne oder 
wenige Blasen entstehen und man deshalb den gan- 
zen Vorgang gut iibersehen kann. 

Der wesentliche Teil der optischen Apparatur be- 
steht aus dem von Gin [12] entwickelten und be- 
schriebenen Funkenkinematographen. Ein Dreh- 
spiegel wirft das Bild des Objektes auf einen halb- 
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kreisformig in die Kamera eingezogenen Film; als 
intermittierende Lichtquelle dient eine Funkenstrecke. 
Eine genaue Beschreibung dieses Gerates findet sich 
bei Gin [12] und soll an dieser Stelle nicht wieder- 
holt werden, da der Kinematograph lediglich in tech- 
nischer Hinsicht etwas verbessert wurde. So war bei- 
spielsweise eine Steigerung der Bildfolgefrequenz auf 
105 000 s+ moglich. Eine schematische Darstellung 
der tbrigen optischen Anordnung findet sich in 
Abb. 1. Die als intermittierende Lichtquelle dienende 
Funkenstrecke wird mittels eines Kondensors zur 
Ausleuchtung des Objektes benutzt, das von einem 
Objektiv tiber den Drehspiegel auf den Film ab- 
gebildet wird. Eine mit dem Licht einer Glihlampe 
erzielte zweite Abbildung des Objektes auf einen 
Wandschirm erlaubt die Messung des Luftvolumens 
einer schwingenden Blase auf folgende Weise: 
Wenige Millisekunden nach Belichtung des Filmes 
wird der Schwinger automatisch ausgeschaltet. Die 
in der photographierten Blase enthaltene Luft bleibt 
dann in Form eines runden Blaschens unter dem 
Schwinger zuriick und wird auf dem zweiten opti- 
schen Wege auf den Wandschirm projiziert und dort 
ausgemessen. 


Reuterlampe 


Funkenstrecke 


Drehspiegel 


Kondensor 


Spiegel Wandschirm 


Abb. 1. Skizze der optischen Anordnung. 


Zum Nachweis der beim Zusammenfallen einer 
Blase entstehenden Stofwellen wurden 
optische Aufnahmen angefertigt. Zwei nach dem 
Téplerschen Verfahren aufgebaute Anordnungen 
wurden benutzt, von denen die eine bei etwas ver- 
kleinertem Abbildungsmafstab die groBte erzielbare 
Schlierenempfindlichkeit ergab, waihrend die andere 
bei Verzicht auf hohe Schlierenempfindlichkeit eine 
vergroBerte Abbildung des Objektes ermoéglichte, die 
sich in einigen Fallen als notig erwiesen hat. 


schlieren- 


4. Aufnahmen von Kavitationsblasen 
und Diskussion der Ergebnisse 


4.1. Gashaltiges Wasser, schwache Anregung 

Es soll zunichst eine Art von Blasen besprochen 
werden, die in gashaltigem Wasser (Leitungswasser 
oder destilliertem Wasser) bei nicht zu starker 
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Schwingung des Nickelrohres entstehen. Eine direkte 
mikroskopische Messung der Schwingungsweite des 
Schwingers ergab fiir diese Versuche Werte zwischen 
0,08 und 0,09mm. Man beobachtet eine milchig 
weibe, in ihren Umrissen verschwommene Blase, die 
etwa halbkugelformig direkt am Schwinger haftet. 
Uber einen ziemlich langen Zeitraum hinweg — 
mehr als 10* Schwingungsperioden — wachst die 
Blase zunachst an, bis sie bei einer bestimmten Groke 
anscheinend heftiger schwingt als vorher und sich 
dabei mit einem dichten Nebel sehr kleiner Blaschen 
umgibt; danach wird die Schwingungsamplitude der 
Blase wieder erheblich kleiner, sie beginnt unter dem 
Schwinger hin und her zu wandern, ohne jetzt ihre 
GroBe noch merklich zu andern, und schlieBlich ent- 
weicht sie seitlich und steigt an die Oberflache auf. 


Abb. 2. Zeitaufnahmen schwingender Luftblasen; Ver- 


groBerung 2-fach. 
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In Abb. 2 werden solche Blasen gezeigt, wie sie dem 
Beobachter mit blo&em Auge erscheinen; zwar han- 
delt es sich bei jeder einzelnen Aufnahme um eine 
andere Blase, die Bildreihe gibt aber eine gute Dar- 
stellung der Lebensgeschichte einer einzigen Blase in 
groben Ziigen. Man erkennt die verschiedenen Sta- 
dien und sieht auf dem zweiten Bild von unten deut- 
lich den Nebel kleiner Blaschen, die infolge der 
langen Belichtungszeit als Striche erscheinen. Die 
Wirbelbildung in der Umgebung der Blase soll 
hier nicht besprochen werden. Die Frage nach der 
Ursache fiir das allmahliche Anwachsen der Blasen 
wird an einer anderen Stelle unter Zuhilfenahme 
weiterer Ergebnisse erortert. 

Der nachste Schritt beim Studium der schwingen- 
den Luftblasen besteht in der Anfertigung von 
Momentaufnahmen mit dem Licht eines einzelnen 


Abb. 3. Momentaufnahmen schwingender Luftblasen; 
VergréBerung etwa 9-fach. 
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Funkeniiberschlages bei, Dunkelfeldbeleuchtung; sie 
ergeben eine scharfe Abbildung, da die Belichtungs- 
zeit klein ist gegen die Schwingungsdauer. Abb. 3 
zeigt zwei solche Blasenbilder, davon eines wieder- 
um mit dem Nebel kleiner Blaschen. Die bereits von 
KornreLp und Suvorov [13] beobachtete stark zer- 
kliiftete Blasenoberflache ist deutlich zu erkennen. 
Wahrscheinlich entstehen also die vielen kleinen 
Blaschen durch Abschniirung von der grofen Blase, 
doch soll auch diese Frage nicht naher erdrtert wer- 
den. 

Um Aufschluftiber die Art der Volumenschwin- 
gung der Blasen zu erhalten, wurden nach den bis- 
her geschilderten Vorversuchen kinematographische 
Bildreihen von schwingenden Blasen angefertigt; 
Abb. 4 zeigt eine solche Bildreihe. Die Blase er- 
scheint hier in der Art eines Schattenrisses als bizarr 
geformtes dunkles Gebilde vor einem hellen Hinter- 
grund. Die Beobachtungsrichtung ist horizontal, die 
obere Bildkante stellt also das mit dem ‘Messing- 
plattchen verschlossene untere Ende des Nickel- 
schwingers dar. 

Wenn man auf jeder Aufnahme einer solchen Bild- 
reihe das Blasenbild ausplanimetriert und aus der so 
erhaltenen Flache F nach der Beziehung R ~~ / F 
einen mittleren Blasenradius bestimmt, so laBt sich 
dieser als Funktion der Zeit auftragen, da man aus 
der gemessenen Umdrehungsgeschwindigkeit des 
Drehspiegels den genauen zeitlichen Abstand zwi- 
schen zwei Aufnahmen feststellen kann. Bei geeigne- 
ter Eichung des AbbildungsmaBstabes ist der mitt- 
lere Blasenradius genau der Radius einer kugelf6rmi- 
gen Blase gleichen Volumens, wenn man voraussetzt, 
daB die Kavitationsblase in bezug auf die Langs- 
achse des Schwingers annahernd symmetrisch ist. 
Diese Voraussetzung ist sicher nicht gut erfiullt, die 
Ergebnisse zeigen jedoch, dah sie offenbar zulassig 
ist. Die in Abb. 5 enthaltene charakteristische Aus- 
wahl derartiger Mefikurven wurde einer Versuchs- 


— 


Abb. 4. Schwingende Luftblase; Bildfolgefrequenz 50000s~1; Vergré8erung etwa 3-fach. 
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Blasenradius R ——= 
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Abb. 5. Radius schwingender Gasblasen als Funktion der Zeit fiir’ verschiedene Blasengré- 
Ben; Erlauterung im Text. 


eihe entnommen, bei der unter genauer Konstant- 
altung der Anregungsenergie (Schwingungsweite 
les Nickelrohres 0,082 + 0,002 mm) zahlreiche Bla- 
en verschiedener Grofe photographiert wurden. Zu 
eder Aufnahme wurde nach der eingangs beschrie- 
yenen Methode der Radius Ry der Blase in ruhendem 
‘ustand bestimmt. Er ist als strichpunktierte Gerade 
n den Schwingungskurven der Abb.5 enthalten. 
‘olgende Ergebnisse lassen sich den Kurven ent- 
ehmen: 


Alle beobachteten Vorginge sind bei einer Mef- 
genauigkeit von +1,5% periodisch in 2.5 kHz. Je- 
doch ist die Schwingungsform von der BlasengréBe 
entscheidend abhangig. Wahrend namlich die kleinste 
Blase (Ruheradius Ry =0,3 mm) deutlich neben dem 
in 2,5 kHz periodischen tiefsten Minimum zwei wei- 
tere Minima erkennen lat, finden sich bei der groB- 
ten Blase (Ry = 1,42 mm) keine zusitzlichen Minima 
mehr; tiber die ferner erkennbare Abnahme der Am- 
plitude mit wachsender BlasengréSe lassen die vor- 
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liegenden Ergebnisse keine Schliisse zu. Die verzerrte 
Sinusform selbst bei den Blasen, deren Schwingun- 
gen nur die Anregungsfrequenz deutlich zeigen, war 
zu erwarten; denn da Gasblasen in Fliissigkeiten 
nichtlineare Schwingungsgebilde sind, da sowohl die 
Federsteife wie die mitschwingende Mediummasse 
veranderlich sind, ist seit langem bekannt. Diese bei- 
den Ursachen einer Nichtlinearitat wirken bei der 
schwingenden Luftblase in einer Richtung und fiihren 
zu einer Verbreiterung der Schwingungsmaxima, 
wihrend die Minima schmaler werden. Eine solche 
Kurvenform ist in allen Beispielen der Abb. 5 gut 
zu erkennen. In welcher Weise diese Nichtlinearitat 
an dem bereits von Brake [19] ausfuhrlich behan- 
delten allmahlichen Anwachsen der mittleren Blasen- 
groBe uber viele Schwingungsperioden hinweg be- 
teiligt ist, lat sich an Hand der Abb. 5 einfach ver- 
folgen. Man erkennt namlich, da die nichtlinear 
schwingenden Blasen sich im Zustand des Unter- 
druckes langer aufhalten als in dem des Uberdruckes, 
mit anderen Worten, da der zeitliche Mittelwert der 
BlasengroBe oberhalb des Ruhewertes Ry, also in 
der Unterdruckphase liegt. Diese Tatsache bewirkt 
eine bevorzugte Diffusion von Luft aus der Flissig- 
keit in das Blaseninnere hinein, also ein allmahliches 
Anwachsen, wie es bereits oben beschrieben wurde. 
Mit abnehmender Amplitude muf das Blasenwachs- 
tum schwacher werden, was den experimentellen Be- 
obachtungen auch entspricht. 


Um die bereits erwahnten teilweise festgestellten 
zusatzlichen Maxima und Minima der Schwingungs- 
kurven zu verstehen, mufi man sich vergegenwarti- 
gen, daB die Schwingung einer Luftblase mit wach- 
sender Amplitude in ihrem abfallenden Teil der 
Ray.etcuschen Zerfallskurve eines kugelformigen 
Vakuums immer ahnlicher wird, solange die in der 
Blase eingeschlossene Luft noch nicht zu stark kom- 
primiert ist. In der Tat zeigen die gemessenen Kur- 
ven in ihrem Hauptabfall zum Teil eine starke Ahn- 
lichkeit mit den jeweils fiir die entsprechende Blasen- 
groBe berechneten Rayercuschen Zerfallskurven, die 
als Strichlinien in Abb. 5 eingezeichnet worden sind. 
Die groBe Steilheit der Schwingungskurve kurz vor 
Erreichen ihres Minimums besagt, da der Druck 
der in der Blase eingeschlossenen Luft wesentlich 
schneller ansteigt, als es der sinusférmigen auferen 
Druckanregung entspricht. Das aber bewirkt eine 
Anregung der freien Eigenschwingung der Blase 
ahnlich einer StoBanregung, derzufolge jedes schwin- 
gungsfahige Gebilde in seiner Eigenfrequenz aus- 
schwingt. Die in Abb.5 vorliegenden Kurven sind 
demnach zu verstehen als eine Uberlagerung aus der 
nichtlinearen erzwungenen Schwingung der Gasblase 
mit der Frequenz des Nickelschwingers und dem 
durch den Druckstof bei dem sehr steilen Abfall der 
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BlasengroBe bedingten Ausschwingen der Blase in 
ihrer Eigenfrequenz. Eine experimentelle Bestati- 
gung dieser Vorstellung wird durch die Tatsache ge- 
liefert, daB die aus dem gemessenen Ruheradius Ry 
berechnete Periodenlange t) der Blaseneigenschwin- 
gung gut ubereinstimmt mit dem zeitlichen Abstand 
zwischen dem Hauptminimum und dem ersten zu- 
satzlichen Minimum bzw. zwischen zwei Zusatz- 
minima. tg ist in die Abb. 5 eingezeichnet. AuBerdem 
ist unter Beriicksichtigung der Gleichung 


1-3/3 x Py 
aR = 326 Hz-em, (3) 


fo Ro =—- 


die zuerst von Minnarrr [17] aufgestellt und von 
Meyer und Tamm [18] experimentell bestatigt wor- 
den ist, die berechnete Abhangigkeit zwischen ft) und 
Ry in Abb. 6 als Gerade dargestellt. Die gemessenen 
Zeitdifferenzen zwischen zwei Minima als Funktion 
von Ry enthalt die Abb. 6 als MeBpunkte. Die ver- 
haltnismaBig gute Ubereinstimmung der Ergebnisse 
mit der theoretischen Kurve weist ebenfalls darauf 
hin, da diese Interpretation der Blasenschwingun- 
gen offenbar zutrifft. In Abb. 7 erkennt man etwa in 
der Mitte zwischen den beiden eine Anregungs- 
periode einschlieSenden Hauptminima ein solches 
zusatzliches Minimum der Blasengrofe. 


3 
3 


oO 
a 


(oe) 
> 


Oo 
ht 


Radius Ro der ruhenden Blasen —=—- 


0 40. = 80 120 160 ps 200 
Periodenlange ty der Blaseneigenschwingung —= 


Abb. 6. Periodenlange der Blaseneigenschwingung als 
Funktion des Ruheradius im Vergleich zur 
Theorie. 


4.2. Gashaltiges Wasser, starke Anregung 


Ein ganz anderes Erscheinungsbild erhalt man, 
wenn man — ebenfalls in gashaltigem Wasser — 
die Amplitude des Schwingers erheblich steigert. Bei 
der maximal moglichen Schwingungsweite von 
0.35 mm kommt es nicht zur Bildung einzelner gro- 
Berer Blasen, sondern statt dessen entsteht eine 
Wolke vieler kleiner Blasen, und ein Strahl winziger 
Luftblaschen wird infolge des Schallstrahlungs- 
druckes (,,Quarzwind“) in das Wasser abgestoBen. 
Abb. 8 stellt den Vorgang so dar, wie er mit bloBem 
Auge wahrgenommen wird. Abb. 9 zeigt einen Aus- 
schnitt aus einer Blasenwolke, in dem die einzelnen 
Blasen gut zu sehen sind. Man erkennt bei dieser 
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Abb. 7. 
3,5-fach. 


Abb. 8. Zeitaufnahme der Kavitation unter dem Schwin- 
ger bei starker Anregung; VergroBerung 2- 
fach. 


kinematographischen Bildreihe, da} die Blasen voll- 
kommen zusammenschlagen; es muf deshalb ein Fall 
von Kavitation vorliegen, bei dem die entstehenden 
Blasen fast nur Dampf enthalten und daher in der 
Uberdruckphase nahezu véllig verschwinden kénnen, 
mindestens bis zu einer Grofe, die mit dem benutz- 
ten Aufnahmeverfahren nicht mehr sichtbar gemacht 


Schwingende Gasblase wahrend einer vollen Periode; Bildfolgefrequenz 105 000s 


1; VergroBerung 


werden konnte. Aus mikroskopisch kleinen Keimen 
bildet sich eine Fille von Blaschen, die noch wah- 
rend derselben Anregungsperiode wieder ganz zu- 
sammenschlagen, wobei gelegentlich ein einmaliges 
schwaches Nachschwingen beobachtet wird. Wenn der 
Unterdruck der nachfolgenden Anregungsperiode 
wieder weit genug angewachsen ist, wiederholt sich 
der Vorgang. Da Anzahl und Verteilung der einzel- 
nen Blaschen sowie ihre Grobe von Periode zu Pe- 
riode sehr verschieden sind, kann man annehmen, 
daB zwischen den Kavitationserscheinungen in den 
einzelnen Schwingungsperioden kein Zusammenhang 
besteht, abgesehen von der Tatsache, da! moglicher- 
weise ein zusammenschlagendes Blaschen einen oder 
mehrere Keime zurticklaSt, die die Bildung neuer 
Blasen in der folgenden Periode begiinstigen. Die 
GroBe eines einzelnen Blaschens als Funktion der 


Zeit — Abb. 10 


eewonnen wurde, zeigt zwar qualitativ in ihrem Zer- 


~, die wieder durch Planimetrieren 


Abb. 9. Kinematographische Bildreihe von Kavitationsblasen bei starker Schwingeranregung; 


50.900s~+: VergréBerung etwa 7-fach. 


Bildfolgefrequenz 
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Abb. 10. Zerfallskurven einzelner Dampfblaschen. 
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fallsabschnitt Ahnlichkeit mit der Rayzereuschen 
Kurve, die aus der Blasengroe berechnete Zerfalls- 
zeit ist jedoch in allen Fallen kleiner als die gemes- 
sene Zerfallszeit. Dieser Vergleich mit der Theorie 
ist allerdings auch nicht mehr zulassig, weil bei der 
starken Anregung die in der Raytereuschen Rech- 
nung gemachte Voraussetzung eines bekannten, festen 
AuBendruckes nicht mehr erfullt ist; denn durch die 
starke Amplitude und die dadurch bedingte heftige 
Gleichstromung werden die Druckverhdltnisse am 
Ort der Blasenentstehung so uniibersichtlich, da8 sich 
keine bestimmten Aussagen mehr machen lassen. 
Aufierdem befinden sich mehrere Blasen dicht bei- 
einander und beeinflussen sich gegenseitig in ihrem 
Verhalten. 

Bei hinreichender Keimzahl und genitigendem 
Blasenwachstum konnen sich unter Umstanden meh- 
rere solcher Blasen zu einem groferen Dampfraum 
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vereinigen, was bei den vorliegenden Versuchen stets 
mit den Blasen in allernachster Nahe des Schwingers 
— Blasen, die die Messingflache beriihren — ge- 
schehen ist. Unter Beriicksichtigung der hydrodyna- 
mischen Verhaltnisse vor dem stark schwingenden 
Nickelrohr kann man sich vorstellen, da sich auf 
diese Weise wahrend einer einzigen Schwingungs- 
periode ein Dampfvolumen bilden kann, dessen Zer- 
fallszeit die Anregungsperiode tiberdauert. Offenbar 


ist KornreLp und Suvorov [13], die in ihren Ver- — 


suchen mit erheblich hoheren akustischen Leistungen 
gearbeitet haben, die. Erzeugung solcher Blasen ge- 
lungen. ; 


4.3. Entgastes Wasser bei starker Anregung 


Ahnliche Versuche wie die oben beschriebenen 
wurden nun in entgastem Wasser gemacht; das Was- 
ser wurde vorher durch eine enge Diise hindurch in 
eine evakuierte Kochflasche gesogen, wobei ihm ein 
erheblicher Teil der gelosten Gase entzogen werden 
konnte. 
Wasser bei maximaler Amplitude des Schwingers die 
Anregungsenergie gerade ausreichend zur Kavita- 
tionsbildung. Denn es verging nach Einschalten des 
Nickelrohres eine Zeit zwischen etwa 5 und 30 Se- 


Abb. 11. Zeitaufnahme der Kavitation in entgastem 
Wasser bei starker Anregung; Vergroferung 
2-fach. 


Abb. 12. Entstehung und Zerfall von Dampfblasen bei Kavitation in entgastem Wasser; Bildfolgefrequenz 


55 000s~!; VergréS8erung 7-fach. 


Offenbar war in dem so vorbehandelten . 
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kunden, wahrend derer keinerlei Kavitation beob- 
achtet wurde — bis sich in irgend einem Augenblick 
plotzlich mit knackendem Einsetzen eines deutlichen 
Zischens eine Blasenwolke unter dem Schwinger ahn- 
lich der friiher beschriebenen ausbildete, die nun bis 
zum Abschalten der Apparatur erhalten blieb. Diese 
Beobachtung unterstiitzt die oben geauBerte Ver- 
mutung, dai beim Zusammenschlagen einer Blase 
Keime zuriickbleiben, die die Entstehung neuer Bla- 
sen in der nachsten Periode fordern. Eine Zeit- 
aufnahme — Abb.11 — zeigt die Blasenwolke in 
diesem Fall. Der Strahl feiner Luftblaschen fehllt, 
was bei der Benutzung entgasten Wassers nicht ver- 
wundert. Kinematographische Aufnahmen dieses Ka- 
vitationsvorganges zeigen einen sehr ahnlichen Ver- 
lauf des Prozesses wie bei gashaltigem Wasser 
(Abb. 12). Die Unterschiede, die wahrscheinlich im 
etwas heftigeren Zusammenschlagen und in erschwer- 
ter Blasenbildung bei entgastem Wasser im Gegen- 
satz zu gashaltigem bestehen, sind zu gering, um mit 
der benutzten Apparatur deutlich hervorgehoben 
werden zu kénnen. 


Die Ahnlichkeit der Kavitationserscheinungen in 
gashaltigem und entgastem Wasser bei geniigend 
starker Anregung weist darauf hin, da das Verhal- 
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ten von Kavitationsblasen — ,,Gasblasenkavitation“ 
oder ,,Dampfblasenkavitation“® — nicht allein vom 


Gasgehalt der Fliissigkeit, sondern auch erheblich 
von der Anregungsenergie bestimmt wird. 


5. DruckstoBe und Materialzerstorung 
bei der Kavitation 


Mit der bereits erwahnten schlierenoptischen Me- 
thode wurde nun versucht, die bei der Kavitation 
entstehenden Druckstéfe sichtbar zu machen. Abb. 13 
enthalt derartige Aufnahmen; eine groBe Zahl von 
StoBwellen ist zu erkennen, und zwar findet man auf 
diesen wie auf allen anderen Aufnahmen, dai das 
Zentrum der Stowellen immer in unmittelbarer 
Nahe der Schwingeroberflache liegt. Das ist damit 
zu erklaren, dafi sich nur hier Blasen bilden, die groB 
genug sind, um geniigend heftig zusammenschlagen 
zu konnen. Da die zerfallenden Blasen oder Teil- 
blasen — beim Zerfall spaltet sich oft eine Blase in 
mehrere Teile auf — nicht alle gleich grof} sind, 
schlagen sie nitht gleichzeitig, sondern nacheinander 
zusammen. In Abb. 14 erkennt man auf drei auf- 
einander folgenden Aufnahmen StoBwellen. Aus der 
Bildfolgefrequenz ergibt sich, da sie wahrend eines 
Zeitraumes von mindestens 80 Us entstanden sind; 


* q Soh: 


Abb. 13. Schlierenoptische Aufnahmen der KavitationsdruckstoBe; Bildfolgefrequenz 12 500s~'; 


2.5-fach. 


VergroBerung 


Abb. 14. Schlierenoptische Aufnahmen der Kavitationsdruckstéfe; Bildfolgefrequenz 35 000s~'; 
2,7-fach. 


VergroBerung 
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das entspricht 1/5 der Anregungsperiode. Diese Tat- 
sache, da8 namlich die StoBwellen nicht streng perio- 
disch, sondern in gewissen Grenzen statistisch regel- 
los entstehen, ist ein Grund fiir das Vorhandensein 
eines Rauschpegels im Spektrum des bei Kavitation 
abgestrahlten Schalls, wie er von Escuz [7] gemes- 
sen wurde. Bei schlierenoptischen Aufnahmen von 
schwingenden Gasblasen, wie sie oben beschrieben 
wurden, konnten keinerlei Stofwellen festgestellt 
werden, obwohl auch solche Blasen einen Kavita- 
tionsangriff am Schwingermaterial hervorriefen. Die- 
ser war allerdings so schwach, da} er erst nach mehr- 
stiindiger Einwirkungszeit der Blase am Messing- 
plattchen wahrzunehmen war. wohingegen eine Ma- 
terialschadigung etwa gleichen Ausmafes bei Dampf- 
blasenkavitation schon nach einer Einwirkungszeit 
von wenigen Sekunden eingetreten war. 

Auch die schlierenoptischen Versuche lieBen kei- 
nen deutlichen Unterschied zwischen Dampfblasen- 
kavitation in gashaltigem und in entgastem Wasser 
hervortreten, obwohl eine StoBwelle um so schmaler 
und hoher sein wird, je weniger Gas die zusammen- 
fallende Blase enthalt. Alle beobachteten Stofwellen 
haben dieselbe Form wie die von Gtiru [12] photo- 
graphierten. In Abb. 15 tritt diese Form besonders 
deutlich hervor; die flache Flanke der Welle befindet 
sich auf ihrer Vorderseite, die steile Flanke auf der 
Riickseite. 

Als Ergebnis dieser letzten Versuchsreihe wird 
festgehalten, daB man in jedem Falle der Kavitation 
einen Materialschaden zu erwarten hat, wenn er 
auch bei der Gasblasenkavitation sehr viel schwacher 
ist als bei der Dampfblasenkavitation. StoBwellen 
als Folge des Blasenzerfalls konnten jedoch nur bei 
der Dampfblasenkavitation festgestellt werden. 

Herrn Professor Dr. E. Meyer sind wir fir die 
Anregung zu dieser Arbeit sowie fiir viele wertvolle 
Ratschlage bei ihrer Durchfithrung zu, besonderem 
Dank verpflichtet. AuSerdem haben wir Herrn Dr. 


Abb. 15. Form einer StoBwelle bei Dampfblasenkavita- 
tion. 
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Kurrze ftir weitere Hinweise und der Deutschen 
Forschungsgemeinschaft fiir ihre finanzielle Unter- 
stiitzung der Arbeit zu danken. 


(Eingegangen am 1. Juli 1956.) 
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EXPERIMENTAL EVIDENCE OF NONLINEAR RELATION 
Pee IN SPECTRAL NUMBER AND ULTRASONIC 


DNDENSITY 
by K. Samau 


Physics Department, Ravenshaw College, Cuttack, Utkal University, India 


Summary 


The relation between ultrasonic intensity and the number of diffraction spectra has been 
shown to be nonlinear for smaller numbers which contradicts SanpErs’ suggestion. With 
the application of the relation, the ultrasonic absorption coefficients of CgH,, CS), and 
CCl, are measured. The expected values were obtained. 


Sommaire 


On a montré que la relation entre l’intensité des ultrasons et l’ordre du spectre de 
diffraction n’est pas linéaire pour les petits ordres, ce qui est en contradiction avec une 
suggestion de Sanpers. Les coefficients d’absorption des ultrasons dans CyH,, CS, et 
CCl, ont été mesurés par application de cette relation. On a trouvé les valeurs prévues. 


Zusammenfassung 


Es ist gezeigt worden, daB entgegen den Vorstellungen von Sanprrs die Beziehung 
zwischen Ultraschallintensitat und Ordnungszahl des Beugungsspektrums fiir kleine Ord- 
nungszahlen nicht linear ist. Unter Anwendung dieser Beziehung werden die Absorptions- 


koeffizienten von Cy,H,, CS; und CCl, gemessen. Man erhilt die erwarteten Werte. 


1. Introduction 


Raman and Nacenpra Natn’s theory [1] on in- 
tensity variation of spectral orders obtained by ul- 
trasonic waves, has been corroborated by the ex- 
periment of Sanpers [2]. Sanpers also considered 
that the observed spectral number n is directly pro- 
portional to v or the square root of ultrasonic in- 
tensity. For confirmation of this idea he had ex- 
tended its application to the experimental result of 
Wyss [3], but the result of the comparison was not 
satisfactory. The experiment performed by Samat 
and others [4] suggested that on reorientation of 
Raman and Natn’s theory a nonlinear relation be- 
tween the ultrasonic intensity or the refractive index 
variation and the observed spectral number would 
be possible. From a critical analysis of Raman’s 
theoretical relations one obtains a graph of n and v 


Fig. 1, A graph between n and v where n is the spec- 
tral number and v=224Aul/2 (Au refractive 
index variation). 


(Fig. 1). The detailed working out of the relation- 
ship is being published elsewhere [5]. 

To test the hypothesis an attempt has been made 
to determine ultrasonic absorption coefficients of 
liquids on the basis of a modified n and v relation- 
ship by observing the visible spectral number ai 
different distances from the source of sound. 


2. Experimental arrangement 


A piezoelectric quartz of natural frequency 1 Mc/s 
was driven at 5 Mc/s. The electric generator con- 
sisted of a full-wave rectifier and an oscillating cir- 
cuit of the Harriey type. The power supplied to 
the quartz was maintained constant. The vibrating 
quartz was fixed inside a rectangular glass vessel of 
a length of 25cm, having optically plane faces. 
The maximum number of symmetrical spectra. ob- 
servable, on both the sides of the undeviated ray 
with a normal illumination by the properly colli- 
mated light from a monochromatic sodium lamp, 
was adjusted. The distance between the quartz and 
the path of the undeviated ray was varied by mov- 
ing the vessel with a mechanical arrangement at- 
tached to the table. The liquids were pure and 
supplied by E. Merck. 


3. Experimental results 


By increasing the distance of the quartz from 
the light path the number of spectra decreases 
gradually. Within the distance of 20cm all the 
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spectra except the central one disappeared in ben- 
zene and carbon disulphide. This assures progres- 
sive waves. In the case of carbon tetrachloride the 
first order is still detectable at the remotest distance 
possible within the cell, indicating that the condition 
of progressive waves is vitiated, there being no ab- 
sorber at the back. By moving the vessel slowly the 
positions of the gradual appearance of spectra were 
read from the scale. 


4. Calculation of absorption coefficients 


The sound intensity /) of a plane sound wave 
decreases in a liquid by passage over a distance x 
to the value 


ieee eres (1) 
where a represents the absorption coefficient of 
sound amplitude. As v® is directly proportional to / 
the relation reduces to the form 

V;=vje ~~ (2) 
(2a) 
where ¢ is a constant. If n, be the order appearing 
at the distance x, we can find the corresponding v, 


or logevg;= —ax+ce, 


from the graphical relation between n and v shown 
in Fig. 1. The graph between the natural logarithm 
of v and x will be a straight line whose slope will 
give the absorption coefficient. 


; 


10g 10 x , LOG 10% 


-10 = aa 
"0 4 8 12 16 cm 20 


,o 


Fig.2. Graphs between x and logv and logn for dif- 
ferent liquids, where x is the distance between 
light path and sound source, and v=2 2 Aw l/Z. 


In Fig. 2, the graphs of logyyv—z are plotted 
and the expected straight lines are obtained. The 
slope values are 0.644, 0.088 and 0.046 for CS,, 
CH, and CCl, respectively. The calculated values 
of a//* (f being the frequency) obtained from the 
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Table I 


Observed | Temperature | PINKERTON’s |RICHARDSON’S 
alf?-10%e.2.8 of the experi- | value of value of 
| “| ment in °C |a/f?- 101". g.s.ja/f?- 10" c.g. 8. 


Liquids 


CH, gos |. 30 300800 
CS, | 5920 || 25 | + 7400 enone 
Col 424) (28 533 1 a0” 


slopes are shown in Table I along with PrvKerron’s 
value [6] and that quoted by Ricnarpson [7]. 

If the linear relation between n and v is taken 
as suggested by Sanpers, the eq. (2) will be in the 
form as * r 
(2b) 
and a plot of log,n against x would give straight 
line. In Fig. 2 logy) n, against x is also plotted for 
these three liquids. They are not straight lines ex- 
cept for the case of CCl, in which case, however, 
the slope of the logyyn—2 line is obtained as 
0.029, giving a/f? x 1017 = 268 c.g.s. which is much 
different from the expected value. 


Ny Cee 


. 
5. Remarks on Wyss’ method . 
i 

The relation between n and v is obtained with a 
specific condition that n is the observable spectral 
number when the last order has a constant threshold 
intensity for visibility. In Wyss’ method of obser- 
vation of spectra at different distances, the various 
intensities were estimated and this is not expected 
to conform to the above condition. Moreover the 
vessel used by him was too small to establish pro- 
gressive sound waves. Thus the application of the 
(n—v) relation to Wyss’ observations for the de- 
termination of absorption coefficient is not ver 
well adapted to the experimentally observed values. 


Acknowledgment 


The author is indebted to Professor A. K. Durra, 
D. Se, F.N.1., for his help and guidance durin 


the progress of the work. hx 
(Received 5th August, 1956.) 


References 


[1] Raman, C.V. and Naru, N.S., Proc. Ind. Acad 
A 2 [1936], 406 —412. 

[2] Sanpers, F.H., Canad. J. Res. 14 [1936], 158— 
170. 

[3] Wyss, R., Helv. Phys. Acta 7 [1934], 406—412. 

[4] Samat, K., Parnark, K.M. and Durra, A.K., Ind 
J. Phys. 29 [1955], 548—557, 

[5] Samaz, K., Sent for Publication. 

[6] Pinkerton, J. M. M., Proc. Phys. Soc. B 62 [1949] 
129 —150. 

[7] Ricuarpson, E.G., Ultrasonic physics. Elsevier 
Publ. Co. Amsterdam 1952, pp. 152 —153. 


